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SECTION  I 


INTRODUCTION 


The  ubility  of  light  to  penetrate  the  sea  to  operationally  useful  depths,  the  ease  with  which  optical 
energy  cun  be  generated  and  modulated,  and  the  ubility  to  direct  this  energy  into  narrow  beams  utilizing 
small  antennas  (optical  systems)  illustrate  the  theoretical  feasibility  of  optical  communication  links  in  the  c 
sea.  For  exuniplc,  when  the  proportionate  sea  path  is  limited  to  lengths  of  hundreds  of  feet  (or  meters),* 
optical  communication  links  between  an  earth  orbiting  satellite  and  submerged  terminals  appear  possible. 

■ Before  such  optical  communication  links  become  practical,  phenomena  reported  in  the  literature 
must  be  quantified  to  reliably  predict  system  performance.  This  Is  necessary  for  systematic  design  and 
cvuluution  of  terminals.  For  example,  optical  beam  spreading,  attenuation,  uir/sea  interface  effects,  und 
background  radiant  intensities  must  all  be  documented  to  be  used  for  component  development. 

With  these  objectives  in  mind,  an  experiment  designed  to  measure  the  propagation  characteristics 
of  light  in  the  ocean  wus  conducted  In  the  vicinity  of  Sunta  Catalina  Island,  California,  between  June  and 
August  1975  (ref.  I ), 

9 F - 

1.1  PURPOSE 

The  purpose  of  this  experiment  was  to  examine  the  propagation  of  light  in  the  sea,  at  the  sca/air 
interface,  und  in  the  air  ubove  the  sea,  This  served  to  quantify  and  catalog  previous  measurements  and  to 
obtain  necessary  new  dutu.  These  data  will  be  used  as  the  basis  for  the  design  of  optical  communication 
links  operating  in  the  conglomerate  environment.  Specifically,  this  experiment  was  designed  to: 

a.  Meusurc  (in  the  atmosphere)  the  radiant  intensity  puttern  produced  by  an  underwater  optical 
projector  at  various  orientations  and  depths  whose  beam  transits  the  uir/sea  interface.  These  measurements 
were  made  in  the  fur  Held  of  perturbations  caused  by  wave  effects. 

b.  Measure  the  undorwuter  rndlant  Intensity  profile  produced  by  a simulated  extended  optical 
source  whose  beam  transits  the  nir/sea  interface.  The  definition  of  an  extended  optical  source  is  one  whose 
projected  dimension  at  the  air/sea  Interface  is  much  greater  than  the  depth  at  which  the  measurement  Is  made. 

Although  the  experiment  was  not  geared  to  any  specific  propagation  model,  the  propagation  model 
that  was  developed  hus  been  verified,  and  provides  a very  useful  tool  for  understanding  the  physical  require- 
ments of  transceivers  operating  In  the  ocean  environment.  The  verification  of  this  model,  which  is  high- 
lighted in  Section  4,  will  provide  a basis  for  the  engineering  design  of  communication  systems  and  a mech- 
anism to  predict  the  performance  of  optical  communication  links  when  penetration  of  the  sea  is  a require- 
ment. 

•Distance  measurements  throughout  this  document  are  mostly  in  feet.  However,  somo  measurements  are  presented  in  metric 
units,  while  others  Indicate  meters  followed  by  the  equivalent  in  feet,  For  those  metsurements  that  do  not  display  a com- 
bins t Ion  of  both,  conversion  factors  ere  provided:  Multiply  Feet  by  0.3048  to  obtain  meters,  ind  meters  by  3.2808  for 
the  value  in  feet. 


1.2 


BACKGROUND  ( ref.  I) 


The  Nliii.lv  aiul  progress  of  ocean  optics  li.iu*  been  closely  connected  with  the  evolution  both  ol’in- 
•'trutnentul  techniques  ;nul  of  motivational  factors  for  investigating  the  optical  character  of  the  seas. 

For  the  earliest  work,  which  dutes  back  to  1885,  only  photogruphic  methods  were  uvullubie  (refs. 

2 and  3),  The  advent  of  photocells  revolutionized  techniques  in  optical  oeeanogruphy  und  led  to  the 
development,  beginning  in  the  1930’s,  of  relatively  sophisticated  instruments  for  measuring  appurent  and 
inherent  optical  properties  of  the  sea, 

Throughout  the  history  of  underwater  optics,  u fundamental  Incentive  for  research  has  been  the 
photosynthetic  activity  of  the  ocean  in  terms  of  planktonic  ubundance  and  primary  organic  productivity 
(ref.  4),  In  particular,  studies  of  the  ecology  of  fish  have  made  wide  use  of  a subjective  measurement  tech- 
nique known  as  the  Scechi  disc  Although  the  properties  of  tlte  Seeelii  disc  have  never  been  standardized, 
Secchl  disc  measurements  represent  tite  most  complete  worldwide  body  of  knowledge  of  submarine  daylight 
and  wuter  transparency  available,  with  more  than  67,000  sightings  cataloged  (ref.  5). 

Another  long-term  motivational  factor  in  ocean  optics  has  been  its  relationship  with  meteorology 
und  the  physical  processes  occurring  at  the  sea  surface  which  determine  the  interchange  of  energy  between 
the  atmosphere  und  tlte  seu,  Tills  work  hus  stressed  modeling  of  radiative  transfer  in  the  sea  and  experi- 
mentally measuring  upwelling  und  downwclling  irrudiunces  in  the  oceans  (ref.  6). 

More  recent  applications  In  underwater  photography,  vision,  remote  viewing/surveillunoe,  and  in 
remote  detection  of  chlorophyll  from  uircraft  or  satellites  have  udded  Impetus  to  the  study  of  underwater 
natural  light  fields  (ref,  7). 

For  the  most  part,  there  exists  adequate  documentation  of  the  ocean  optical  propagation  proper- 
ties (ref.  8)  utilizing  the  sun  as  a source,  There  ure,  however,  new  applications  of  optics  where  datu  required 
to  predict  system  performance  in  the  ocean  are  either  scarce  or  nonexistent.  Examples  of  such  optical  sys- 
tems In  the  oceun  ure  those  which  operutc  through  the  air/water,  for  exumple,  an  optical  communication 
link  between  a satellite  und  a submerged  terminal. 

1.3  SCOPE  OF  WORK 

To  describe  opticul  systems  which  operate  in  the  environment,  quuntitutive  datu  ure  required  on 
underwater  rudlunt  intensity  distributions  from  natural  backgrounds;  radiance  distributions  from  artifi- 
cially collimated  sources  on  ••  xis,  off  axis,  and  operating  through  the  uir/watcr  interface;  and  the  optical 
properties  of  the  uir/wuter  Interface. 

Quantitative  descriptions  of  naturul  background  rudiunee  (or  special  rudiunce)  and  the  rudiunce 
distributions  ure  essentlul  for  determining  noise  levels  in  submerged  optical  equipment.  Although  absolute 
or  quantitative  data  have  not  be  :n  reported,  relative  datu  have  been  presented  by  several  workers  (ref.  8). 

The  only  work  which  specifically  studied  the  effect  of  propagating  a light  beam  through  the  air/water 
Interface  (Ohio  State  University)  employed  a small  diameter  beam  (conipured  to  the  surface  wave  struc- 
ture) and  examined  the  results  in  the  near  field  of  the  wave  optics  (ref.  9).  The  data  needed  to  understand 
most  systems  that  operute  through  the  surface  relate  to  beum  diameters  of  several  water  wave  lengths  and 
the  effects  observed  In  the  far  field  of  both  upward  and  downward  directed  energy, 

Some  Investigutors  Itave  published  results,  both  experimental  und  theoretical,  describing  underwater 
light  fields  from  Immersed  light  sources  (refs,  4,  10,  11,  12.  13,  und  14).  Unfortuniwely,  it  is  not  possible 
to  reconcile  the  vurlous  datu  sources  to  obtain  u meaningful  description  of  the  geometrical  beam  spreading 
expected  in  a collimated  light  bourn  for  predicting  the  performance  of  u system.  Little  work  has  been  reported 
on  the  optics  of  waves  In  terms  of  scattering,  reflection,  transmission,  und  their  effect  on  the  radiant  inten- 
sity pattern  thus  produced. 


I 

I What  is  needed  is  a workable  model  that  cun  be  used  to  describe  the  propagation  of  light  beams 

over  a water  path  and  through  the  air/sea  interface.  The  objective  of  the  experiment  herein  described  was 
to  obtain  the  measurements  necessary  to  verify  a model.  The  model  used  is  presented  in  Appendix  A. 

I The  results  of  this  experiment  will  find  immediate  applications  among  electro-optical  system 

designers  working  on  other  programs  that  must  operate  in  this  environment.  It  is  currently  necessary,  in 
at  least  one  udvunccd  development  program,  to  overdesign  and  construct  equipment  to  compensate  for 

I uncertainties  in  the  hydrosol  environment.  This  approach  has  proven  expensive,  and  the  equipment  de- 

signs are  complex  and  lurge.  Acquisition  of  the  necessary  optical  oceanographic  data  and  development  of 
a predictive  ocean  optics  model  will  result  In  substantial  savings  in  these  programs,  and  will  provide  a 
primary  mechanism  to  enhance  reliability  at  lower  system  costs. 


SECTION  2 


FORMULATION  OF  THE  EXPERIMENTAL  MODEL 


The  busic  theory  describing  the  propagation  oflight  through  any  multiple  scattering  medium  is 
radiative  transport,  References  on  this  topic  are  listed  in  Appendix  A Unfortunately,  closed  form  solutions 
exist  only  for  certain  special  cases.  The  solution  most  closely  matching  the  underwater  environment  is  the 
forward  scattering  solution,  in  which  the  scattered  energy  is  concentrated  in  a tight  cone  around  the  direc- 
tion of  origin  of  the  ray.  In  this  case,  the  mutual  coherence  function  is  known  to  have  the  general  form 


which  specializes  to  the  function 


M(p,z)  ■ e"*az)exp  [- 


when  an  empirically  derived  scatter  function  (ref,  15)  is  used.  The  0 ^ is  a measure  of  the  scattered  cone 
angle; s is  the  scattering  coefficient,  a is  the  absorption  coefficient;  p is  the  correlation  length;  X is  the  wave- 
length; und  z is  the  distance  propagated.  An  off  axis  correction  has  been  offered  (ref.  16)  to  account  for 
wide  angle  scattering.  This  solution,  which  is  presented  as  Appendix  A,  forms  the  basis  for  the  Optical 
Satellite  Communications  (OPSATCOM)  model.  In  this  section,  we  will  give  some  elementary  insight  into 
the  model  and  a discussion  of  how  the  model  was  adapted  for  reducing  the  data  taken  at  Santa  Catalina 
Island. 

It  is  shown  (ref.  17)  that  the  Fourier  transform  of  the  mutual  coherence  function  is  a function 
whose  rclutlve  amplitude  is  proportional  to  the  direction  of  arrival  of  the  intensity.  Thus,  by  examining  the 
mutual  coherence  function,  information  can  be  obtained  concerning  the  angular  distribution  of  the  arriving 
intensity.  For  example,  a constant  mutual  coherence  function  implies  a point  source,  whereas  an  impulsive 
mutual  coherence  function  is  generated  from  radiation  coming  from  all  directions  to  the  receiver.  It  is  also 
shown  that  equation  (2-1)  has  the  following  limiting  form: 


(pIpoT  ^ 

Mlp./.)  *e~aze  *■  : 0~  p~ 


sz  » 1 , 


p-(a+s)z 


; (f)2»V». 


; sz  « 1 


where 


p°  is;  3*- 

Thus,  we  see  that  for  short  distances  a point  source  will  retain  its  imuging  properties.  At  large  distances,  the 
mutual  coherence  function  is  predominantly  Gaussian  In  shape  with  an  asymptotic  value  of  exp  -{a  + s)z, 
whereby,  the  residual  Imaging  term  is  reduced  to  this  asymptotic  value.  Since  M(p,  z) « exp  -(az)  for  p - 0, 
we  see  that  the  total  power  at  u point : is  always  much  greater  than  that  contained  in  the  residual  (unscat- 
tered) term.  Specifically,  it  is  exp  (sz)  times  greater.  However,  this  abundant  power  is  associate^  mainly 
with  the  Gaussian  portion  of  the  mutual  coherence  function,  and  consequently  appears  to  come  from  a 
source  subtending  a large  field  of  view.  This  holds  true  whenever  the  receiver  is  in  the  scattering  medium. 

In  the  model  used,  the  energy  at  uny  point  is  constructed  from  all  the  rays  converging  there.  Thus, 
the  solution  is  derived  In  the  transform  domain  by  examining  the  divergence  of  a zero  cross  section  colli- 
mated beam  which  can  be  described  as  a spatial  impulse  response.  Once  a boundary  has  been  defined,  the 
solution  at  any  point  is  obtained  by  convolving  the  spatial  impulse  response  over  the  boundary.  The  spatial 
impulse  is  derived  as  a sum  of  two  terms.  The  first  term  corresponds  to  the  residual  image  (or,  in  this  case, 
the  point  source)  which  is  merely  reduced  by  the  factor  exp  -(a  + s)z.  The  second  term  in  the  spatial  im- 
pulse response  is  obtained  by  subtracting  the  residual  value  and  assigning  the  remaining  weight  to  the 
Gaussian  term,  which  becomes 
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The  model  takes  the  random  surface  into  account  by  assuming  that  the  point  source  is  Gaussian 
distributed  and  then  adds  the  appropriate  variance  to  the  angular  spread.  This  procedure  is  correct  on  an 
ensemble  basis.  However,  it  must  be  noted  that  when  the  point  image  is  discernible,  this  spread  is  traced  by 
the  motion  of  the  point  cuused  by  the  dynamics  of  the  surface.  Thus,  if  the  radiation  were  frozen  in  time, 
it  would  most  likely  be  seen  as  one  or  several  spots  coming  from  a direction  whose  angle  of  arrival  was 
Gaussian  distributed  in  time,  but  fixed  at  any  instant.  Although  we  have  not  done  so  in  this  case,  it  can  be 
represented  as  a point  or  sum  of  points  with  a time  varying  mean.  For  the  scattered  part  of  the  radiation, 
this  ensemble  averaging  has  already  been  performed  by  the  multiple  scattering  medium  and  is  of  no  concern, 

Another  addition  to  the  model  is  identified  as  the  glow  field  component.  This  term, as  described  in 
Appendix  C,  is  included  to  take  account  of  the  transition  region  existing  between  1 and  10  scattering 
lengths  (sz),  In  essence,  what  is  done  is  to  subtract  both  the  residual  image  and  the  scattering  field,  and  to 
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fit  a second  Guussiun  term  to  the  remainder.  This  gives  a better  fit  to  the  function  in  equation  (2-2) 
through  the  transition  region.  A final  modification  to  the  mojiel  is  described  in  Appendix  B,  This  modifi- 
cation accounts  for  the  fact  that  all  the  coefficients,  a , ,v.  and  02,  are  actually  functions  of:  in  any  real 
environment.  We  are  therefore  able  to  layer  the  model  to  accommodate  uny  variations  in  depth.  The 
presence  of  the  experimental  test  barge  cun  also  be  accounted  for  when  convolving  the  impulse  response. 
By  knowing  the  exact  orientation  and  coordinates  of  the  boundary,  this  area  is  given  zero  weight  during 
the  convolution. 


SECTION  3 


EXPERIMENTAL  METHODS  AND  APPARATUS 


3.1  METHODS  UTILIZED 

In  general,  the  atmospheric  radiance  pattern  of  an  underwater  optical  projector  was  measured  by 
overflying  the  surface  spot  illuminated  by  the  projector  witli  u tracking  optical  receiver.  The  orientation, 
field  of  view,  und  depth  of  the  projector  were  correlated  in  consonance  with  the  measured  oceanographic 
parameters  and  with  the  radiance  pattern  measured  during  the  surface  spot  overflights, 

Tiie  sun  wus  used  to  simulate  a point  source,  and  a hemispherical  scanning  receiver  was  employed  to 
meusure  the  radiance  profile  at  various  receiver  depUis  and  source  orientations. 

Both  the  uplink  and  downlink  data  obtained  in  open  ocean  waters  were  concurrently  calibrated  in 
terms  of  the  physical  characteristics  of  the  water.  Tills  calibration  was  accomplished  by  the  Visibility  Labo- 
ratory of  Scrlpps  Institution  of  Oceanography  whore  measurements  were  made  to  determine  the  extinction 
coefficient  (a),  scattering  coefficient  (s),  absorption  coefficient  (a),  and  the  volume  scattering  function 
(0(d)). 

The  collection  of  this  empirical  data  has  permitted  verification  of  the  model  presented  in  Section  2 
and  Appendix  A.  Using  this  model  it  will  be  possible  to  predict  link  parameters  which  influence  the  inter- 
cept contours  of  air  or  space  to  undersea, and  of  undersea  to  air  or  space  communication  links. 

3.2  EQUIPMENT  DESCRIPTION 

The  equipment  used  to  gather  data  during  this  experiment  can  be  separated  into  three  distinct 
areas:  the  underwater  equipment  package,  the  aircraft  equipment  package,  und  the  surface  support  plat- 
form, For  u more  complete  description  of  the  above  equipment,  refer  to  Volume  II,  Sections  2,  3,  and  6, 
respectively. 

3.2.1  UNDERWATER  EQUIPMENT  PACKAGE 

In  general,  the  underwater  equipment  package  is  a submergible  equipment  platform  which  is  con- 
trolled from  the  surface. 

Attuched  to  it  are  u dye  laser  (see  Volume  II,  Section  1 ) and  the  underwater  radiance  scanner  (see 
Volume  II,  Section  2). 

3.2.1 .1  THE  DYE  LASER.  The  dye  laser  (figure  3-1),  which  wus  mounted  on  the  underwater 
platform,  was  housed  in  a watertight  container  and  controlled  from  the  surface.  It  could  be  lowered  on  the 
platform  to  a predetermined  depth  and  activated  at  will.  The  mount  allowed  360®  azimuth  rotation  with 
0°  to  50°  variances  in  zenith  angle.  The  azimuthal  rotation  was  continuous,  while  the  zenith  angle  was 
controlled  discretely  in  10°  steps. 


Tlu;  laser  driver,  which  controlled  the  activation  of  the  laser,  also  output  a trigger  voltage  that  acti- 
vated a surface  mounted  radio  transmitter.  Tins  radio  transmitter  was  used  to  synchronize  time 
with  a tracking  receiver  that  was  mounted  in  an  aircraft  and  flown  over  the  experimental  test  site  to 
measure  the  upward  radiance  pattern. 

The  pertinent  specifications  of  the  laser  were: 


Pulse  repetition  rate 
Peak  pulse  amplitude 

Pulse  width 
Center  wavelength 
Bandwidth 

Beumwidth  (in  water) 


20  pulses/ sec”* 

5 kW  (±6%  relative 
±10%  ubsolute) 

0.7 5 psec 
5,2 14  A 
46.5  A 
1.23s 


3.2.1 .2  UNDERWATER  RADIANCE  SCANNER.  The  underwater  radiance  scanner  (figure  3-2), 
which  is  essentially  un  electronic  video  camera,  image  dissector  tube,  was  mounted  on  the  underwuter 
platform.  By  scanning  its  focal  plane,  it  discretely  sampled  the  composite  field  of  view,  including  the 
upper  hemisphere.  The  scanner  was  mounted  Inside  a watertight  container  (figure  3-3),  with  an  acrylic 
transparent  dome  and  a flshcye  lens  (figure  3-4),  The  flsheye  lens  was  Installed  so  that  its  focal  plane 
was  the  photoemissive  surface  of  the  Image  dissector  tube. 

Appendix  K contains  a program  for  the  implementation  of  calibration  taken  by  the  underwater 
radiance  scanner, 

The  platform  could  be  lowered  to  any  desired  depth.  Computer  control  from  the  surface  support 
platform  allowed  the  field  of  view  to  be  scanned  in  segments  of  approximately  1°  by  1°.  Each  field  of  view 
segment  was  uniquely  identified  and  the  measured  radiance  was  recorded  as  raw  data. 

A narrow  spectral  filter  of  93  A bandwidth  centered  at  5,200  A was  incorporated  Into  the  optic 

train, 


3.2.2  AIRCRAFT  EQUIPMENT  PACKAGE 

The  aircraft  equipment  package  consisted  of  a calibrated  tracking  optical  receiver,  Its  electronic 
control  circuitry,  and  a digitul  data  recording  system.  A complete  description  of  the  aircruft  receiver  sys- 
tem is  contained  in  Volume  II,  Section  3. 

The  trucking  optical  receiver  (figure  3-5)  was  designed  to  measure  the  radiance  distribution  of  an 
underwater  laser  at  vurious  depths  and  orientations.  Measurements  were  made  while  the  receiver  passed 
overhead  in  its  aircraft  plutfonn  (figure  3-6)  and  trucked  the  optical  source,  the  dye  laser.  The  receiver  was 
manually  controlled  by  a joystick  pointing  mechanism  that  was  pointed  at  the  surfuce  of  the  water  above 
the  laser  source.  When  sufficient  energy  was  detected,  the  receiver  automatically  tracked  the  laser  signal 
(figure  3-7)  and  recorded  on  magnetic  tape  the  signal  intensity,  and  the  receiver  pitch  and  roll  angles.  Using 
these  data  and  known  aircraft  altitude,  one  can  determine  the  cuts  that  were  taken  of  the  desired  upward 
radiance  profiles.  Changing  the  laser  depth  or  angular  orientation  and  repeating  the  above  procedure  yields 
a family  of  radiance  profiles,  as  a function  of  the  transmitter  configuration  and  the  water  parameters  which 
were  measured  simultaneously. 

The  operational  characteristics  of  the  aircraft  optical  tracking  receiver  were: 


•a  nsitivity 
f /# 


4.7  X 10”9  w/cm"2 
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Figure  3-4.  Closeup  view  of  acrylic  dome  (with  fliheye  lens  of  the  underwater  radiance  scunner  visible), 


3-6 


‘ " * 1 i ta 





Figure  3-5.  Trucking  optical  receiver  In  laboratory  teat  mount. 
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Figure  3-6.  Tracking  optical  receiver  mounted  In  PBY  aircraft. 


Figuro  3>7.  Aircruft  over  flight  of  the  experimental  teat  site  (with  tho  tracking  optical  receiver  automatically 
tracking  the  lubmorged  laier  tuurco). 
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Lens  diameter 

6" 

Spectral  bandpass 

5.240  ±38  A 

Data  field  of  view 

1.4° 

Visual  field  of  view 

15° 

Maximum  roll  unglc 

±22.5° 

Maximum  pitch  angle 

H- 

V© 

O 

o 

Roll  rate  correction 

> 5°  sec" 

Pitch  rate  correction 

> 5®  sec"1 

Joystick  slew  rate 

> 15°  sec-1 

A radio  link  between  the  surface  support  platform  and  the  aircraft  provided  timing  informa- 
tion that  was  synchronous  with  the  tracking  optical  receiver. 

3,2.3  SURFACE  SUPPORT  PLATFORM 

Tiie  surface  support  platform  (figure  3*8)  provided  stubility  for  conducting  the  experiments  in  an 
exposed  ocean  environment  under  nonsevere  weather  conditions.  A complete  description  is  presented  in 
Volume  11,  Section  6.  It  hud  a mechanism  in  the  form  of  an  assembly  which  could  be  raised,  lowered,  or 
rotated  to  provide  vertical  access  to  depths  of  approximately  164  feet,  and  upon  which  tiie  underwater 
radiance  scanner  and  the  dye  laser  could  be  mounted  (figure  3*9).  The  platform  also  provided  the  means  to 
traverse  un  optical  detector  ulong  u horizontal  path  slightly  beneath  the  surface  of  the  water.  Necessary  power 
supplies,  personnel  accommodations,  a control  equipment  shelter,  oceanographic  instruments,  support  hut.  and 
operating  equipment  to  facilitate  experimental  data  gathering  operations  were  provided.  Also  mounted  on  this 
platform,  were  the  following  environmental  sensing  instruments,  which  ure  described  in  Volume  11.  Section  7, 

a)  Solar  monitor  - a narrow  Held  of  view  (7°)  calibrated  receiver,  filtered  identically  to  the 
underwater  radiunce  scanner,  which  was  used  to  monitor  direct  sun  rays. 

b)  Deck  cell  - a lumbertlan  collecting  thermopile  filtered  at  5,220  A center  wavelength  with  an 
820  A bandwidth  that  was  used  to  monitor  ambient  background  levels, 

c)  Wind  speed  and  direction  measuring  instruments. 

3.3  EXPERIMENTAL  TEST  SITE 

The  test  site  was  located  in  waters  adjacent  to  Santa  Catalina  Island.  The  criteria  for  tills  selection 
are  described  in  ref,  l. 

This  test  site  was  near  NELC  and  provided  appropriate  water  depths.  Water  clarity  was  comparable 
to  that  of  open  ocean  areas.  Laboratory,  dockage,  and  personnel  support  facilities  were  available  at  the 
Santa  Catalina  Murine  Biological  Research  Station,  operated  by  the  University  of  Southern  California. 

Figure  3-10  shows  the  test  site,  the  approximate  location  of  the  Marine  Biological  Research  Station,  and 
the  aircraft  overflight  path  that  was  used  in  conducting  the  uplink  measurements. 

3.3.1  AZIMUTHAL  ALIGNMENT  OF  THE  UNDERWATER  RADIANCE  SCANNER 

In  order  to  align  the  azimuth  of  the  underwater  radiance  scanner  (URS)  with  the  earth’s  coordinate 
system,  a precise  mark  on  the  URS  was  aligned  with  a prominent  projection  of  the  coust  of  Siuita  CataJinu 
Island  at  Fisherman  Cove  so  that  the  Y-axis  was  oriented  north-south.  To  calibrate  this  alignment,  the  URS 
was  operated  out  of  the  water  at  a known  time,  Tills  Imaged  the  sun,  slightly  out  of  focus,  on  the  image 
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I'lguro  3-8.  lixporlmotuul  ton  barge  (with  180  foot  guide  pipe  suipendod  from  right-hand  comer  and  undorwuter 
platfonn  In  raised  position). 
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plane  which  was  then  scanned  by  the  URS  image.  dissector  tube.  The  location  of  the  peak  radiance  in  the 
image  plane  yielded  the  angulur  locution  of  the  sun  with  respect  to  the  URS.  By  simultaneously  determin- 
ing the  angular  orientation  of  the  sun  with  respect  to  the  known  geographic  location  of  the  URS,  it  was 
possible  to  determine  the  alignment  of  the  URS  in  relation  to  true  North  and  the  zenith. 

A summary  of  these  measurements  follows: 

(a)  URS  location,  33°27'I4"  N,  1 18°28'55"  W. 

(b)  Time  of  observation  2224  Greenwich  mean  time  24  June  1975. 

(e)  Calculated  Observed 

azimuth  262.5°  263.0° 

elevation  56.2°  54.72° 

This  measurement  confirms  that  the  angular  alignment  was  within  the  experimental  goal  of 
±3°. 


3.3.2  GENERAL  WATER  CHARACTERISTICS 

Tlie  test  site  was  located  in  waters  adjacent  to  the  Southern  California  coast  generally  classified  as 
coastal  waters.  This  water  lias  its  maximum  transmission  peaking  in  the  480  to  550  nm  wavelength  region 
of  the  visible  spectrum.  Typical  attenuation  lengths  in  terms  of  meters/ln,  the  distance  in  which  the  intensi- 
ty is  reduced  by  a factor  of  exp  (-1),  are  shown  in  figure  3-11  which  is  extracted  from  ref.  10.  For  the 
experiment  described  here,  a laser  source  of  wavelength  5,214  A was  chosen  and  the  URS  was  spectrally 
filtered  about  u center  frequency  5,200  A wavelength. 

The  test  site  was  somewhat  protected  from  the  prevailing  seas  by  the  northeastern  projection  of 
Santa  Catalina  Island.  During  the  experiment,  the  sea  state  ranged  from  a flat  calm,  with  swells  of  1 to  2 
feet  from  the  northwest,  to  seas  of  5 to  8 feet  with  patches  of  wind-generated  whltecaps  covering  approxi- 
mately S percent  of  the  ocean  surface. 

3.4  EXPERIMENTAL  SCENARIOS 

The  experiment  made  three  measurements  to  describe  the  propagation  of  light  in  the  ocean  environ- 
ment: Downlink,  uplink,  and  F(0).  Of  these  three  measurements,  approximately  99  percent  of  the  experi- 
mental time  was  spent  makhy  the  downlink  and  uplink  measurements  which  were  considered  the  most 
important.  The  F(0)  data  is  presented  in  Appendix  M.  Following  is  a brief  description  of  the  experimental 
scenarios  that  were  used  to  collect  the  data  analyzed  that  is  presented  in  Section  4. 

3.4.1  DOWNLINK  MEASUREMENTS 

The  sun  was  used  as  the  downlink  source  because  as  a stable  source  of  sufficient  inten- 
sity, it  adequately  represented  a satellite  transmitter,  The  use  of  an  airborne  laser  as  the  source  was  rejected 
early  in  the  experiment  for  the  following  reasons.  To  satisfy  the  simplified  range  equation,  the  radius  of 
the  surface  spot  is  required  to  exceed  half  the  depth  of  the  receiver  (see  Appendix  A,  figure  10).  If  a spot 
of  such  dimensions  (c.g.  164  feet),  is  projected  from  an  airborne  laser  at  reasonable  altitudes  of  3,280  feet 
or  less,  five  wavefront  departs  substantially  from  a plane  wave.  In  addition,  flying  and  pointing  a laser  from 
an  aircraft  is  a difficult  and  expensive  task. 

The  solar  intensity  was  monitored  continuously  at  the  sea  surface  to  correct  the  data  for  any  chang- 
es in  source  intensity.  Recorded  data  were  normalized  to  these  readings,  simultaneously  taken,  so  that  the 
results  could  be  expressed  in  terms  of  loss  per  stcradian  as  a function  of  depth. 
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Figure  3-12  depicts  a downlink  recording  event.  Using  the  sun  us  the  optical  source,  the  solar  cell 
(or  deck  cell)  wus  used  to  measure  the  solar  intensity;  Scripps  Institution  of  Oceanography  was  at  the  test 
site  providing  water  calibration  measurements  (a  and  s calibration),  and  the  UR5  was  mounted  to  operate 
vertically  on  the  180  foot  guide  pipe  suspended  from  the  southern  corner  of  the  experimental  test  barge. 

Dutu  wus  recorded  ut  approximately  extinction  length  intervals  down  to  a depth  of  1 50  feet.  The 
URS  was  used  to  record  a hemispherical  mapping  of  the  radiant  intensity  profile  produced  by  the  solar 
source.  The  tilt  ungies  of  the  underwater  platforms  were  monitored  to  determine  both  vertieul  and  azi- 
mutlial  references. 

Measurements  were  conducted  under  predominantly  cloudy  conditions,  however,  some  were  under 
clear  conditions.  These  are  described  in  detail  in  Section  4.  Figure  3-13  represents  the  zenith  and  uzimuth- 
al  angle  trajectory  of  the  sun  us  recorded  at  the  test  site  during  the  experiment.  This  trajectory  permitted 
zenith  angle  opportunities  of  from  approximately  10s  to  greater  than  90°  over  an  azimuthal  range  of  ap- 
proximately 50°  to  310°  referenced  to  true  North. 

The  suspension  of  the  underwater  platform  and  Its  guide  pipe  left  a clear  field  of  view  to  the  sun  at 
all  angles  during  the  day,  and  provided  unobstructed  views  of  sunrises  and  sunsets. 

Experimental  data  was  recorded  in  the  following  way: 

(a)  The  URS  was  lowered  to  a known  depth. 

(b)  The  URS  was  exercised  to  set  its  gain  characteristics,  wliich  were  then  recorded  in  the  test  log. 

(c)  The  deck  cell  (solar  monitor  instruments)  readings  were  recorded. 

(d)  Tiie  URS  wus  programmed  to  execute  a preset  quantity  of  scans  and  was  then  activated,  thus 
automatically  recording  the  data  in  digital  format  on  magnetic  tape  for  later  processing. 

(e)  The  time  of  each  scan  was  both  manually  and  automatically  recorded  individually  or  by  series. 

(0  The  URS  was  then  lowered  or  raised  to  a new  depth  and  the  above  procedure  was  repeated. 

3.4.2  UPLINK  MEASUREMENTS 

The  dye  laser  was  used  as  the  underwater  source  for  the  uplink  tests.  The  tracking  optical  receiver 
was  mounted  in  an  aircraft  to  record  data  while  flying  over  this  source, 

The  short  pulses  of  the  underwater  laser  required  the  use  of  a peak  reading  circuit  to  hold  the  level 
of  the  pulse  until  the  next  pulse  occurred.  The  output  of  this  circuit  was  recorded  on  one  channel  of  the 
9-track  digital  tape  recorder.  A timing  signal  was  recorded  for  reference  purposes.  As  the  signal  was 
received,  the  pointing  angles  of  the  tracking  optical  receiver  were  recorded.  These  data  were  later  re- 
duced and  are  presented  in  Section  4,  paragraph  4.2  of  this  document. 

Figure  3-14  presents  a pictorial  representation  of  an  uplink  data  gathering  event.  This  event  differs 
from  the  downlink  in  that  the  dye  laser,  which  was  mounted  on  the  underwater  platform,  was  activated  at 
various  depth  and  zenith  angle  orientations.  As  the  aircraft  overflew  the  test  site,  the  tracking  optical  re- 
ceiver measured  the  radiant  intensity  produced  by  the  laser  source.  All  background  and  water  calibration 
measurements  that  were  conducted  during  the  downlink  events  were  repeated  during  the  uplink  episodes 
with  the  exception  of  simultaneous  operation  of  the  URS.  A radio  link  between  the  laser  con- 
trol point  and  the  aircraft  passed  time  gating  signals  to  syneltronize  the  optical  receiver  and  laser 
transmitter. 
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3.4.3  W MEASUREMENTS 

The  F(0)  measim  ments  recorded  the  spatial  impulse  response  of  the  dye.  laser.  In  order  to  make 
this  measurement,  an  optical  receiver  (described  in  Volume  II,  Section  4)  wus  transited  just  beneath  the 
water  surface  along  a specially  constructed  l-beam  attached  to  the  side  of  the  barge  (see  figure  3-14).  This 
receiver  was  mounted  so  that  it  could  be  stopped  at  any  position  along  the  optical  bench  and  pivoted 
through  160°  (±80°  about  the  nadir)  in  the  plane  which  contained  the  underwater  platform  suspension  pipe 
and  the  F (0 ) optical  receiver. 

Data  were  gathered  as  follows; 

(a)  The  dye  laser  was  placed  in  a preselected  angular  orientation  with  respect  to  the  vertical  and 
was  lowered  to  a predetermined  depth  so  that  it  transmitted  in  the  plane  that  contained  the  laser,  the  pipe, 
and  the  F(b)  optical  receiver. 

(b)  The  F(0)  optical  receiver  was  positioned  at  a known  horizontal  position  away  from  the  laser. 

(c)  The  angular  position  of  the  F(0)  optical  receiver  was  set,  the  dye  laser  was  activated,  and  the 
signal  level  as  received  was  recorded. 

(d)  After  the  angular  position  of  the  F(0)  optical  receiver  had  been  exercised  through  its  range, 
the  horizontal  position  was  changed,  and  the  sequence  was  repeated. 

(e)  When  a complete  sequence  of  measurements  had  been  made  that  sufficiently  exhausted  the 
receiver  angular  and  horizontal  position  possibilities,  the  laser  was  lowered  to  a different  depth,  its  angular 
orientation  was  changed,  or  a combination  of  both.  This  different  configuration  then  formed  a new  set  of 
experimental  knowns,  and  steps  (b)  through  (d)  were  repeated. 


SECTION  4 


DATA  AND  RESULTS 


The  experiment  consisted  of  two  major  portions:  uplink  and  downlink  measurements.  In  addition, 
several  levels  of  calibration  were  employed.  The  calibration,  which  is  discussed  in  detail  in  Volume  II,  con- 
sisted of  equipment  calibration,  water  calibration  at  the  fundamental  level,  and  water  calibration  at  the 
system  level.  The  equipment  calibration  was  used  to  determine  the  absolute  levels  of  the  uplink  and  down- 
link measurements  in  order  to  obtain  quantitative  results.  This  resulted  in  absolute  accuracies  of  3 to  5 dB. 
The  water  calibration  was  necessary  as  an  input  to  model  development.  These  measurements  were  per- 
formed by  Scripps  Institution  of  Oceanography  with  the  results  appearing  in  Appendix  D and 
Volume  II,  Section  5. 

4.1  DOWNLINK  MEASUREMENTS 

The  downlink  measurements  that  were  performed  by  the  URS  are  described  briefly  in  Section  3 
and  in  detail  in  Volume  11,  Section  2,  This  Instrument  was  used  to  make  tolar  measurements  in  a 95  A 
bandwidth  around  5,200  A.  This  scenario  closely  describes  that  of  a satellite  except  that  there  is  no  modu- 
lation and  the  sun  Is  not  a point  source,  but  rather  subtends  an  angle  of  0,5°.  Neither  of  these  deficiencies 
mattered  as  far  as  the  goals  of  this  experiment  were  concerned.  The  data- recording  equipment  was  pre- 
viously described  in  Section  3,  Two  data  acquisition  programs  were  used  ind  are  described  in  Appendices 
B and  H,  and  in  Volume  11,  Section  2.  The  data  taken  with  the  quick  scan  program  are  summarized  in 
figure  4-1 . Representative  curves  are  shown  in  figures  4-2  and  4-3,  For  eiich  data  point  plotted,  there  is  an 
accompanying  hemispherical  radiant  pattern,  togother  with  an  integration  of  the  power  in  this  pattern  as  a 
function  of  the  field  of  view  (figure  4-2).  Notice  the  power  loss  in  the  upper  left  hand  quadrant  of  figure 
4-3  resulting  from  the  shadow  cast  by  the  supporting  pipe.  The  origin  of  the  coordinate  system  corresponds 
to  the  zenith,  while  edges  are  90°  from  the  zenith.  The  patterns  are  truncated  whenever  the  noise  level  is 
reached  so  that  ull  the  plots  do  not  display  the  horizon.  The  datu  points  In  figure  4-1  correspond  to  the 
asymptotic  value  of  the  integrated  curves  in  figure  4-2  which  correspond  i;o  the  total  Irrudiarice  (intensity) 
incident  upon  the  radiance  scanner,  All  curves  are  normalized  to  the  surface  irradiancc  so  that  we  have 
units  of  loss/steradlan  for  radiance  und  loss  for  the  integrated  power. 

4.1.1  DATA  CALIBRATION 

Notice  that  at  40  feet  and  below  anomalies  exist  In  the  duta.  Thi  s is  better  observed  by  focusing  on 
a single  set  of  duta  taken  us  the  radiance  scanner  was  submerged  (sec  figure  44A).  Observe  the  10  dB  dis- 
continuity occurring  at  40  feet.  It  is  our  opinion  that  there  was  an  inter  mittent  malfunction  regarding  the 
gain-determining  mechanism  which  was  depth  dependent.  (The  underwater  connector  failed  late  in  the 
experiment  and  might  have  been  the  trouble,)  We  support  this  claim  in  'he  following  manner.  If  10  dB  is 
subtracted  from  each  of  the  data  points  below  40  feet,  and  the  resultant  curve  compared  to  a K-meter 
measurement  of  Irradlance  taken  simultaneously  by  John  Shannon  of  the  Naval  Air  Development  Center, 
then  the  comparison  In  figure  44B  can  be  made.  Notice  that  the  greateiit  deviation  is  only  3 dB  which 
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occurs  at  the  shallow  depths.  R.  W.  Austin  of  Scripps  Institution  of  Oceanography  indicated  (in  private 
conversation)  that  based  upon  nis  extensive  experience  in  K-meter  measurements,  the  behavior  in  figure 
4-fB  is  what  lie  would  expect;  namely,  a lower  value  for  the  radiance  scanner  at  shallow  depths  with  an 
asymptotic  convergence  at  the  deeper  depths.  He  also  expected  the  results  at  the  deeper  depths  to  be  inde- 
pendent of  whether  it  was  a sunny  or  an  overcast  day.  A suggestion  has  been  made  at  NELC*  that  the 
10  dB  correction  should  be  made  by  increasing  the  value  of  the  shallow  depth  measurements.  It  is  believed 
that  this  would  be  incorrect  for  two  reasons.  First,  this  would  result  in  a systematic  10  dB  discrepancy 
between  the  radiance  measurements  and  the  K-meter  measurements.  Since  the  latter  are  taken  with  a single 
instrument,  this  would  imply  a 10  dB  error  in  calibration.  Second,  and  more  compelling,  is  the  fact  that  the 
greatest  loss  occurring  at  150  feet  on  June  20  and  24  was  19,3  dB  without  the  lOdB  correction.  However, 
the  absorptive  loss  alone  when  projected  from  measurements  of  the  absorption  coefficient  taken  by  Scripps 
is  at  least  22.5  dB  for  these  two  days.  Allowing  for  a 10  percent  error  in  the  measurement  of  the  ubsorption 
coefficient  still  results  in  an  absorptive  loss  of  at  least  20.5  dB.  Examining  consecutive  scans  of  the  quick 
scan  program  at  a 1 50  foot  depth  shows  a peak  deviation  of  only  0.3  dB  in  irradlance.  Consequently,  the 
proposed  correction  would  violate  the  laws  of  thermodynamics.  What  is  important  to  point  out  is  that  the 
only  difference  between  the  K-meter  and  the  radiance  scanner  is  that  the  former  has  a cos  0 pattern  while 
the  pattern  of  the  latter  is  sin  8/8.  Hence,  within  a maximum  2.16  dB  factor,  both  measurements  are  identical. 

Finally,  figure  4-4B  plots  the  results  of  the  propagation  model  with  the  water  calibration  data 
inserted  for  the  same  time  and  day.  This  appears  to  give  confirmation  to  within  5 dB.  It  is  interesting  to 
point  out  that  the  model  seems  conservative;  so  that  projection  to  Jerlov  II  water  made  in  figure  1 1 of 
Appendix  A should  be  accurate. 

By  contrast,  our  confidence  in  the  repeatability  or  accuracy  of  the  Automatic  Hemispherical  Scan 
(AHS  I program  is  low.  This  can  best  be  explained  with  reference  to  figure  4-5  where  the  irradiance  vaiues 
obtained  from  the  AHS  are  displayed.  The  measurements  made  at  50  feet  were  recorded  in  22  consecutive 
scans,  Tlie  irradiant  loss  ranged  from  0.7  dB  to -12,7  dB,  or  peak  to  peak  of  13.4  dB.  This  range  in  irrHdiance 
values  supports  the  earlier  conjecture  of  an  equipment  malfunction.  Furthermore,  it  took  1 hour  and  15 
minutes  to  make  these  measurements  since  each  scan  took  approximately  4 minutes.  (By  contrast,  the 
quick  scan  progruni  on  June  20  took  26  scans  in  one  minute  at  150  feet  with  0.3  dB  peak  to  peak  deviation.) 
As  a consequence  further  discussion  is  concentrated  on  the  quick  scan  measurements  with  the  inclusion  of 
the  10  4B  correction  as  proposed.  Table  4-1  shows  a complete  duta  summary  of  the  experiment;  table  4-2 
reports  each  of  the  quick  scan  measurements;  and  table  4-3  is  a complete  report  of  the  Automatic  Hemis- 
pherical Scan  measurements. 

4.1.2  INTERPRETATION  OF  DATA 

As  evidenced  by  tables  4-2  and  4-3,  the  bulk  of  the  data  were  taken  on  June  '24  and  June  26.  As  it 
was  decided  not  to  use  the  Automatic  Hemispherical  Scan  data  because  oflarge  variances,  and  since  most 
of  the  June  26  data  were  taken  with  the  AHS  program,  we  will  concentrate  on  the  June  24  data.  These 
data  consist  of  two  distinct  portions.  The  portion  between  1 1 1 1PDT  and  1 156PDT  was  taken  during 
heavy  overcast  conditions  (refer  to  paragraph  4.1.1).  The  airport  at  Santa  Cutalinn  Island  reported  a cloud 
thickness  of  800  to  1,500  feet.  We  consider  this  set  of  data  to  he  well  calibrated.  Figure  4-3  displays  the 
hemispherical  radiance  pattern  of  the  data  at  150  feet.  Figures  4-6A  through  O represent  cuts  of  the 
radiance  profiles  taken  through  the  sun  angle  (Appendix  G),  Also  plotted  are  cuts  through  the  radiance 
peak,  and  the  model  developed  with  the  Scripps  data  inserted.  Since  it  was  an  overc.st  day,  the  sun  could 
not  he  considered  as  a point  source,  Instead,  an  initial  spread  (0J,)1  of  45°  was  used  to  simulate  the  isotropic 
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Figure  4-6C.  Radiance  profile  through  tun  angle. 
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nuture  of  the  source.  (A  cos  0 puttern  is  approximately  460°  with  a contraction  of  3/4  because  of  the 
watei  index  ) This  correction  seems  tc  be  better  behaved  at  the  deeper  depths.  Because  little  data  were 
obtained  on  u sunny  day,  it  wun  difficult  to  make  a comparison  between  the  model  and  the  data.  However, 
tv>o  good  runs  were  tnude  on  June  20,  one  ut  150  feet  and  the  other  at  20  feet.  At  1 50  feet,  a deviation  of 
approximately  5 dB  at  the  peak  of  the  radiance  pattern  was  noted.  However,  away  from  the  peak  there 
wus  good  agreement.  This  wux  interpreted  to  mean  that  the  shallow  water  correction  was  truncated  too 
quickly.  The  model  was  rerun  with  the  original  function  for  A,  in  equation  (€- 1 2).  A noticeable  improve- 
ment wus  seen  as  evidenced  in  figure  4-7 . For  the  20  feet  dutu,  both  A and  A*  produced  identical  results 
which  is  in  close  agreement  with  the  measured  datu  in  figure  4-8.  The  model  was  also  rerun  at  1 SO  feet 
with  the  new  value  of  A for  the  June  24  data  (figure  4-60)  and  u slight  Improvement  was  noticed.  Because 
of  the  lack  of  data,  It  was  decided  to  discontinue  model  reruns  However  it  it  believed  that  u shallow  watoi 
correction  can  be  empirically  determined  to  further  reduce  the  discrepancy,  What  is  more  important  for 
system  design  is  the  integrated  power  us  a function  of  the  field  of  view,  This  is  displayed  in  figure  4-9  A 
through  N for  the  June  24  data  together  with  upper  and  lower  bounds  derived  in  Appendix  A.  Notice  that 
there  is  a 0 to  5 dB  discrepancy  throughout  the  range,  with  the  larger  discrepancies  occurring  at  the  asymp- 
totic values,  Also  note  that  the  3 dB  field  of  view  is  approximately  ±30°.  This  is  true  even  on  a sunny  day 
us  seen  in  figure  4-10.  Here  the  discrepancy  between  datu  and  model  is  also  smaller.  We  also  point  out  that 
the  ocean  roughness  which  did  not  vary  greatly  wus  estimated  using  data  taken  on  June  23.  Here,  the 
radiance  scanner  was  submerged  to  1 foot,  and  1 00  consecutive  scans  were  taken.  As  pointed  out  in 
Appendix  A,  the  ensemble  average  of  all  these  scans  should  approach  a Gaussian  distribution  whose 
variance  is 

■ | 1 ~p|  var[R]  . 

The  ensemble  average  Is  plotted  In  figure  4-1 1 and  the  resulting  value  for  var(R)  is  ,044. 

The  second  portion  of  the  June  24  measurements  concerns  uncalibrated  data  taken  with  the  sun  ut 
or  belo  w the  horizon.  These  measurements  were  taken  to  simulate  the  effects  of  ship-submarine  line  of 
sight  and  over-the-horizon  scatter  modes  of  communication.  Figure  4-12  portrays  a radiance  profile  taken 
by  the  radiance  scunner  ut  50  fret  with  tire  sun  precisely  on  the  horizon.  Notice  the  large  solid  angle 
around  the  Snell's  angle  over  which  the  intensity  is  upproxlmutciy  constant.  In  figure  4-13,  there  is  u 
radiunce  profile  with  the  rudluncc  scunner  submerged  to  20  feet,  the  surface  Irradtunce  set  equal  to  unity, 
and  the  sun  15°  below  the  horizon.  Note  the  concentration  of  power  at  the  S nil’s  angle.  In  order  to  ob- 
tain a rough  calibration,  the  radiunce  scunner  wns  taken  out  of  the  water  and  the  sunset  viewed  directly 
(figure  4-14),  Now  the  concentration  of  power  falls  at  the  horizon.  Tire  total  integrated  power  In  figure 
4-14  wus  only  13  dB  greuter  than  the  integrated  power  In  figure  4-13,  Furthermore,  the  peuk  power  on  the 
horizon  in  figure  4-14,  was  only  6 dB  greater  than  that  of  figure  4-13,  Thus,  it  is  concluded  that  the  lot*  in 
power  und  the  Inability  to  enlibrute  were  due  to  the  absence  of  5,200  A In  a red  sunset  and  not  due  to  loss 
through  the  wuter,  To  emphasize  this,  irradtunce  vs.  zenith  angle  were  plotted  ut  each  depth  for  all  of 
our  data  (Appendix  L).  The  most  onlightening  of  these  curves  is  shown  in  figure  4*15.  The  results  of  the 
model  were  also  plotted  to  substantiate  the  invurlance  of  loss  witli  the  zenith  angle  at  a depth  of  100  feet. 
These  results  are  most  important  for  any  system  communicating  to  a submarine  in  the  blue -green  portion 
of  the  visible  spectrum. 

4.1.3  CLOUD  PENETRATION 

Although  no  attempt  was  made  to  record  quantitative  data  on  tight  penetration  through  clouds,  u 
pyrheliometer  was  used  to  continuously  monitor  the  lrradiance  at  the  ocean  surface.  The  output  from  this 
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device  wus  recorded  on  u strip  chart  with  a few  representative  samples  shown  in  figures  4-1 6A  through  C. 
Figure  4-1 6(’  is  a trace  for  a sunny  day.  By  using  this  curve  as  u reference,  it  was  possible  to  obtain  mean- 
ingful statistical  data  on  cloud  penetration.  Thus,  for  example,  figure  4-1 7A  is  a distribution  of  points 
taken  every  hour  between  1 100  and  1 500.  The  value  of  each  point  is  divided  by  the  corresponding  value  of 
the  sunny  day  curve  at  the  same  time  of  day,  This  is  designated  as  the  transmission  and  is  always  less  than 
or  equal  to  one.  The  transmission  vulues  are  plotted  on  the  abscissa  and  the  relative  frequency  (number  of 
times  a particular  value  occurs  divided  by  the  total  number  of  values)  is  plotted  on  the  ordinate.  This  bar 
churt  is  un  estimute  of  the  probability  distribution  for  transmission.  It  is  bimodal,  that  is,  an  impulse  occurs 
between  .9  and  1 .0,  and  a bell-shaped  portion  between  0 and  .9.  Of  a total  of  36  points,  21  represent  clear 
conditions  or  a cloud  -free  line  of  sight  probability  of  .58,  The  cumulative  probability  was  superimposed  to 
demonstrate  tills  better.  However,  36  points  do  not  represent  a significant  sample  so  the  same  procedure 
was  repeated  every  10  minutes  and  plotted  in  figure  4-17B.  Notice  the  same  basic  shape,  but  more  filled  in. 
The  same  data  were  reprocessed  every  10  minutes  between  0900  and  1700,  as  plotted  in  figure  4-17C,  which 
gave  a more  continuous  sampling  of  the  same  time  period.  Notice  that  there  is  no  basic  change  to  the  curves 
with  the  exception  that  the  probability  of  a cloud-free  line  of  sight  is  now  .62.  Therefore,  now  that  274 
points  have  been  obtained,  It  represents  a significant  set  of  data  and  some  conclusions  may  be  reached,  For 
example,  the  probability  of  having  less  than  a 10  dB  loss  is  ,986,  or  98.6%.  Even  the  conditional  probability 
(the  probability  when  there  is  a cloud,  referring  only  to  the  bell-shaped  part  of  the  curve)  of  having  less  than 
a 10  dB  loss  given  an  overcast  condition  is  96.2%.  Although  only  a spot  sample  was  taken  of  all  clouds,  it  is 
clear  that  transmission  through  clouds  with  nominal  values  and  at  high  probabilities  is  feasible.  This  conten- 
tion is  substantiated  since  the  mqjor  set  of  water  penetration  measurements  were  made  in  overcast  conditions 
so  thut  the  penetration  of  diffuse  light  into  the  water  is  understood.  The  extension  of  these  results  to  spot 
beam  transmission  and  other  types  of  clouds  would  be  desirable  in  view  of  the  potential  impact  on  system 
availability. 

4,2  UPLINK  DATA  REDUCTION 

The  experimental  procedures,  descriptions,  and  scenarios  for  the  uplink  measurements  are  outlined 
in  Section  3.  The  equipment  used  and  the  calibration  employed  are  described  in  Volume  II,  Section  3.  The 
goal  of  the  uplink  measurements  was  to  determine,  as  best  as  possible,  those  purumeters  which  arc  critical  in 
determining  uplink  performance.  Specifically,  there  were  three  parameters  which  received  the  most  atten- 
tion: overull  link  loss,  beumwidth  (antenna  gain),  und  beam  direction.  Furthermore,  it  was  desired  that 
models  be  developed  which  could  adequately  describe  these  parameters,  The  model  used  in  this  report  is 
the  one  developed  in  Appendix  A.  All  the  data  presented  are  in  one  format;  link  loss  vs.  zenith  angle. 

These  represent  a slice  of  the  transmitted  pattern  as  seen  by  a receiver  passing  through  the  beam  toward  the 
source.  The  program  used  to  reduce  the  uplink  data  is  described  in  Appendix  E.  Because  of  equipment  dif- 
ficulties, aircraft  failures,  und  u narrow  operational  window,  the  only  data  deemed  acceptable  were  obtained 
on  July  21  and  22,  1975,  We  also  point  out  that  for  logistical  reasons,  the  uplink  measurements  were 
always  given  first  priority  over  the  downlink  measurements,  which  contributed  to  some  spottiness  in  the 
latter. 

In  figures  4-1 8A  through  F,  data  taken  with  the  luser  pointing  in  the  zenith  direction  are  displayed; 
e.g.  zenith  ungle  is  0°.  There  are  several  items  that  should  be  pointed  out  with  regard  to  the  data.  First,  the 
dynamic  range  of  the  receiver  was  only  10  dB.  so  thut  the  dutu  went  from  the  noise  flow  to  saturation  very 
quickly.  This  required  some  adjustment  on  each  pass;  for  example,  figures  4-1 HB  and  C were  repeated  pusses 
taken  with  different  gain  settings.  Notice  also  that  at  the  higher  zenith  angles  the  noise  floor  sturts  to  in- 
crease. This  is  due  to  the  secant  squared  correction  that  was  used  to  account  for  the  difference  in  path 
lengths  at  the  various  zenith  angles  to  a constant  altitude  aircraft.  This  effect  can  also  be  observed  when  the 
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Figure  4-18A.  Radiance  profile  through  angle  of  aircraft  (Run  No.  16, 22  July  1975). 
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Figure  4-I8C.  Radiance  profile  through  angle  of  aircraft  (Run  No,  19, 22  July  1975). 
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Figure  4-18D.  Radiance  profile  through  angle  of  aircraft  (Run  No.  20, 22  July  197S). 
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Figure  4*lBli.  Radiance  profile  through  angle  of  aircraft  (Run  No.  21, 22  July  1975). 
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Figure  4-I8F.  Radiance  profile  through  angle  of  aircraft  (Run  No.  22, 22  July  1975). 
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data  are  saturated.  On  most  data  plots  there  is  a constant  line  going  through  zero  degrees  zenith,  Because 
the  aircraft  could  not  always  pass  directly  overhead,  a minimum  zenith  angle  resulted.  Since  a sign  reversal 
occurs  in  the  gimbul  readings,  a plus/minus  indication  occurs  on  adjacent  pulses,  Along  the  abscissa,  an 
indication  of  when  the  receiver  was  in  track  has  been  noted.  A deviation  from  the  abscissa  indicates  track. 
The  filtering  of  the  data  is  described  in  Appendix  E.  Filtering  of  approximately  one  second  was  used  lor  a 
few  reasons.  First,  this  was  the  shortest  integration  time  that  would  eliminate  the  grass  in  the  datu  due  to 
noise  without  altering  the  results.  Second,  the  aircraft  was  traveling  at  80  mph  which  is  approximately  1 17 
feet  per  second.  At  3,000  feet  this  is  2.2°,  while  at  2,000  feet  this  is  3,3°.  Since  the  beamwldth  cf  the  sourer 
was  2°,  this  is  approximately  the  maximum  resolution  inherent  in  the  experiment.  Finally,  the  rrvdts  of 
the  model  were  overluid  using  the  data  taken  by  Scripps. 

In  figures  4-18  A through  F,  the  model  ranges  from  several  dB  too  high  to  several  dB  too  low,  and  is 
consistent  with  the  downlink  data.  There  also  appears  to  be  some  conservatism  at  the  deeper  depth,  which 
is  appropriate  for  system  design.  However,  there  is  some  uncertainty,  that  is  more  difficult  to  explain. 

Notice  that  figures  4-18  B and  C which  were  taken  within  15  minutes  of  each  other  and  represent  the  same 
scenario  arc  10  dB  different,  with  the  model  falling  halfway  in  between.  One  could  possibly  envision  some 
sudden  change  in  the  environment  to  expluln  this.  Whut  is  more  plausible,  however,  is  to  look  for  other 
causes,  For  example,  the  sensitivity  of  the  receiver  was  changed  between  the  two  runs  by  inserting  a neutral 
density  filter.  Although  it  was  accounted  for  in  the  calibration,  this  might  be  suspect.  Also,  the  percentage 
of  track  was  different  in  the  two  runs.  And  finally,  there  could  be  some  dynamic  effects  caused  by  the 
ocean  surface  that  might  have  caused  the  difference.  What  can  be  concluded,  however,  is  qualitative  concur- 
rence. This  same  qualitative  concurrence  is  also  maintained  us  the  angle  of  the  laser  source  is  varied. 

In  figures  4-19  A through  E,  the  laser  angle  Is  changed  to  12.5°  off  the  vertical,  (The  refracted  angle 
predicted  from  Snell’s  law  is  also  plotted.)  In  figures  4-20  \ through  D,  tills  angle  Is  Increased  to  32.5°  and 
in  figures  4-21  A through  J,  to  42.5°.  In  every  case  where  two  scenarios  were  repeated,  there  were  several 
dB  variation,  with  the  model  residing  In  between.  Furthermore,  the  data  clearly  Indicate  a spreading  of  the 
beam  away  from  the  Snell's  angle  and  toward  the  zenith.  This  has  been  predicted  by  the  model  and  Is  In 
qualitative  agreement  with  the  data, 

The  above  datu  have  been  selected  so  thut  only  the  better  runs  are  presented.  Tables  4-4  A through 
D represent  complete  listings  of  all  the  data  with  some  accompanying  comments.  Approximately  50  hours 
of  aircraft  time  were  employed,  The  data  taken  on  July  24  were  lost  due  to  a tape  recorder  malfunction. 

4.3  CONCLUSIONS  AND  RECOMMENDATIONS 

Wliile  it  cannot  be  said  that  all  the  goals  or  the  experiment  were  reached  In  a quantitative  manner, 
the  experiment  was  nevertheless  highly  successful,  Use  of  the  radiative  transport  theory  as  a principal  tool 
In  predicting  system  performance  in  the  ocean  environment  has  been  clearly  established,  Although  some 
approximations  were  used  in  the  OPSATCOM  application,  It  Is  also  clear  that  a more  gloi  a)  application  Is 
well  within  the  state  of  existing  knowledge.  The  basic  parameters  on  which  the  radiath  transport  theory 
are  based  are  well  known.  However,  the  mechanism  for  extracting  these  parameters  Is  still  an  art,  and  some 
advances  in  this  direction  were  also  made  (Appendix  D).  Some  of  the  regression  curves  derived  are  shown 
In  figures  4-22  through  4-24.  The  most  interesting  one  from  a syjtem  point  of  view  is  figure  4-24,  This  Is 
the  regression  of  against  s.  Notice  thut  in  the  clearer  waters,?^  is  larger  than  In  the  turbid  waters.  This 
of  course  implies  that  In  the  latter  case  the  particulates  are  large  and  concentrate  the  scattering  in  a forward 
direction.  This  In  turn  will  somewhat  offset  the  deleterious  effects  of  a large  value  of  s upon  system  per- 
formance, On  the  other  hand,  In  clear  water  wo  must  recognize  that  the  medium  is  less  forward  scattering, 
which  will  degrade  system  performance  to  some  degree. 
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Figure  4-I9A,  Radiance  profile  through  angle  of  aircraft  (Run  No.  12, 21  July  197S). 
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Figure  4-19B,  Radiance  profile  throu^i  angle  of  aircraft  (Run  No.  1,21  July  1975). 
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DEPTH  « 18.3  METERS  ( 60  FEET  ) 

R/C  RLTJTUOE  * 914  METERS  f 3000  FEET  ) 
LRSER  RNGLE  * 12. 5 DECREES 
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Figure  4-19E.  Radiance  profile  through  ingle  of  aircraft  (Run  No.  6, 21  July  1975). 
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OEPTH  - 9.1  METERS  ( 30  FEET  ) 

A/C  ALTITUDE  - 914  METERS  ( 3000  FEET  } 
LASER  ANGLE  - 32,5  DEGREES 
PERCENT  TRACKING  -93.8  „ 

PEAK  -2.75  X 10-7  METER-2  AT  54.9  DEGREES 
HALF- POWER  REAMWIOTH  ■ 20,1  DEGREES 


Figure  4-20A,  Radiance  profile  through  angle  of  aircraft  (Run  No.  14,21  July  1975), 
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ZENITH  RNGLE,  DEGREES 

Figure  4-20B.  Radiance  profile  through  angle  of  aircraft  (Run  No.  16, 21  July  1975). 
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DEPTH  * 15.2  METERS  ( 50  FEET  ) 
ft/C  ALTITUDE  = 914  METERS  ( 3000  FEET  ) 
LASER  fiNGLE  « 32.5  DECREES 
PERCENT  TRACKING  « 95.0  , 

PEAK  * 1.35  X 10°  METER*2  AT  37.8  DEGREES 
HALF-POWER  BEAMWIDTH  - 19.2  DEGREES 


Figure  4-20C.  Radiance  profile  through  angle  of  aircraft  (Run  No.  17, 21  July  1975). 
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DEPTH  = 9.1  METERS  ( 30  FEET  ) 

A/C  ALTITUDE  * 610  METERS  ( 2000  FEET  ) 
LASER  ANGLE  = 42.5  DEGREES 
PERCENT  TRACKING  « 74.6  „ 

PEAK  * 1.35  X 10"7  METER'2  AT  56.5  DEGREES 
HALF-POWER  BEAMWIOTH  * 11.5  DEGREES 


Figure  4-21A.  Radiance  profile  through  angle  of  aircraft  (Run  No.  12, 22  July  1975). 
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DEPTH  « 9.1  METERS  [ 30  FEET  1 
R/C  RL71TUDE  = S10  METERS  ( 2000  FEET  ) 
LRSER  RNGLE  * 42.5  DEGREES 
PERCENT  TRRCKING  - 96.7  „ 

PERK  = 6.3  X 10‘7  METER“2  RT  36.5  DEGREES 
HRLF-POWER  BERMWIDTH  « 1.8  DEGREES 
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Figure  4-2 IB,  Radiance  profile  through  angle  of  aircraft  (Run  No.  1,22.  July  197J). 
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DEPTH  * 12.2  METERS  ( 40  FEET  ) 

R/C  ALTITUDE  * 610  METERS  ( 2000  FEET  ) 
LASER  ANGLE  * 42.5  DEGREES 
PERCENT  TRACKING  - 91.3  _ 

PEAK  « 3.72  x 10"7  METER"2  AT  65.4  OEGREES 
HALF-POWER  8EAMW3DTH  - 1.6  DEGREES 


Figure  4-21C.  Radiance  profile  through  angle  of  aircraft  (Run  No,  2 22  July  1*>75). 
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OEPTH  - 12.2  METERS  ( 40  FEET  ) 

A/C  ALTITUDE  = S10  METERS  ( 2000  FEET  ) 
LASER  ANGLE  = «.s  DEGREES 
PERCENT  TRACKING  „=  77.4 

PEAK  * 7.59  X 10"6  METER"2  AT  58.4  DEGREES 
HALF-POWER  8EAMWIDTH  = 19. 4 DEGREES 


Figure  4-2 1 D.  Radiance  profile  through  angle  of  aircraft  (Run  No.  1 1 , 22  July  I '>75). 
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DEPTH  * 15.2  METERS  f 50  FEET  ) 

P/C  RLTITUDE  = 510  METERS  f 2000  FEET  ) 
LPSER  ANGLE  = 42.5  DEGREES 


Figure  4-2 1 E.  Radiance  profile  through  angle  of  aircraft  (Run  No,  3,  22  July  1975). 
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DEPTH  * 15.2  METERS  ( 50  FEET  ) 

A/C  ALTITUDE  « 610  METERS  ( 2000  FEET  ) 
LASER  ANGLE  « 42.5  DEGREES 


Flguri  4*2 1 F.  Radiance  profile  through  angle  of  aircraft  (Run  No.  10, 22  July  1975). 
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DEPTH  « 21.3  METERS  ( 70  FEET  ) 

R/C  RLTITUDE  * 610  METERS  f 2000  FEET  ) 
LftSER  RNGLE  - 42.5  DEGREES 
PERCENT  TRRCK1NG  = 89.9  9 

g PERK  • 1.32  X 10“7  METER"2  RT  42. 2 DEGREES 

5 HRLF-P0NER  8ERMWJ0TH  « 18.1  DEGREES 
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Figure  4>2  111.  Radiance  profile  through  angle  of  aircraft  (Run  No.  $,  11  July  1975). 
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DEPTH  * 24.4  METERS  ( 80  FEET  ) 
ft/C  ALTITUDE  * 610  METERS  ( 2000  FEET  ) 
LRSER  ftNGLE  * 42.5  DEGREES 
PERCENT  TRACKING  - 86.5  . 

PERK  * 9.77  X 10"8  METER"2  RT  43.1  OEGREES 
HALF-POWER  8EAMWI0TH  = 36.5  DEGREES 
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Figure  4-211.  Reditnce  profile  through  angle  of  aircraft  (Run  No.  6, 22  July  1975), 
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DEPTH  = 27.4  METERS  ( 90  FEET  ) 

/VC  ALTITUDE  =610  METERS  ( 2000  FEET  1 
LASER  ANGLE  = 42.5  DEGREES 


TABLE  4-4 A.  PBY  DYE  LASER 
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Had  trouble  starting  APU  an  this  trip.  After  connector  rattlrog  and  relay  Upping,  it  west.  Intermittently  dropped  out  three  tunes  and  finally  held. 


TABLE  4-4C.  PBY  DYE  LASER  TRACKS. 
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With  the  confirmation  of  radiative  transport  as  the  underlying  theory  for  these  system  applications 
them  are  some  suggestions  that  can  be  made  for  future  work: 

a.  General  Geanup.  This  is  basically  a firming  up  of  the  general  theory  to  include  all  scenarios 
that  the  Navy  might  envision. 

b.  Monte  Carlo  Techniques.  These  are  primarily  an  augmentation  of  the  general  theory  to  include 
thoae  cases  that  are  not  conducive  to  analytic  computation. 

c.  Cloud  Propagation.  This  it  an  extension  of  (a)  and  (b)  that  includes  any  scenarios  in  the  total 
environment  that  are  of  interest  to  the  Navy. 

d.  Model  Range.  Establishment  of  the  theory  automatically  implies  an  understanding  of  scale 
change.  This  means  that  system  scenarios  can  be  studied  in  a suitably  designed  test  tank  similar  to  the 
manner  that  antennas  are  designed  on  an  antenna  range. 

e.  Experimental  Test  Site.  In  much  the  same  manner  that  a test  tank  can  scale,  so  too,  can  a test 
site.  Having  such  a facility  as  a convenient  adjunct  to  system-related  design  activities  would  allow  for  total 
Interaction  between  theory,  experimentation,  and  testing. 
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APPENDIX  A 

OPTICA  I C MMU  NIC  ATION8  BETWEEN  UNDERWATER  AND 
AhO\  c.  SURFACE  (SATELLifS''.’  TERMINALS 

SHERMAN  KARP 


A buna  - A multipl*  wittering  modal  it  uird  arvO  extended  to 
vhuictnri**  lh»  channel  between  undtrwstn  and  airborne  tiatellttn) 
terminals  at  optical  frequencies,  Tht  effect!  of  th*  air/ wi  Inicrfse*  an 
also  iiuluilrtl  with  approximate  solutions  tccuratr  for  elevation  angWl 
ibovt  45* . Tht  results  u«  presented  In  Mims  of  i radiance  function 
which  It  related  lo  ih«  transform  of  Iht  spatial  covastanc#  fimvilon 
(mutual  cohrrtnc*  function).  Th*  print uy  Iohm  an  shown  to  bt  t 
result  of  tht  wiiti  absorption  coefficient  mi  not  Iht  extinction  co* 
•fflcltnt.  Tht  wittering  lowai  can  bt  laolsted  fiom  Iht  abioipllon 
loattt  and  for  certain  cates,  whtrt  Iht  receivet  la  Imbtddtd  In  tht 
scattering  modlum,  can  bt  completely  recovered,  New  compontnu  may 
bt  required  lo  achieve  this  performance.  The  afftcU  of  ocean  roufhnaia 
art  ahown  lo  have  a minimal  effect  upon  Iht  aubiurtact  reception 
while  cautlns  potalble  beam  iteerins  of  tubtutfact  tianamlwion.  Al- 
though lubitantlal  loiaei  art  experienced,  duplex  operation  can  bt 
achieved  it  modeii  data  rater, 


I.  INTRODUCTION 

THE  ACCEPTANCE  of  optical  communication!  for  uie  In 
operational  ayatemi  hat  been  severely  bumpered  by  our 
Inublllt)  to  udequulcly  compensate  lot  channel  effects  Induced 
by  the  environment,  C'uniequetuly . most  ll  not  ull  of  the  pro- 
jected system  gains  are  quickly  nullified  when  rudimentary 
meusuras  of  system  margin  arc  added  to  the  link  budgets  to 
account  for  these  effects,  ll  It,  therefore,  extremely  Important 
that  environmental  effects  be  accurately  accounted  lor,  und 


platform  In  a multiple-scattering,  absorbing,  and  dispersive 
medium,  through  a random  surface  end  Into  a medium  with 
a different  Index  which  it  nunscittering  and  nondltperslve. 
These  two  problems  are  nonreciprocal.  Thus,  It  It  necessary  to 
decompose  them  Into  their  fundamental  elements  and  to 
Individually  Identify  and  characterize  the  contributing  factors. 
To  this  end  this  discussion  It  divided  Into  four  purta.  The  first 
part  involves  the  actual  propagation  effects  encountered  while 
traversing  a multiple  scattering  medium.  The  three  system 
parameters  which  can  be  Identified  are  the  attenuation,  the 
beam  spreading  and  the  apparent  source  size.  These  in  turn  are 
related  to  the  absorption  coefficient,  the  scattering  coef- 
ficient, and  the  volume  scattering  function.  The  second  part 
uddresses  the  problem  of  transmitting  through  a rundom  sur- 
face characterized  by  a slope  distribution,  Tht  effect  on 
scintillation  wilt  be  discussed  In  addition  to  beam  pointing 
und  beam  broadening.  The  third  und  fourth  parts  will  address 
link  calculations  from  the  satellite  plutform  to  the  submerged 
platform  and  from  the  submerged  platform  to  the  satellite 
platform,  respectively.  It  is  estimated  that  the  model  presented 
can  ba  verified  to  within  several  decibels  over  most  operational 
scenarios  envisioned. 

II.  THE  UNDERWATER  CHANNEL 


systems  designed  to  best  exploit  these  channels  In  u most 
advuntugeoui  manner,  The  most  difficult  channel  that  the 
optical  communications  engineer  has  to  deal  with  is  the  mul- 
tiple scattering  channel.  Such  a channel  exists  when  propagat- 
ing through  clouds,  log  , water,  etc,  1 1 )-|.t) . In  this  paper  we 
will  extend  the  model,  which  has  been  Independently 
developed  by  llcggestad  |4|  and  Arnuslt  |5|  for  multiple 
scattering  mcdlu,  and  upply  It  to  compute  the  effects  we 
would  encounter  while  traversing  a satellite  to  underwater 
channel,  In  doing  so.  we  will  try  to  validate  the  use  und 
Interpretation  of  the  model  by  applying  ll  to  experimental 
data, 

In  Its  most  general  form,  the  problem  of  optlcul  com- 
munications between  u satellite  and  a submerged  platform  cun 
be  described  as:  I)  u problem  In  communications  from  u 
platform  In  u nonscaiterlng,  nondltperslve  environment, 
through  u random  surface  and  Into  a medium  with  u different 
index,  which  Is  multiple-scattering,  absorbing,  und  dispersive: 
and  conversely.  2)  a problem  In  communications  from  u 


Over  the  punt  two  decades  there  has  been  an  Interest  In 
understanding  the  behavior  of  light  while  propagating  through 
water,  With  the  advent  of  the  laser  this  Interest  Intensified 
when  viewed  In  the  context  of  operational  equipment. 
Although  there  have  been  numerous  measurements  (6|,  |7) 
and  many  empirical  curves  derived  to  lit  tht  data  |7|  the 
latter  ure  of  limited  use  for  extrapolating  system  petformunco, 
For  this  paper  we  will  use  a model  which  hus  been  developed 
Independently  by  two  separate  authors  |4J,  |5).  While  this 
model  Is  derived  for  smull  ungic  forward  scattering  It  appears 
to  be  fulrly  uccurule  out  to  *45“  providing  the  optical  thick- 
ness is  not  too  (urge.  Fortunately  the  runye  of  validity  Is  within 
the  operational  runges  envisioned.  The  mudel  describes  the 
radiance  transfer  while  (ruverstng  ihe  multiple  scattering  region. 
Tills  region  Is  characterized  by  three  variables. 

II  Absorption:  The  absorption  coefficient  of  the  modlum, 
u,  Is  the  umounl  of  energy  absorbed  by  the  medium  per  unit 
length  of  propagation.  This  loss  Is  uttenuutlon  und  goes 
directly  into  heutlng  und  otltet  Irreversible  processes. 


Puper  upprovsd  by  the  Associate  Editor  lot  Space  Communication 
or  th*  IEEE  Communications  Society  lor  publication  without  ore! 
presentation.  Manuscript  received  March  11,  1975:  revised  July  35, 
1973, 
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2)  Scattering  Function  Ft 0V  Multiple  scattering  medtu  are 
characterized  by  icatterlng  centers  which,  In  the  case  of  ocean 
water,  appear  to  be  both  from  plankton  and  from  molecular 
scattering.  The  volume  scattering  function  Is  defined  as  the 


secondary  radiation  puttern  created  by  a plane  wave  Wuverslng 
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a small  enough  volume  to  that  only  tingle  scattering  occurs. 
This  represents  the  iversge  scattering  distribution  of  all  the 
scattering  centers.  There  does  not  appear  to  be  a .treat  deal  .if 
variation  in  the  general  shape  of  /'(0)  although  the  average 
width  does  change 

Si  Jtii'riny  Coc/jUn'M  It  we  normalize  /'((M.  then 
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Is  the  normalized  version,  s Is  the  Mattering  coefficient  with 
i“l  interpreted  as  the  average  distance  between  scatterings. 
Ainush  assumed  a form  lor/10)  as 
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lleggestad.  on  the  other  hand,  defines  a modified  variance 
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and  by  equating  0*  with  I /A*  the  two  models  aje  Identical, 
This  is  true  for  large  6.  We  will  use  the  notation  07,  It  Is  also 
common  to  define  an  extinction  coefficient  o,  defined  by 
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Thus,  to  completely  characterize  the  environment  It  would 
be  necessary  to  have  measuring  equipment  for  a,  i,  and  /(0). 
An  alternate  procedure,  und  less  desirable,  would  he  to  mea- 
sure two  parameters  und  scale  the  measurements  to  the  third. 
Althuugh  feasible,  this  would  assume  the  validity  of  a model 
und  the  confidence  to  extrapolate  from  It.  With  llicsn  param- 
eters In  mind  wc  will  now  present  the  model. 

For  convenience  we  assume  that  the  transmitted  beam  Is 
Gaussian  und  hui  the  form 
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Some  explanation  of  the  Interpretation  of  (7)  Is  now 
warmntud.  First  notice  lliul  If  we  hud  u receiver  ut  point 
tr.z)  und  added  equally  the  contributions  coming  from  all 
ungias,  we  could  Integrate  over  the  vurluhles  6„  6$  and  obtain 

the  result1 
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That  is.  it  huR  a Gaussian  distribution  both  In  In  spatial  cross 
section  und  Its  iuy  direction.  Next  we  assume  the  geometry  in 
l:lg,  I.  whore  the  source  Is  at  (0,0),  the  observer  Is  u\  |r,:l  and, 
as  we  will  see,  the  apparent  source  Is  at  (0,zo>-  lb  terms  of  the 
observation  point  (r.rl  we  have  us  the  triinifer  In  Intensity 
/t«,D 


/i»,D 


Thus  the  total  energy  has  a distribution  which  is  Gaussian  in 
the  z-ptane,  centered  at  r » 0,  und  Is  a result  of  scattering.  The 
itumlurd  deviation  of  this  spread  it{R\*l2)vi  ■ (s»,sffa/6)U  + 
)(/  + I'))  1 /'*.  If  In  addition  we  could  collect  all  the  scattered 
rudlullon  in  the  z-plane  (a  largo  collector)  we  would  Integrate 
uverz  und  obtain 


/■*  Ja 


/(r)  Jr  • 


02) 


und  identify  this  at  an  Irretrievable  loss  which  we  see  Is  due  to 
ubaorptlon. 

1 Mire  If  the  ray  Is  coining  from  the  direction  • *»  then  the 

reveivet  Is  pointed  in  the  ( >)  - ( -fi,,  os)  direction,  Hence  then  Is 
only  a sign  difference  between  the  two. 
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' • -uypoje  we  observe  the  source  at  the  point  (i>,  x)asa 
function  of  angle,  Fig.  2. 

Notice  I*  'he  unit  vector  In  the  direction  O0,  the  angle 
out  of  the  r-;-plane,  and  y.  Is  the  unit  vector  representing  the 
uiigulur  tilt  up  from  the  x-uxls  0r  in  u plane  described  by  r.i. 
rims  fur  any  r.  the  maximum  always  occuis  in  the  er-planc 
U\,  - 0)  at  a till  angle  of  l)m  Alternatively,  for  a lived  tilt 
angle  (I,  lire  maximum  occurs  when  the  receiver  is  off  the  axis 
a distance  /-m.  The  net  result  of  Inn  It  ii.'cipieiuthms  is  that  the 
source  uppeurs  to  be  locuted  at  the  point  x0,  Fig.  I,  where 

Furthermore,  the  source  will  huve  an  uppareni  extent  (sine)  in 
diumetei  (twice  the  standard  deviation)  of 


(15) 

In  radiant  (field  of  view).  Consequently,  any  system  should 
account  for  the  spatial  filtering  that  may  occur  when  optical 
elements  are  used. 

Finally  we  can  Identify 


as  the  ratio  of  the  initial  beam  cross  section  to  the  cross  section 
at  (r,x)  which  should  be  much  loss  than  one  and 


us  the  rutlu  of  the  Initial  beum  spread  to  the  beam  spread  at 
(r.ij  which  should  ulso  be  much  less  than  one.  Thus  we  cun  set 
/ * V « 0 when  collimated  beams  are  used. 

Strictly  speaking,  the  model  used  here  is  only  valid  for 
small  angle  forward  realtor,  This  Is  true  because  In  the  deriva- 
tion, the  approximations  ain  5—0  and  cos  0 - I are  used, 
llowovor,  these  approximations  are  only  off  by  10-20  percent 
at  30-40°  und  hence  the  model  should  degrude  gracefully  at 
lurger  angles.  Some  modification  has  to  he  made,  however,  to 
use  this  at  large  unglei.  This  la  due  to  the  fact  that  the  absorp- 
tion and  scattering  pa t ha  are  longer  by  the  factor  :|sec  0 I ] 
ai  the  unglo  0.  Thls_can  be  easily  accounted  for  by  changing  x 
in  i sec  0 s/*3~+  ^ wherever  x occurs,  Then  we  will  Inter- 

pret (7)  to  be  the  transfer  In  Intensity  from  the  source  to  a 
sphere  of  radius  x,  With  the  latter  Interpretation  in  mind  we 
will  now  show  the  Justification  uf  using  this  model  und  then 
point  out  the  remaining  verification  naadad. 

Consldar  the  geometry  In  Fig,  3.  A collimated  aourca  amlta 
radiation  along  the  x-axii.  The  medium  ti  characterized  by  the 


l lg.2.  Source  si  ti  function  of  angle. 
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Fig.  3,  Geometry  used  by  Duntley. 


ratio  a/a.  Since  a “ x + rr,  s/a  = (a/a)  - I , The  unit  of  length 
t»  N * <u  extinction  lengths,  Thus.  N extinction  lengths  cor- 
respond by  the  relationship 


N » as  ■■ 


a / a\ 

a W 


•biorpllon 


to  |AV(a/tf)J  absorption  lengths  and  since 

[a/a  --  1 1 

A/,!.,,  * k « ((*  + a)  ~ a\:  *»  a:  -a:  -•»  ,V  ! (Id) 

a/a 

to  ((a/tf)  - 1 1 Ha/a)  scattering  lengths. 

With  this  geometry  und  purunttieilzetion.  Duntley  |7|  has 
made  extensive  measurements  of  the  power  collected  us  u 
function  of  N for  various  values  of  4'  Two  tepresenlailve 
samples  urc  shown  In  Fig.  4,  For  this  case,  (II)  Integrates  to 

f -_£*)] 

['  "M’(,.'TSTi;)] 

•'"(  «,) 

f - V*  1 

, | -exp  ; (20) 

■ w a/a  ~ I 


L 1 a/a  , _ 


h a/a  - 
(2/3kf,.V  


As  pointed  out  by  Duntley  and  as  observed  In  (12), 
for  large  values  of  ♦ we  would  expect  the  relationship 


‘.DM Ml. I Ill'll  ANl  | 


cot  I it.  ion  Mil  it  'lAUWm- 

rig.  5.  Duntloy  m*uiur«m«nli  or  anomalous  modal  behavior. 

Atph«/U  «<  2.49. 

v 'V»biorpHon  to  hold,  However.  tliU  Is  not  happening  even 
though  the  curve  saturates  lor  'I'  > 60°.  We  cun  compensate 
with  our  model  by  reculling  that  the  effective  source  broudens 
and  ihll'ti  location  at  N Inureatsa.  Duntloy  wai  able  to  observe 
and  meaiure  the  former  phenomenon  although  he  was  not 
able  to  explain  It  (see  Fig.  5).  In  the  actual  measurements  an 
Integrating  sphere  was  used  In  the  collecting  optics.  This  has  a 
spatial  response  of  cos  0.  Consequently  ut  lurge  values  or  iV  one 
would  expect  to  start  to  observe  sputlul  filtering  of  the  source. 
This  Is  precisely  what  we  see  for  N greuter  than  8.  This  was 
corrected  for  as  follows.  It  wus  usiumed  that  the  difference 
between  the  e~N»biorption  line  and  the  W » 100°  curve  was 
due  to  spatial  filtering.  Therefore,  ut  every  value  of  y this 
difference  was  added  to  each  of  the  curves  (on  a log  scale)  and 
replotted.  Thn  model  was  then  calibrated  ut  the  largest  value 
of  N and  the  smallest  value  of  sh,  where  It  should  be  most 
accurate,  und  P was  calculated.  The  model  wus  then  plotted 
on  the  revised  curves,  Fig,  (t.  Although  the  agreement  Is  not 
perfect,  It  is  remarkably  close.  The  vulldlty  of  this  culibrutlon 
should  be  checked  at  some  point  by  comparing  the  results  of 
an  Integrating  sphere  with  those  of  u hemispherical  coverage 
lens  (flshcyc), 

III,  THU.  AIR/SHA  INTERFACB 


""M "|iNi,>,i.iii'.>i  In  this  section  a geometric  optics  model  will  be  developed  ' 

(hi  to  determine  the  effect  of  surface  Irregulatlties  on  heum 

1 ip.  4.  Total  power  mi  spherical  eap.  (at  Alpha/, ■!  2.49.  spreading,  pointing  und  scintillation  when  traversing  the  J ' 

(hi  Alpha//)  4,  boundary,  To  do  so,  we  consider  the  following  model  of  an  i 
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In,  We  do  ilili  In  the  following  manner,  Given  u sample  R from 
the  ut  of  poiilble  ilopei,  ■>'  U well  dellned,  That  is.  the  prob- 
ability of  V conditioned  upon  7 und  R it 

p(y'/R,y)  * 4 j\’  - ^iln-'1  *ln  (7  + R)j 


(2J) 


We  arrive  at  the  angular  dlitrlbuilon  of  7'.  given  7.  by  averag- 
ing over  the  varluble  R.  Thus, 


p{y‘h)mJ  ilRpn(R)p(y'ly,R) 

•n /”  <JRpH(R) 

•6^7'  - Itn-1  ^iin(7  + «)^ 


(24) 


where  PuiR)  ts  the  probability  density  of  R. 

Rigorously,  this  It  merely  the  change  of  variables  in  the 
density  ) from  R to  7'  using 


R « tan-*1 


rt|.  6. 


Revised  model  curves.  « t~Wli.it),  tb)e  **  - 

t~Wli), 


-,sln  7 - sin  7' 


cos  7 ~ — cos  7 
ft 


1 R(y'.y)> 


(2-M 
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L 


Thus, 


iJfUt'.y) 


mH 

• 

*PH 

tan 

, sin  ■> 
M 

sin  ) ' 

cos  7'  - 

n 

...  C0Sj  7 

» 

n 

-;COI(7  - 7')  - 1 


I + 


(5H 


COB  (7  71) 


Knowing  piy'/y)  we  can  compute  the  uveruge  spreuding  uml 
offut  of  1 ray  incident  at  the  angle  7.  Tltii  becomes 

average  offset  « j y'plyiy)  dy  ■»  y 

“ j | •in-1  £ ~ »iu  (7  + A )J  - H J />*tAl  i/A 
- Jain*  1 sin  17  4 A|jpH(AMA  — 7i. 

De  lining  7*  us 

71  - J y'^y'huty' 

” jjslir1  £ ~iln|7  + rtl^j » A J pH[R)dR  t-Hi 


the  rm»  spread  becomes 

A.  <7* -?*)«* 


t:i»i 


There  are  tome  practical  llmlluiiotn  to  these  results  which 
require  modification,  Tig.  8, 

A ray  of  light  with  zenith  angle  7 will  never  inteicept  a 
wave  whose  slope  Is  greaici  than  it 1 2 7 because  of  wave 

obscuration  However,  the  my  will  itill  peneiiau  the  inter- 
face with  probability  one.  Consequently,  the  limits  ot  integra- 
tion for  A arc  cel  at  | v!2.  nl2  y|  uml  Ihe  densin  />#lfl) 
should  be  modified  to  that  of 


/. 


era  y 

pH\R)dR 


Pk(A’),  -nl2<R<r2  y 

0,  elsewhere,  l.Ul) 


Tlf  results  in  would  then  he  modlllcd  In  lepluclng 

PtiUi)  with  Ph\R)  In  general,  the  results  prescuied  can  be 
simplified  by  only  considering  lliosc  values  of  (7  l- A) <45°. 
This  corresponds  to  the  major  operational  icqulromenu  and 


(20)  gives  good  engineering  insight  iniu  the  bchuvlor  of  a ray  going 
through  the  air/iea  interface.  Fur  this  case 


P\y'h)-Pn 


« 

-7-7 

n 

1 

■ ■!■■■ 

H 

| 

n 

1 __ 

L «'  J 

1 r 

If 

>("?)* 


n 

I -- 
n 


var  (A) . 
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Notice  tlial  I l 11/111  |a  < I for  the  air/sea  Interface  (Index  of 
wutci  - I ..Cl,  Index  of  air  * I)  and  consequently  tint  ray 
spreading  is  appreciably  leas  than  die  slope  spreading  of  the 
ocean.  In  addition  the  surface  adds  the  contribution  (I  - 
rt/’r'V?  to  the  normal  bending  due  to  Snell's  law. 

In  all  cases,  the  modal  used  represents  nn  optical  beam  of 
zero  cross  section  and  zero  divergence.  The  exiling  beam  also 
has  zero  cross  section  and  zero  divergence  hut  Is  being  steered 
by  the  roughness  of  tire  surface-  If  this  surface  is  the  ocean. 
lAy'ly)  represents  the  average  time  history  ol  the  beam  direc- 
tion 7',  with  p(>’, '7)1/7 ' the  probability  dial  li  is  pointing  with- 
in a dy'  Interval  ol  the  y'  direction  ui  any  Instuid  I11  dine.  This 
apparent  beam  wander  would  cause  severe  scintillation  m un 
operating  system  As  the  cross  section  of  the  beotp  increases 
the  refracting  surface  can  no  longer  be  considered  locally  flat 
and  different  portiuns  of  die  beam  are  refracted  at  different 
atiglci  Consequently.  II  would  be  possible  to  uveruge  out  the 
ocum  wander  In  the  direction  of  die  mean  7 by  spreading  the 
beam  over  a linger  portion  of  the  surface.  If  the  area  of  the 
beam  is  A.  and  the  correlation  length  of  the  surface  statistics  is 
then  there  are  approximately  ,d/(ir(/,/2)la  Identically  dls- 
idbuicd  independent  paths  slmllur  to  diversity  paths.  If  we 
lint  her  assume  a depth  : such  that  the  beam  cross  section  Is 
greater  than  .1. 


.1 


>.-». 


1.121 


then  nil  of  the  paths  will  overlap  at  die  icccivet.  This  can  be 
analyzed  In  die  following  mumicr.  I'iisl  wc  notice  dial  the 
probability  of  hiving  the  beam  within  an  nils  dcv.allon  about 
the  mean  is 
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/">  +A 

I ph'hMy'  (33) 

-V  u 

win,  li  Ini  the  Gaussian  density  becomes  0.<>s.  Ihux.  even  il  we 
Ii.kI  mi  lime  variations  the  beam  would  only  he  within  a devil- 
lion  oi  the  mean  <>K  percent  of  ilie  lime  Now  suppose  we  pick 
N independent,  identically  distributed,  paths  to  the  receiver 
and  transmitted  (l/A/’lilt  of  the  power  p,  in  each  path.  Since 
the  paths  are  identically  distributed,  the  average  direction  of 
the  sum  is  still  y.  Now,  however,  the  variance  of  ihe  sum 
becomes  A*/M  or  a standard  deviation  for  the  sum  of  A/y/fi' 
about  the  mean  y.  If,  for  example,  we  set  ,V  = 2$.  und  assume 
the  central  limit  is  approximately  valid,  Ihe  probability  that 
the  beam  is  within  Is  now  0.999994  Since  the  correlation 
length  Is  approximately  the  separation  between  Independent 
spatial  Nyqulst  samples  we  sec  that 


and  can  be  used  accordingly.  Furthermore,  it  can  be  shown 
that  the  scintillation  will  reduce  the  average  xignal-to-rulse 
ratio  by  the  factor 


tc  i 

I + S*IN 

A verification  of  these  results  and  the  relationship  in  (32) 
would  be  warranted. 

Finally,  we  can  interpret  the  function  p[y'ly\  as  a beam 
spreading  factor.  Thus,  if  we  have  a propagating  beam  of  ihe 
form  ,r\0,r)  in  (7),  or/fy.M  then  the  output  beam  after  travers- 
ing the  surface  will  be 

/li '.»•)-  f>iyp(y'lyV\y.r)  (-'(>» 

J 

oi  an  average  over  all  input  ray  directions  weighted  by  the 
relative  Intensity,  Notice,  that  we  have  net  restricted  the 
results  to  which  medium  corresponds  to  air  and  which  to 
water,  When  going  from  air  io  water  set  n - I,  «'  •*>  1.33  and 
when  going  from  water  lo  alt  set  n - 1.33  and  «’  - 1.  Then 
the  computation  of  the  beam  moments  after  traversing  the  sur- 
face yields 

fy')\y.r)dy 

j liy'.rfdy' 
j(y'  -y)*f(y\r)dy‘ 

var  (Vl  (37) 

fK,»r' 

The  results  derived  in  this  section  were  performed  for  e 


one-dime nslonul  surface  To  extend  them  to  a two-dimensional 
surface  is  straightforward  if  we  restrict  ourselves  to  Cartesian 
coordinates.  The  variable  R would  then  become  the  pair  R -■> 
(x,y)  and  the  one-dimensional  results  would  carry  over  lo 
each  of  the  orthogonal  coordinates.  The  interpretation  would 
i hen  be  one  of  projecting  the  true  slope  distribution  onto  the 
Cartesian  coordinate  system,  Although  simple  in  theory  the 
actum  compulations  are  difficult.  If  we  use  the  linearization 
implicit  in  (31),  this  problem  is  greatly  simplified.  For  this 
reason  we  will  restrict  our  analysis  to  this  assumption.  To 
make  the  calculations  for  large  zenith  angles  a more  rigorous 
assessment  of  the  surface  geometry  must  be  performed  [fi|, 

IV,  SATELLITE-TO-SUBMERGED  PLATFORM 

The  computation  of  satelllte-to-submerged  plntform  power 
budget  is  aided  by  a brief  discussion  of  the  geomerty  It  is 
assumed,  for  a varie'y  of  reasons,  that  we  will  project  a spot 
on  the  ocean  approximately  one  mile  In  diameter.  Thus,  if  we 
transmit  Pt  watts  of  radiation,  from  the  zenith  the  full  ungle 
of  the  beam  will  he  approximately  (1/22  000)  rad  * SO 
jirad  = Os  The  power  density,  intensity,  at  the  surface  will 
be  approximately 

p 

_ --=  4.62  X lO-^w/Jlf4.  (38) 

rr(830)a 

If  the  surface  Is  illuminated  at  an  ungle  y from  the  zenith  then 
the  power  density  will  be 

„ 4.62  X 1 0-  'P,  cos  7 w/Af*  - (39) 

rr|830)a 

Now,  however,  the  circular  spot  has  elongated  Into  an  ellipse 
with  minor  axis  H3QA/  and  major  axis  (830/cos  y)M.  We  wifi 
use  the  symmetry  along  the  major  axis  of  the  ellipse  to  |>lck  a 
convenient  coordinate  system.  We  will  cull  tills  the  .t-axis,  The 
minor  axis  will  define  the  y-axls  of  the  coordinate  system  and 
the  depth  of  the  ocean  will  constitute  the  z-a\ls.  Fig,  9. 

In  pructice.  the  spot  will  have  a nonuniform  illumination. 
Wc  will  jiecount  for  this  by  defining  a normalized  Intensity 
/frouful  w,'tch  is  then  multiplied  by  the  factor  in  (38)* 
Notice  that  the  angle  y Is  always  measured  between  the  .t-.v 
plan,  and  t line  in  the  x-:  plane.  This  will  allow  us  to  use  the 
second  form  in  (7).  At  any  location  Uo.T’o)  in  the  x-y  plane 
an  elementary  surface  element  .toy,>  contributes  an  amount 
Hx0,v0)dxodyo’  A ruy  with  this  intensity  passing  through  the 
alr/sea  interface  yields  the  value 

/( Y ' ' l.r  - (v  • • v0 ))  ~ l{x0.y0  Ux0dy0p(  y'ly)  (40) 

as  the  Intensity  on  the  water  side  of  the  boundary.  At  this 
point  we  will  consider  a functional  form  for  p{y  ly)  to  aid  In 
ihe  computation.  From  experimental  result#  1 10]  il  can  be 
assumed  ilut  ply  I y)  is  Gaussian  Futlhemtorc,  If  we  keep  to 
the  angles  such  that  (31 1 Is  valid,  then* 

1 To  be  cotract,  Ike  ttsnmiulon  coefficient  it  th«  boundary  should 
■lu  ba  Included  u a function  of  tngie  (9) . 

1 To  uaa  correlated  Gaunter  variable!  would  only  be  the  refineniant 
of  an  approximation, 
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l-'lg.  0.  Axes  of  coordinate  system. 
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where 

n ' „ J*1  + f >/TO1  1 

w “ : + mi'  + r>  Lv.-  + (.v„  - .v)*  + O'lr-TPj 


pir'iy) 


For  this  approximation 
/(  Y ’ ; <-v  — .v0 — >‘o )) 

Hx0ly0)dx0dya 

2n  I--,  (var  |/<  | ) 


r'~  - t + 

a 


b..  ..  . ■ M- 

II  2 

I --  fvar  |/e] ) 

f :+  mi'  + i^)] 

W J 

(41) 

L 3 J 

<V2  = *[*»  + (.v0  -*)*  + O’o  -.Tla)  in0* 

[I  + 2 ('*  + b/'  + .V'l'l 

: + Mi'  + n J 


(-»!]• 


vur  |rt| 


n * 

] I - - vur  |/<| 
ii 


Inserting  this  into  (7),  and  using  the  ol'f-uxls  correction  to 
the  intensity  we  And  at  u point  below  the  surface  that  the 
contribution  from  the  Intensity  /<jet),.s'„  > in  an  area  dx0Uy0  at  7.v 
the  point  (.v0..vo.O)  to  the  Intensity  ut  (x,.v,2)  Is 

M(x,y,:\ y'x7y')  ^ 

/(.v0,.v0  )dx0dy0  I , „ , 


,«0 * [ r a + (,v„  — .v  )a  + (.t'o  -.ft*  | 1/2 

f* *0 

jflI(ra  + (,v0  v)2  + (v0  - ,t')a  | a,a 


■r-H)*- 


(rr  Ufi'R,’)1 

~yx)2  + (7y‘  “7y)a  + 0m"* 

, J (x  1*0 -X  I (7x‘  -7»> 


exp  - a|ea  + (,v0  — x)a 


m | Vu0  - x)5  + O' 
(y  l.t'o  ~y  I (t/  ~7y) 
\/(*o  - x >a  + O’o  -V)a 


**  “ - * + ai"-‘  ^0^  I ^ + (X°“*),,W 

e,v  = sgn  : f y = sgn  (46) 

In  (4.1)  we  have  Int  educed  the  new  variables  7/  and  7y‘.  The 
former  Is  the  angle  measured  In  the  x-r  plane  while  the  latter 
(4,1)  Is  t he  unglc  perpendicular  to  the  x-J  plane.  This  set  of  variables 
results  from  a rotation  of  coordinates  of  the  variables  Sr  and 
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0#  by  the  transformation 


C‘o 

?A 


t . t'.»_  ' I , 

%/('■()  - i't5  +• Tt(>  VP  >v 

_bvu  ~_XJ , 

<47) 

n/U’o  ~.»lr+  <*o  ~*)v 

and  represents  the  viewing  angle  at  the  receiver. 

The  necessity  lor  this  rotution  arises  from  the  fact  that  6r 
lies  in  the  plane  described  by  the  points  (*o..VoM-*.>)  and  the 
refracted  angle.  Consequently,  it  is  necessary  to  project  the 
angular  contributions  onto  a common  set  of  coordinates 
before  integration.  Since  the  transformation  is  unitary,  the 
variables  tx‘  and  7V'  are  still  normalized  to  one.  Finally,  we 
see  that 


AJ(x,y,z)dx0  dy„, 


It  is  evident  thut  even  with  the  simplifying  assumptions  used, 
the  model  is  complicated,  Therefore,  in  terms  of  an  experi- 
ment, it  Is  important  to  pick  a geometry  such  that  we  can 
make  further  simplifying  assumptions.  For  example,  an  experi- 
ment using  the  sun  at  zenith  would  Ituve  y * 0,  /(JCo,.v0 ) " 
constant.  ,v  « y » 0.  We  could  also  pick  a culm  day  so  that  we 
can  assume  Rx  = Ry  <*  V‘  =*  0.  If  In  addition  we  collect  over 
a sphere  with  u llshcye  lens,  we  have 

/<*„,. Vo,:  I 

• | 

“ fj  rrff  ‘ue*p  ' +-*«a  +VIW* 


WUy*> 

II 

] .V  y - “ > y * 0 

Ri*  3 (2/3 )sOa  Is®  + ,*0a  + .v0al  ',/a 


u'*0  dy0 


Setting. v0  - t'o  cos  py0  = r<>  sin  p and  assuming 


(Notice  that  for  * ■ 0 the  approximation  I*  = 0 docs  not 
hold.) 

The  power  collected  at  depth  z will  be  merely  ,4/(0, 0.z) 
where  A is  the  si/e  of  the  collet  ting  aperture.  Consequently, 
a measurement  of  /t0,0,z)//(xOl.Vol  over  many  extinction 
lengths  would  indicate  the  validity  of  extrapolating  the  model 
lo  great  depths. 

We  have  now  presented  three  separate  methods  for  comput- 
ing the  power  lost  to  a depth  z when  the  source  is  at  the  zenith 
and  no  other  effects  are  considered.  By  order  of  expected 
accuracy  they  are  (4(>).  (50),  and  (7)  when  the  heuin  radius  r0 
is  considered  large,  In  Fig,  10  we  plot  (50)  as  a function  of  the 
upper  limit  of  Integration.  Notice  that  In  all  cuses  convergence 
occurs  when  the  radius  Is  approximately  z/2  for  P o o.oi,  As 
Pancreases  from  0.01  -0. 1 1 the  effective  surface  area  in- 
creases and  the  total  contribution  decreases.  A calculation  of 
(4(>)  was  also  made  and  the  result  was  within  a few  percent  of 
that  calculated  by  (50)  for  P » 0,01.  Finally,  when  we  use 
(7)  with  r0  large,  it  can  be  easily  shown  that /(O,0,z)//(x0l>»0) 
is  merely  e "ai,  This  Is  plotted  together  with  theprevlous 
results  in  Fig.  1 1 as  a function  of  r At  300 M and  P ■=*  0.01 
the  difference  was  only  3 dB  This  result  Implies  that  the 
diffuse  reflection  coefficient  |6]  when  measured  at  the  zenith 
is  approximately  the  absorption  coefficient. 

In  a practical  system,  one  will  encounter  background  noise 
arising  from  the  sky  and  the  sun.  When  this  occurs,  the  use  of 
a 4 it  sierudlan  collector  will  admit  an  unacceptable  amount  of 
noise  lino  the  detector  circuitry.  For  these  cuses  it  cun  be 
shown  thut  to  optimize  the  received  slgnul-lo-noise  ratio  a 
spatially  matched  lllter  should  be  used.  Simply  stated,  the 
matched  filter  will  take  two  forms  depending  upon  whether 
we  have  blue  sky  or  the  sun  (or  both).  To  eliminate  a source 
such  as  the  sun  the  niter  reduces  to  an  obscuration  covering  the 
Held  of  view  subtended  by  the  sun  to  the  receiver.  For  an 
extended  source,  the  (liter  tukes  on  the  angular  distribution 
subtended  by  the  source  to  the  receiver.  In  practice  this 
reduces  to  an  obscuration  which  only  passes  thut  portion  of 
the  field  In  which  the  major  portion  of  the  source  subtends. 
Muthemuticully.  we  would  Integrate  (48)  over  the  variables 
7*'  and  7),'  with  the  Integration  boundary  determined  by  the 
receiver  field  of  view.  Then  (43)  cun  be  rewritten  as 


Mxk\\:\  7v'-7v') 

/(v0.v0W.v0c/i'o 


t IT (V*,  ')» 
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we  have 

»xn  — I tflz®  4-  r«a  I ^ ^a ,+ — -—* 

/(0,0,z)  fm  1 (2/3)jfla[sa  +r0a 
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If  wo  assume  (lull  the  nveivei  will  ho  punning  ill  the  rol'ruoied 
angle  (/;/»(')>  we  can  porloim  ilio  in lop.i ;il n >n  nvet  u I'm itc 
1'iokl  of  view  between  sny  (mw  1)  A anil  I «■'»’)■)  ' A ll  wo 
ivil'iiim  l lie  inlegialioii  nvoi  a oi'tio.  wo  I mil  I li  al  sumo  till- 
iKiiln  wmilil  arise  in  living  in  obtain  .1  ilnsoil  hnin  miIiiIiiiii 
llowevei,  in  leleimig  In  I II!  I.’  wo  soo  lli.ii  nppoi  and  lowei 
hounds  oan  oiisily  he  obtained  in  closed  hum.  I lie  resultant 
leceived  pi.) wo  1 nvoi  the  finite  field  nl  view  li  mu  ho  nhlained 
lining  combinations  nl  the  function 


/t.v 


•'•••• '2|:-  JJ  I"*1’  o|:J 

>0'>«  L 

*- 


+ (.V  V0la  f (.»'  '1 


where 
a n 


C,i:  I l»  ’V  I _ I j 

\fix~~ a:o“)*'  + u*  .»’«»*„  <vl 


fy  l.t'n  I P,n  ' | 


x/f#  -/0T+(.r  - 5o>s  j/v 

For  small  fields  of  view  , O'  01111  ho  loplaooil  hy 


(1 
(52) 


(551 


(54) 


•» 

p (sco  0 ■■■  I ).  see  0 = 1,115 

1 

ami  1 he  maximum  timo  dilToionoo  is 
h‘  l.55r 

M t>  |n.l  I5|  • 5 X 10  (55) 

c 2c 

Ai  : - MMM,  ±T  - l>()  ns.  If  the  primary  contrihui Ions  come 
from  Iwloo  ilie  disk  diameter.  A7" » 524  ns. 

Il'pposiud  has  cnmptiied  the  impulse  response  of  the 
medium  from  which  he  evaluates  the  delay  spread  us  the  l/r 
point  < litis  value  lakes  the  form 


A/  “ 


[i  + :(i„>/l)a, 

(s/I -')i+“(:«  ~sfl) 


(5(5) 


hoi  ; * ,5004/.  this  \ iokls  W ns. 


V.  SlIBSlIRfAl't-TO-SATELLm;  BUDGETS 


The  pari  of  (he  system  must  dtl'tlcult  to  model  has  been  the 
suhsuifacc-to-satclllte  uplink.  This  difficulty  cun  host  be 
understood  by  showing  why  the  two  arc  no  I reciprocal.  On 
the  downlink  a one  mile  spot  projected  from  22  000  miles 
represents  an  antenna  gain  of 


4tr(22  000)* 


00  d B 


On  the  uplink,  howcvci,  if  we  go  through  one  scattering  length 
ol'wutci  the  Hearn  solid  angle  will  be  approximately  01  or  un 
nmcimu  gain  of  4 rr  li- . Since  6*  - 10'  12  111 ' * the  gain  is 

only  21-51  dll.  The  gain  then  goes  down  as  10  log  A',cill 
in  seuiiering  lengihs,  Because  of  the  paucil)  in  gain  It  would 
be  necessary  10  operate  the  system  closer  to  the  surface  on  the 
uplink  than  for  the  downlink.  If  the  scattered  radiation  us 
described  by  the  Heggestud-Artuish  approximation  Is  used 
exclusively. 

Fur  this  portion  of  the  link  ll  is  necessary  lo  investigate 
the  radiation  In  greater  detail.  To  do  so  ll  Is  helpful  to  use  the 
norniuh/ocl  version  of  the  mutual  coherence  function  ( MCI7) 
IIIJ-II5I  (spatial  covariance  function).  FVi  the  scattering 
function  described  In  (51  this  becomes 


7(p)  -1  exp 


"''h'oV 


Vi  i-  <A0(ipff3 


1 

J’ 


(57) 


On  a final  note,  l(  Is  possible  lo  obtain  an  estimate  of  the 
pulse  spreading  by  referring  lo  Pig,  15,  lor  the  /.cult It 
geometry.  If  the  primary  contributions  come  from  the  disc 
with  diameter then  the  maximum  path  dlffccneo  is 


Notice  Ihul  nl  p - 0 iliis  is  normull/ed  to  unity,  and  we 
ussumed  a Gaussian  source  wtilt  an  apeilure  eipial  to  nr„- 
focused  a(  Infinity  In  normal  system  design  we  are  commonly 
Interested  In  the  henmwidth  of  Ihc  antenna  defined  a'  the  5 


'..IMttIRti 

(U) 


ih) 

Fig.  10.  liquation  (30)  plotlid  «» a function  ol'  the  upper  limit  of 
Integration.  (•)  X *•  100  m.  (b)  / - 200  m. 


JB  points,  Since  MCI1  It  the  transform  of  the  nnKulur  distribu- 
tion of  the  source  us  seen  by  the  receiver,  we  cun  ci|ulvulcntly 
do  tine  u "coherence  length"  pr  us  it  ctmipurublr  mcusurc  of 
untennu  ctdllmution.  Thus,  the  greater  the  coherence  length, 
the  closer  the  source  appears  to  approximate  an  Impulse  In 
angle  (point  source),  by  setting  the  MCF  equal  to  e "°'H1 
( 3 dB)  and  solving  for  p,,  we  cun  Investigate  the  behavior  of 
the  radiation  as  it  traverses  the  scattering  medium.  The  expres- 
sion for  pf  becomes 


P,  - 0.693 


0,693 

Z< 

,s 


p,  =■  X 


[|l.3H6(Jil)"0.4H)/fl?)»(g 
ln[tZ  - 0.693] 


0.093 

/>  

J (58) 


which  Is  shown  In  Fig,  14  as  a function  of  l for  an  Initial 
divergence  of  I0“3  rad,  and  for  the  water  ptoperties  defined 
by  i ■ 0,06,  a **  0.04,  and  ?■«  0,01.  Notice  that  for  a distance 
'/.  “ 0.693/»  the  beam  propagatea  at  It  would  In  vacuum,  and 
the  correlation  length  Increaiei  aa  the  beam  diverge*.  However, 
the  scattering  mechanisms  abruptly  taka  hold  at  this  distance 
and  the  cohetence  length  decreaaei  dramatically  In  a very  short 


distance,  and  rapidly  approaches  the  vuluc 


defining  the  Heggcstud-Arnush  approximation,  Since  thii 
behavior  Is  dependent  upon  the  scattering  properties  of  the 
water.  It  is  instructive  to  define  the  albedo  [2|  w us  the  ratio 
j/(j  + a)  und  the  extinction  length  us  N m aZ  =»  (a  + ,v)Z  and 
replot  Fig.  14  for  vurlous  wuter  parameters  in  Fig,  15.  finis, 
we  see  thut  we  rapidly  lose  the  guln  (or  imaging  capability) 
of  the  medium  us  we  truverse  a few  scattering  lengths,  which 
can  vurj  In  terms  of  the  extinction  length,  This,  however,  is 
not  the  whole  story. 

If  we  observe  the  MC'F,  (57),  lor  large  values  of  p,  we  ob- 
serve the  usymptotlc  value  of  < '•*.  From  Fourlet  transform 
theory  we  know  that  this  corresponds  to  u point  sou  me  which 
relates  to  the  unscattered  portion  of  the  beam.  And,  while  the 
power  associated  with  this  portion  of  the  beam  is  significantly 
leu  than  that  usoclated  with  the  scattered  radiation  (lower  by 
*~'2),  It  nevertheless  retains  the  full  gain  of  the  original 
source.  Consequently,  It  cen  be  shown  that  for  the  uplink 
geometry  a receiver  locited  out  of  the  icetterlng  media  at  a 
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Fig,  14,  Coherence  length  11  a function  of  propagation  path. 


great  distance  from  the  source  will  always  collect  more  power 
from  thli  unscaltered  component  that  from  the  scattered 
component  |12).  It  It  therefore  possible  to  comtder  the  up* 
link  radiation  u compoted  of  two  additive  Gaussian  terms. 
The  first,  retaining  ill  the  geometric  properties  of  the  ndlited 
source,  but  attenuated  by  the  factor  *-<•+*)*,  and  the  second 
consisting  of  the  scattered  portion  of  the  radiation  as  con* 
sldered  for  the  downlink.  Either  of  these  terms  may  be  used  to 
develop  a system,  but  the  resulting  systems  will  huve  vastly 
differing  operating  scenarios  due  to  the  difference  In  coverage 
and  pointing  requirements,  Consequently,  we  will  consider  the 
general  problem  of  a Gaussian  beam  propagating  up  through 
the  alr/sea  interface  and  determine  the  effects. 

We  will  first  make  the  computations  outlined  in  (36)  for 
the  Heggestad-Arnush  approximation  and  then  show  how  the 
unscattered  result  follows. 

To  compute  the  surface  Irrudlanee  profile  upon  passing 
through  the  Interface  on  the  uplink  It  is  only  necessary  to  In- 
sert (7)  Into  (36).  This  yields 

A y's)  ■ y‘i  r>/(  y/)  (bo) 

[we  assume  a collimated,  zero  cross  section  Illuminating  source 
(/  y =*  0)1 . Finally,  to  determine  the  angular  distribution 
of  the  beam  (tV)  Is  integrated  over  the  surface  to  yield 

/Tr‘)*  //lr'/)«fr-  dr.  (ei) 


Istlrwtlen  lM|Un  (10*  Water)  N,, 

Pig.  I J.  Coherence  length  u>  a function  of  extinction  lengths  for 
various  values  of  the  albsdo  and  »4. 


t 


To  perform  the  Integration  In  (60),  several  factors  must  bo 
taken  Into  account.  First  recall  (hut  we  hove  defined  the 
medium  containing  the  source  to  have  the  Index  «,  antheenllh 
angle  y.  Thus.  If  we  wont  a more  figurative  description  we 
should  turn  the  coordinate  system  In  Fig.  9 upside  down  to 
yield  Fig,  16. 

Next  recall  thut  we  projected  ihe  true  slope  statistics  of  the 
surface  onto  the  x and  .v  coordinates,  However,  the  scattered 
beam  has  circular  symmetry  with  regard  to  the  angular  diver- 
gence, Consequently.  It  Is  again  necessary  to  rolute  the  axis  of 
the  angular  coordlnetes  by  the  traniformatlon  In  (47).  which 
wtU  allow  ut  to  perform  the  Integration  In  (60)  In  C artesian 
coordinates.  Taking  the  latter  remarks  Into  eccount  uml  again 
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assuming  Gaussian  slope  stuiislics  wolds  the  function  /It /I 
with 
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/ 


*v.-[4  v;u  lA> I J 

exp  1 v5  l ,r2  ‘ 

I A 1 

■[(*'- ;■•*'('■ -?)*■)’ 

(•  ;>)‘] 


vai  |W) 


/l  y .M  “ 


/>/ 


- exp  | ti\r:a  i va  * i 5 

l ■ 


/'<  Y .M' 


1 

7 exp 
'v  1 

• exp  ■ fwv'e5  + .v-  + f-1 


/’i 

rK,’* 


HI 


,v’ 1,11 

0*,/  (*  ')*' 

*.v  l.vio,,,’  Y J i n 

' V^'  + y1  I Ix/v^T  yi  [ J 


y/x 


//  'a 

1 O' 


HO) 


where  we  huve  assumed  thut  the  source  Is  located  at  i.v,i  *»  0. 
,\  0 » 0,  :l  and  the  surl'acv  radiance  profile  Is  over  Con- 

sequently, the  angular  distribution  of  the  emerging  beam  Is 
obtained  by  Integrating  yit'.fl  uvet  the  variables  (.Y,.t'l.  where 
■>,, ' and  are  now  the  angles  projected  b>  the  source. 

As  the  surface  roughness  goes  to  zero  In  (<0|.  /hy.Y* 
approaches  the  delta  function  fiiY-m  //’ >Y>-  * "r  this  case  the 
integration  cun  he  performed  ovet  fiA.  0,  to  yield 
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(O.l) 

We  cun  ulso  perform  the  integration  m HO)  Ionite  trurrow 
beam  eusc.  For  this  euse  wo  can  extend  the  Integration  from 
» to  « vteldlng 


Pot  the  case  where  the  ocean  roughness  l*  absent  to.M  we 
see  that  the  exiting  beam  Is  centered  around  the  linearized 
Snell's  angle.  («/«')> ,v.  hut  somewhat  steered  toward  the 
zenith.  Physically,  this  Is  due  to  the  fact  that  the  scattered 
paths  which  lie  closest  to  the  zenith  traverse  a shorter  distance 
and  consequently  are  uhsoihed  less  which  shews  the  beam, 

We  can  Integrate  the  contribution  at  an  angle  y over  'he 
entire  surface  us  Indicated  in  (hi  I.  This  results  In  the  function 
/Iy’i  which  Is  lire  angular  power  distribution.  To  compute  an 
uplink  budget  one  needs  to  integrate  l\y  l over  the  solid  angle 
subtended  by  lire  collecting  uperture  At  a distance  R large 
enough  for  the  function  y( y I to  be  constant  over  the  collect- 
ing upeilurc.  this  solid  angle  Is  merely  .4  R*  Consequently 
the  collected  power  Is 


l>,  A Y ) 


Ra  ' 


fl'li  | 


t<'4! 


where 


An  unporlunt  case  to  notice  Is  when  the  beam  is  exiling  the 
wain  at  the  zenith  and  the  surface  is  smooth.  I hen  tt'.'l 
miegtuted  over  the  hemisphere  yields  the  result  in  (4'M,  mult i- 
plied  hy  (lie  iiomiallzcd  angular  distribution 


(t'7| 


If  we  look  at  the  .1  dl)  contour,  the  beam  half-angle  is  i n/n  l \ 


( Oh<)T,  and  i lie  effective  gain  ii 


(68) 


The  link  loss  for  lliis  ease  becomes  (?/  * yy'  «•  0)/,,  where 
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To  obtain  (ha  remits  for  the  unfettered  beam  we  would 
use  (6)  In  piece  of  (7).  This,  however,  merely  requires  the 
substitutions 
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with  0m  ■ rm  «•  0.  With  these  substitutions  (62H64),  (67)- 
(69)  become  (71M76),  respectively, 
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We  aun  point  out  that  for  large  zenith  angles,  the  linear* 
Izatton  tied  to  derive  these  results  is  not  valid,  and  t more 
detiileujoiiysis  is  required. 

Condoling  (68)  and  (75)  and  letting  cos  X ■ 1 - (X*/2) 
we  see  tut  the  ratio  of  the  gain  in  the  unscattered  to  scattered 
beam  or  ait  Is 


jV,£ 

(A/i„V 


(77) 


where  w uve  assumed  that  the  beam  Is  focused  at  infinity 
(do  “ .‘ffil  and  we  are  transmitting  from  N,  scattering 
lengths,  'Hocc  the  unscattered  beam  base  N«  times  the  power 
of  the  tottered  beam  we  see  that  the  Inequality 


(78) 


determine  which  portion  dominates. 
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APPENDIX  B 

LAYERED  WATER  CORRECTION 

When  a and  j Me  function*  of  z,  It  bgcomes  necetiary  to  determine  the  effect  upon  the  model  and 
to  cornet  for  It  when  necessary.  Write  s(z)(M(z)  u 

«<a)7*<x)  ■ *1®?  I st«  (*i_l»  zl)*  (B*l) 

I^t  u»  flret  consider  the  effect  upon  the  angle.  For  constant  parameters,  the  angular  vMianoe  grows  as 
*0*2*0.  Since 


us  the  exponential  absorption  loss  to  replace  ax. 

The  spatial  spread,  on  the  other  hand,  goes  us 

s U2  z3 

which  ii  actually  SI  z2.  Now 

(fl  + Aft)(z  + Ax)2  • flz:  + Anz2  + 2flzAz  + AftAz2+0<Az2) 
or 

A (O  z2)  ft  z2  + AO  z2  ♦ 2 fl*  A » + 0(A »2)-.0-?.-  « 3 g(z)  z2  J2(l)  . 
A z A z 

Hence,  the  spatial  spread  can  be  modeled  as 

•2 


f 3s(z)z2t^(z)dz  * ijd|2  (z|3  - z3j-i). 

Jo  1*1 


which  becomes 


‘ N 
l*j 


»|82  (*i3-*3|-l) 


sec3  0 


(B-12) 


to  account  for  the  slant  angles.  Better  approximations  of  the  integrals  could  be  obtained  with  linear  Inter- 
polation, but  was  not  considered  necessary. 
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APPENDIX  C 


SHALLOW  DEPTH  CORRECTION 


In  this  appendix,  we  will  develop  un  empirical  correction  to  the  model  in  Appendix  A which  will 
apply  in  the  region  of  1 to  10  scattering  lengths.  There  will  actually  be  two  corrections  presented;  one  for 
illustrative  purposes  and  the  other  for  data  reduction. 

In  figure  14  of  Appendix  A,  we  tee  that  the  correlation  length  p of  the  model  used  is  always  less 
than  an  empirical  (It  to  the  volume  scattering  function.  Since  the  measure  of  radiant  spread  is  Inversely 
related  to  the  correlation  length,  this  Implies  an  overestimate  of  the  radiant  spread  und  hence,  a conserva- 
tive model.  We  know,  however,  tiiat  the  correlation  length  defined  has  the  form 

* * . 


which  can  be  rewritten  os 


*e  “ PH-A 


TTm\ 


We  can  define  un  effective  value  for  0^  ns 


(Tim 

.V 


and  rewrite  pr  us 


*>e  * P|l-A  " 


; sz>  .693 


F is  plotted  In  figure  C-l. 

This  correction  gives  the  value  for  the  correlation  length.  However,  since  the  radiant  pattorn  used 
tomuins  Gaussian,  it  turns  out  to  be  only  slightly  better  than  tile  uncorrected  function,  hut  falls  fur  short  of 
reproducing  the  comet  correlation  function  und  hence,  correct  radiant  pattern.  A better  approach  would 
be  to  try  u functional  tit  to  the  actual  correlation  function.  In  addition,  since  there  already  was  a Gaussian 
solution,  we  will  try  to  make  tills  functional  fit  by  adding  another  Guusslun  term  In  u judicious  maimer. 

The  function  that  we  are  trying  to  fit  is 





sz 


- 1 


(C-5) 


where  the  exp  (-  az)  multiplier  was  omitted.  It  was  already  pointed  out  that  the  asymptotic  value  of  this 
function  is  exp(-sz)  so  that  we  are  trying  to  lit 
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By  observing  several  values  of  this  function  (figures  C-2  through  C-S  *.  the  scattered  term  in  equation  (09) 
is  always  visible  for  small  values  of  p.  This  term  will  always  have  the  form 


(t)2) 

Now,  if  the  following  is  set 

P * 3 Pq  * Pc0  ■ ~ 


2tt v szd*1 


(C-7) 


(C-8) 


then  tills  por'  ion  of  the  coherence  function  at  most  contributes  1 percent  to  the  total.  Consequently, 
whatever  remains  of  equation  (C-6)  at  this  value  of  p can  be  considered  to  be  one  value  for  a new  function. 
If  this  new  function  takes  the  form 


then  we  ore  clearly  solving  for  the  point  where 


Similarly,  a second  fit  at  (he  value  p - 3 pc0  yields 


Solving  for  A and  B yields 


(C-9) 


(C-IO) 


(C-ll) 
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Figure  C-3 
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By  defining 


^ 5 

°cff  " 3i 


(C-13) 


we  are  able  to  use  our  previous  Gaussian  solution  without  any  changes.  Ignoring  absorption,  the  relative 
weights  for  the  three  terms  become: 


exp  (-sx)  for  the  residual  image, 


A for  the  third  or  glow  field  term  using  dfa 
(l-A-e“M)  for  the  scattered  term. 

Values  for  A and  B are  plotted  in  figures  C-6  and  7.  Notice  that  0^  > and  the  glow  field  makes  its 
greatest  contribution  In  the  region  of  1 to  10  scattering  lengths.  The  fact  that  A starts  to  increase  after  10 
scattering  lengths  is  attributed  to  the  fact  that  after  this  point  the  Gaussian  fit  for  the  third  tenn  is  no 
longer  valid  for  a number  of  reasons.  As  a consequence,  a modified  value  for  A was  used  which  takes  the 
form 


A'  - .235  (s?.)  o' 


(CM  4) 


This  is  also  plotted  in  figure  04.  Finally,  we  have  also  shown  the  values  for  A and  B when  the  Becond  cali- 
brate point  Is  taken  at  p • 5pco.  This  is  also  plotted  in  figures  C-6  and  7 and  shows  no  appreciable  difference 
other  than  the  calibrate  ut  p ■ 3pco. 
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VUtBlLlTY  LA  MORATORY 

ICIUBPI  INSTITUTION  Of  OCSANOOKAf  KY 


TECHNICAL  MEMORANDUM 


SAN  MXOO,  CAUfOBMU  lit  IS 

23  dune  1976 
ML-76-005t  - REVISED 


SUBJECT:  Ocean  Optical  Proparties  — 520  Nanometers 

Santa  Catalina  Island,  Lat:  33°27.2'N,  Long:  118'29.0'W 


1.  Near  surface  daytime  optical  measurements  of  the  water  properties  at 
the  .Santa  Catalina  Island  SATCOM  test  site  were  limited  to  the  beam  trans- 
mittance and  Its  derivative,  the  volume  attenuation  coefficient,  a.  This 
limitation  was  Incurred  because  of  the  effects  of  ambient  daylight  on  the 
Instruments  used  for  measuring  scattering.  Measurements  of  all  required 
properties  were  obtained  at  night,  however,  and  correlations  between  them 
were  explored  for  a variety  of  water  conditions.  Curve-fitting  techniques 
were  applied  to  these  night  data,  and  empirical  relationships  were  derived 
allowing  the  prediction  of  the  required  but  unmeasured  daytime  scattering 
and  absorption  properties  of  the  water  from  the  measured  beam  transmittance. 
The  fits  of  the  nighttime  data  to  the  derived  analytic  expressions  was 
generally  good-to-excellent,  with  correlation  coefficients  of  0.946  or 
better.  These  expressions  may,  therefore,  be  used  with  confidence  to  pre- 
dict the  unmeasured  properties. 

2.  NOTE:  A window  correction  factor  of  0.9794  has  been  applied  to  tfll 
transmittance  values  obtained  In  the  field.  Any  preliminary  data  Issued 
prior  to  12  April  1976  were  not  so  corrected.  Any  original  transmittance 
or  alpha  profiles,  l.e.,  those  obtained  In  the  field  with  the  x-y  plotter 
or  data  logger,  must  have  this  correction  applied.  All  computer-generated 
output  has  now  been  corrected. 

3.  Using  the  expressions  obtained  from  the  nighttime  June  and  July  data, 
daylight  transmittance  profiles  have  been  used  to  obtain  estimates  of  the 
scattering  and  absorption  properties  of  the  water  column  to  a depth  of 

50  meters.  These  values  are  listed  for  the  various  depths,  Z(m),  at  which 
the  transmittance  data  were  sampled.  The  following  Is  a description  of 
the  information  supplied  In  the  computer  listing  titled,  "Ocean  Optical 
Properties:" 

a.  The  transmittance  on  the  computer  listing  under  the  column  T(l/M)  Is 
the  field  data  multiplied  by  the  correction  factor  (see  paragraph  2 above). 


• !•  < 
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b.  The  column  headed  ALPHA*  is  obtained  by  taking  the  natural  logarithm  of 
the  reciprocal  of  the  corrected  transmittance,  l.e., 

ALPHA*  * ln(l/T) 

c.  The  true  volume  attenuation  coefficient  o (column  headed  ALPHA)  Is 
greater  than  the  apparent  volume  attenuation  coefficient  ALPHA*  by  a small 
amount  which  accounts  for  the  scattering  Included  within  the  acceptance 
angle  (1.5  mill Iradlans)  of  the  transmlssometer. 

d.  The  scattering  coefficient,  a,  (column  headed  SCAT)  Is  calculated  from 
a linear  relationship  of  the  form 

MB  1 *0 

where  a and  a were  slightly  adjusted  from  the  constants  derived  from  the 
straight  line,  least  squares  curve  fit.  Specifically,  from  the  regression 
of  a vs  o (disregarding  two  data  points  with  obvious  difficulty) 


• « 0. 9 58a-0 . 101 , 


0.0315 


and  after  adjustment 


■■  0.97  a - 0.0977,  o§  ■ .0324 

The  reason  for  the  adjustment  was  to  provide  a better  agreement  between  the 
absorption  coefficient  calculated  from  the  relationship,  a ■ a - and  other 
estimates  of  absorption  arrived  at  by  Independent  methods.  As  Is  Indicated 
by  the  standard  deviations  between  the  data  points  and  the  values  predicted 
by  the  two  relationships  (l.e.,  os),  there  Is  only  a slight  decrease  In  the 
precision  with  whlcn  the  adjusted  equation  predicts  the  observations  of  a 
from  a. 

e.  The  absorption  coefficient,  a,  listed  In  the  column  ABS,  Is  obtained  by 
subtracting  the  values  In  the  SCAT  column  from  the  ALPHA  column  since 

a ■ a-a. 

f.  The  volume  scattering  function  a(0)  at  3,  6 and  12  mllllradlens  was 
calculated  using  expressions  obtained  from  the  regression  of  a (0)  vs  a. 

g.  The  normalized  second  moment  of  the  volume  scattering  function?,  l.e., 


2 TT  2 

0 s [ a to)  . o . slnO  do, 

• n 


Is  listed  In  the  column  headed  THTA*  2 BAR.  The  values  listed  were  calculated 
using  the  expression  obtained  from  the  regression  vs  h. 
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4.  The  narrow  angle  volume.  scattering  function  data,  vtzj  VSF3,  VSF6,  and 
VSF12,  wire  obtatned  w.tth  the  ALSCAT  instrument  which  performed  these  measure* 
mants  concurrently  with  determination  of  the  beam  transmittance,  These  data, 
therefore,  are  all  obtained  from  the  same  sample  of  water.  The.  volume  scat* 
torlng  function  data  at  angl6s  from  10°  to  170°  were  obtained  with  a second 
Instrument,  the  General  Angle  Scatter  Meter  (GASM),  which  was  lowered  to  the 
same  indicated  depth  as  ALSCAT  but  was  about  4 meters  away  horizontally. 

Some  vertical  separation  may  also  have  occurred  due  to  errors  in  the  indicated 
depth.  This  mey  have  been  es  much  as  2 meters  on  some  occasions.  Operational 
considerations  did  not  always  permit  data  from  the  two  Instruments  to  be  ob- 
tained at  the  same  time  although  that  was  the  intention.  In  general,  the 
operational  procedure  In  July  was  Improved  over  that  In  June,  and  the  time 
difference  between  the  two  measurements,  l.e.,  ALSCAT  end  GASM  was,  with 
three  exceptions,  reduced  to  less  than  8 minutes.  These  time  and  spatial 
separations  between  the  ALSCAT  & GASM  data  were  usually  of  little  consequence 
with  the  possible  exception  of  those  circumstances  where  the  measurements 
were  obtained  In  the  turbid  scattering  layer  that  occurred  frequently  at 
depths  of  from  IS  to  40  meters.  As  this  layer  was  sometimes  very  localized 
vertically  and  would  move  up  and  down  with  the  passage  of  the  coastal  current 
or  with  the  propagation  of  Internal  waves,  temporal  and  spatial  coincidence 
of  the  measurements  assumed  greater  significance  for  these  cases.  Despite 
these  concerns  there  appears  little  difference,  as  determined  by  the  corre- 
lation coefficient  r,  with  the  tightness  of  fit  of  the  observations  to  the 
resulting  analytic  expressions  obtained  for  the  18  nighttime  data  sets  In 
June,  the  30  sets  In  July  and  the  48  combined  sets.  The  expressions  derived 
from  the  regressions  are  listed  below. 

Expressions  Coupling  Measured  Ocean  Optical  Properties 
(Derived  from  48  sets  of  nighttime  data) 


a - 1.0188a' 1,000  r - 1.000 

s - 0.970  a-  .0977 

a ■ a - 3 


o(3  )-  1398.4s  '• 110,4  r » 0.9597 

0(6)-  353,2s  1 • 1 78 * r - 0.9800 

o(12)-  218.8s  1,2,01  r « 0.9872 

ST  - 0.0283s  r - 0.9467 

5.  The  scattering  properties,  as  determined  from  the  combined  ALSCAT  and 
GASM  data  for  the  48  nighttime  measurements,  are  submitted  herewith.  Each 
sat  consists  of  a plot  of  log-volume  scattering  function  vs.  log-angle  and 
two  pages  of  computed  properties.  The  "Iterated  Data"  differs  from  the 
"Data  Read-In"  by  a correction  applied  to  "Sigma"  valuer,  for  the  three 
smallest  angles,  viz,  0(0)  at  3,6  and  12  mil  11  radians  to  account  for 
multiple  scattering  In  the  one-meter  measurement  path  length  used  In  ALSCAT. 
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The  correction  to  '"ALPHA"  is  that  required  to  account  for  the  inclusion  of 
scattering  In  tne  acceptance  angle  of  the  transmissonmter.  The  column  headed 
"Integral'7  is  the  portion  of  the  total  scattering  coefficient  due  to  angles 
from  zero  to  the  listed  angle.  Thus 

, r° 

"Integral"  (0)  s 2*  I o(0‘)  slno'  do'  . 

Similarly,  the  "Normalized  Integral"  is  the  decimal  fraction  obtained  by 
dividing  the  value  listed  in  the  "Integral"  column  by  the  total  scattering 
coefficient  e.  q 

"Norm.  Integral"  (0)  s o(o')  slno'  do*  , 

Jo 

6.  The  values  of  a,  •,  and  a obtained  from  the  nighttime  measurements 
described  in  paragraph  5 above  have  bean  entered  on  the  curves  of  these 
properties  generated  from  the  regression  expressions.  The  times  of  the 
measurements  ere  printed  beside  the  data  points  to  aid  In  evaluating  this 
comparison  betwean  the  measurements  and  the  values  predicted  from  a trans- 
mittance profile  obtained  on  the  same  night  but  at  a somewhat  different  time. 
It:  will  be  noted  that  the  predicted  values  agree  well  with  the  measured 
values  considering  the  changes  that  are  found  In  repetitive  measurements 
performed  within  relatively  short  periods  of  time  In  a dynamic  ocean. 

7.  Secchl  disc  observations  were  obtained  on  six  days  In  June.  The  dis- 
appearance depths  varied  from  14  to  22  meters.  The  observations  are  listed 
In  the  attached  table.  No  analysl^of  these  data  will  be  attempted."  No 
single-value  parameter,  e.g.,  « or  Tat  a single  depth,  would  be  expected 
to  correlate  as  well  with  these  observations  as  would  an  average  or  "effec- 
tive" value  obtained  by  some  weighted  Integral  of  the  parameter  ovor  the 
depth  Interval. 


•SECCHI  DISC  OBSERVATIONS 


mi 

IIME 

PEPTH 

6/18/75 

1230 

18  m 

1330 

21  m 

6/19/75 

1030 

15  m 

1130 

15  m 

1230 

14  m 

1330 

15  m 

1430 

18  m 

6/20/75 

0950 

19  m 

1054 

19  m 

1205 

15  m 

1250 

15  m 

1345 

15  m 
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TIME 
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6/24/75 

1125 

22  nr 

6/25/75 

1008 

17  m 

1014 

17  m 

1103 

17  m 

1125 

17  m 

1204 

17  m 

1255 

15  m 

1303 

15  n 

1439 

18  m 

1448 

18  m 

1554 

17  m 
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15  m 
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8.  The  errors  in  the  optical  properties  presented  have  been  the  subject 
of  considerable  stud/.  Unfortunately  no  single,  simple  statement  can  be 
offered  regarding  the  magnitude  of  the  errors  as  they  depend  upon  the  pro* 

forty  under  consideration  and  the  value  or  magnitude  of  that  property, 
urtnermore,  the  optical  nature  of  ocean  water  is  highly  variable  both 
spatially  and  temporally  — particularly  so  in  coastal  waters.  Any  values 
of  ocean  optical  properties  presented  as  a result  of  measurements  in  this 
dynamic  systan  must  be  understood  to  iiave  inherent  "errors"  due  to  sampling 
in  space  and  time.  As  mentioned  in  Paragraph  4 above,  the  sanpling  pro- 
blem was  examined  with  regard  to  the  nighttime  data  sets  and  was  not  found 
to  affect^  significantly,  the  observed  correlations  between  the  several  pro- 
perties. The  prediction  of  the  set  of  optical  properties  pertinent  to  the 
SATOGM  tests  by  the  application  of  the  regression-derived  relationships  to 
the  daytime  transmittance  profiles  circumvents  any  concern  with  simultaneity 
of  data  acquisition.  The  remaining  question  is  "did  the  transmittance  pro- 
file change  significantly  from  that  which  existed  at  sane  critical  time  of 
concern?"  This  question  must  be  examined  for  each  situation.  Frequent  pro- 
files ware  taken  in  order  to  provide  a measure  of  the  rates  of  change  ox 
the  optical  properties  and  with  the  intent  of  bracketing  the  SATCCM  observa- 
tions. 


Absolute  errors  in  the  measurement  of  beam  transmittance  were  less 
than  *0.4t.  The  resulting  relative  errors  in  the  volume  attenuation  coef- 
ficient would  be  between  ±2%  for  the  clearest  water  and  tit  for  the  most 
turbid  water  encountered  in  the  upper  50  meters  at  the  site  in  June  and 
July  1975. 

Absolute  errors  in  the  determination  of  the  volume  scattering  'function, 
o(0),  over  the  angular  range  of  measurment  are  somewhat  greater,  probably 
ranging  up  to  101  in  seme  circumstances.  More  significant  to  the  calculation 
of  the  scattering  coefficient  is  the  fact  that  the  volume  scattering  function 
was  not  measured  between  0.7*  and  10°,  necessitating  interpolation  of  the 
function  between  these  values.  As  up  to  501  of  the  scattering  coefficient 
may  be  due  to  the  scattering  in  this  angular  range,  the  magnitude  of  the 
coefficient  is  significantly  affected  by  the  nature  of  the  interpolation. 
Furthermore;  since  the  scattering  coefficient  accounted  for  something  be- 
tween 511  (for  the  clearest  waters  encountered)  and  891  (for  turbid  water) 
of  the  volume  attenuation  coefficient,  the  magnitude  of  the  absorption 
coefficient  --  as  determined  by  the  differences  between  a and  o becomes 
particularly  sensitive  to  the  method  of  interpolation  of  the  volume  scat- 
tering function.  For  example,  when  a ■ 0.6m  , s/a  * .81  and  a/a  * .19. 

Under  these  conditions  a 10%  change  in  s results  in  a 431  change  in  a. 

If,  us  In  the  application  of  these  water  properties  data  to  the  &ATCCM  sys- 
tem analysis,  the  absorption  data  is  particularly  significant,  we  need  a 
better  capability  for  determining  absorption  either  by  more  complete 
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and  accurate  determination  of  the  volume  scattering  function  or,  more 
desirably,  by  a direct  measurement  of  the  absorption  coefficient  itself. 

t 

The  determination  of  the  normalized  second  moment  of  the  volume 
scattering  function,  07,  is  inherently  less  sensitive  to  the  value  of  the 
VSF  at  small  angles  than  is  the  scattering  coefficient.  07  is,  in  fact, 
primarily  dependent  upon  the  VSF  in  the  angular  range  measured  by  the 
General  Angle  Scattering  Meter,  i.e.,  for  G > 10°. 
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Figure  D-l 6.  Ocean  optical  properties  (sheet  2 of  2). 
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Figure  D-23.  Ocean  optical  properties  (sheet  1 of  2 ). 

D-51 


SANTA  CATALINA  IS.  LAT  33-27. 2 N LONG  118-29.0  W 


H'JNTCIVM'l 


ijKIkj™ 

■rmt-M 

■■79  > im 

0.226 

0.108 

266.2 

96.3 

33.1 

0.076  < 1 

0.710 

0.362 

0.366 

0.238 

0.108 

280.2 

102.3 

37.3 

0.072  1 1 

0.66 

0.689  0.373 


tw 


NI'J 

nr.Tii 


0.110  362.9  1)6.6 

0.378  0.269  0.109  320.6  118.0 


I'nui 

M-H 


I'lUU 

■cut 


L'rvMJ 

hi 


mm 

till 


IUI 

III 


0.062 

0.066 


h LTTMi'n'^VJ 


11.0 

0.316 

0.663 

0.676  l 

12.1 

0.333 

0.626 

0.636  l 

■ 0.383 

0.668 

0.113  ! 

■ 0.383 

0.668 

0.113  ! 

19.9 

0.399 

0.312 

0.319 

0.606 

0.113  ! 

23.0 

0.613 

0.686 

0.693 

0.380 

0.112  < 

9.1  0.588  0.332  0.339  0.623  0.116  336.6 


0.036 

0.069 


191.7 

177.6 

72.1 

66.6 

0.031 

0.033 

unm 

■m 

0.030 

171.6 

66.2 

0.036 

t I- W jSHMfir  n 

Figure  D-23.  Ocuan  optical  propertiei  (iheet  2 or  2). 


D-52 


OCEAN  OPTICAL  PROPERTIES  - 520  NH 

WITA  CML2M  a.  UT<  SW7.I  n;  UM>  Ul-tf.o  N 
2ULU73  mm 


Figure  D-26.  Ocean  optical  propartiei  (iheet  1 of  2 ). 


D-57 


OCEAN  OPTICAL  PROPERTIES  - 520  NH 

SANTA  CATALINA  IS.  LAT*  39-27.2  N!  UM>  W-B.O  N 
2U0UN75  1331WT 


i 17.4 

0.661 

0.414 

0.418 

0.2  {4 

0.164 

336. 

.9 

113.8  4 

1 18.7 

0.661 

0.414 

0.416 

0.254 

0.164 

336. 

.9 

M3. 8 4 

K 21.7 

0.664 

0.409 

0.414  0.251 

0.163  • 

1 23.2 

0.663 

0.411 

0.415  0.252 

0.164  2 

i 

35.0 

0.704 

0.351  ( 

1 

36.6 

0.713 

0.33B  1 

I't'9  il'lilil  MM V d:(«I  #!M 


SANTA  CATAUNA  ISi  MT  »-2n*  * ' LONO.  • '116*29.0 . H 


■ iii  i—  in  i — n 1 r i 


W • J f.  -> 

0.552 

0.371 

V • * IV 

0.181 

494,0 

0.477 

0.304 

0.173 

404.2 

■I  i 1 MUIJ 

SSIBuBB 


'2 

0.398 

0.402 

0.241 

0.161 

318.8 

107.0  4 

>6 

0,406 

0.411 

0.248 

0.163 

328.6 

110.7  4 

6 0.075 

2 0.074 


QHES  QHS  \ SQC  BBCES 


XLLB 

■■rrira 

i 354 

0.203 

0.152 

267.6  ( 

.342 

0.193 

0.149 

254.5  1 

i»rr4~ii— «gT4'iLi 


musK  2XjS  SUU 

IK  ! IS 

EoElaDIBBuQO 


51.1  0.765 

«AI)SS  READ  ‘ 


0.266 

PLOTTER 


3 0.100 

6 0.103 


• i u.iu* 

21. S 0.106 


OCEAN  OPTICAL  PROPERTIES  - 520  NH 

SANTA  CATALINA  8.  UTS  SMT.2  1/  UMB<  liKf  .0  Si 

nuum  usmT 


!] 


1 


Figure  D*30.  Ocean  optical  properties  (sheet  1 of  2). 
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Figure  D>36.  Ocean  optical  properties  (iiheet  2 of  2). 
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Figure  D-38.  Ocean  optical  properties  (nheet  2 of  2). 
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Figure  P-39.  Ocean  optical  properties  (iheei  2 of 
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Figure  D-42.  Oceun  optical  properties  (sheet  2 of  2 ). 
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Figure  D-51.  Ocean  optical  propertiei  (iheet  2 of  2). 
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Figure  D-59.  Ocean  optical  propi  rtiei  (sheet  l of  2 ), 
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Figure  D-61,  Ocoun  optlcel  properties  (sheet  1 of  2 ). 
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Figure  D-64.  Ocean  optical  propertlei  (sheet  1 of  2). 
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Figure  D-64,  Ocean  optical  propertlei  (sheet  2 of  2 ). 
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Figure  D-68  Oeem  optical  propertiei  (sheet  2 of  2). 
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Figure  D-69.  Ocean  optical  properties  (sheet  1 of  2). 
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Figure  D-71 , Volume  Kettering  function  (sheet  1 of  3). 
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Figure  D-71.  Volume  icattering  function  (iheet  2 of  3). 
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Figure  D-71 . Volume  scattering  (Unction  (sheet  3 of  3). 
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Figure  D-72.  Volume  Kettering  function  (iheet  2 of  3). 
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Figure  D-72.  Volume  icattering  function  (iheet  3 of  3). 
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Figure  D-73,  Volume  ecitterlng  function  (iheet  1 of  3). 
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Figure  0*73.  Volume  Mattering  function  (sheet  2 of  3). 
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Figure  D-73,  Volume  scattering  function  (iheet  3 of  3), 
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Figure  D-74,  Volume  scattering  function  (sheet  2 of  3). 
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Figure  D*74.  Volume  Kit  taring  function  (iheet  3 of  3). 
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Figure  D*75.  Volume  Kettering  function  (sheet  3 of  3). 


Figure  D-76,  Volume  nattering  function  (theet  1 of  3). 
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Figure  D-76,  Volume  scattering  function  (iheet  2 of  .1). 
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Figure  D-76,  Volume  Mattering  (Unction  (iheet  3 of  3). 
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Figure  0-77.  Volume  Mattering  function  (iheet  2 of  3). 
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Figure  D-77.  Volunu  Wittering  function  (iheot  3 of  3), 
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Figure  D-78.  Volume  loettarlng  function  (iheet  1 of  3). 
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Figure  D-78,  Volume  scittering  function  (sheet  3 of  3). 
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Figure  D-79.  Volume  tottering  function  (theet  3 of  3). 
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Figure  D-81 . Volume  Mattering  function  (iheet  2 of  3). 
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Figure  D-81 . Volume  weltering  function  (rheel  3 of  3). 
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Figure  t>82.  Volume  luttering  function  (iheet  1 oi'.'l). 
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Figure  D-83.  Volume  Mattering  function  (iheet  3 of  3). 
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Figure  D-84,  Volume  icetf^ng  function  (iheet  1 of  3). 
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Figure  D<84.  Volume  Kettering  function  (iheet  3 of  3). 
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Figure  D-85.  Volume  Katterlng  function  (sheet  3 of  3). 
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Figure  D-86,  Volume  tctHtrlng  function  (iheet  2 of  3). 
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Figure  D-86.  Volume  Kettering  function  (iheet  3 of  3), 
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Figure  D-87.  Volumo  icatterlng  function  (sheet  3 of  3). 
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Figure  D-88.  Volume  scattering  function  (iiheet  2 of  3). 
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Figure  D-S8.  Volume  nattering  function  (ihect  3 of  3). 
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Figure  D-89.  Volume  Kettering  function  (iheet  2 of  3). 
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Figure  D-89,  Volume  Kettering  function  (iheet  3 of  3). 
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Figure  D-90.  Volume  Kitterlng  function  (iheet  2 of  3). 
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Figure  D-90  Volume  scattering  function  (iheet  3 of  3). 


D-207 


MwaagnamgaHa 


Sir.MAt  0.0  DEGREES  ) « 279.7 


SLOPS!  3 MILLIARD)* 

S UP  TO  O.L  DEGREES* 

-1.412 

2.A084B-DI  NORMAL IZBD*  1.79697S-02 

-MX!  MUM  PART  1 CL 
| EXPECTED  K/ALPH 

e 01 AMBTER  ( 
'A.  0.929 

MICRONS)*  106.0 

6 EXPECTED  01  PAUSE  ATTENUATION  COBPPICIBNT  - K . 

920  17JUITS  2136  STA . CAT* IS  - SC*  MID  3 OH 


V-UI'IS-W.U'H 
kmiC'ii 


u»i<n  ii'H'ii 


n.i.n 

k.M  .m  T3«i  Hniw.]  n>u 


2.19721*03  1,23896*01  1.I027E  02  3.39831*03  2.84336-02 


fiMriCin 


irvixiiCMi 


inrl.nwrl 


Lll.HM-D'M 


r 

3.48241-03 

1.99336*01 

1.04141 

01 

6.33038-03 

3.0S36C-02 

31 

i • 

4.38411*03 

2.51196-01 

7.3(378 

01 

8.48716-03 

6. 60671-02 

41 

fillNO  1 


knn*iatpi 


r»iif 


6,94838-03 

3.96116-01 

3.92736 

01 

1.32938-02 

1.03488-01 

61 

8.74748-03 

3.01196-01 

2.83718 

01 

1.62)08-02 

1.26348-01 

71 

mmrm  er.> 


,186 
1.743)8*02 
-o 


rT9  Mil  -0*1  !]-■>)• 

in  n r ii  ri— is  i jin 


1.00008  00  1.10398  01 


imuii'i'i 
[•■mm  ■>!•■ 


../HUJ'H 


1.99338  00  4.33318  00  4.67348-02  3.63938-01 


ifuinDu 

rrcRi 


44 

2.79368-02 


mmiQfi 


imuinii 

mviiinrii 


623 
2.17471-01 


ln>nua<ii 


rniiiiK'in 

mv*T<i 


CTT'IU  ■■I'l 


rmiiKiii 

Tfiinmi 


6.94836*02 


■urn 


r.to 

1.38648-01 


IfltWCU 


,61806-0 

3.49078-01 


.16236  00  1 

3.98116  00  1.06311  00  7.04846-02 


mii  j cm  i 


6.30966  00  3.39328-01  8.41878*02 

7.94336  00  2.08036-01  9.02138-02 


KMOTKll 


m-irmi 


roirrai 


2.00008  01  2.33288-02  1.09628-01 


lt mrm'ii 
(imrciii 


3.48678-01 


(•LTlIOdl 


6.33348-01 

7.02246-01 


«fT.Vr*!  0*1 1 


.11868-01 

8.33128-01 


1.04726  00 
1.13458  00 


muiwi 


i • 

1.39636  00 
48 


VUWB1 

Lfcl4-mT4 


1.63816  00 
6336  or 


n vm  a-r-J 


1.91996  00 
2 


■naaii*!?] 

m.tmwT.i 


6.00006  01 
6.30008  01 


■fflU'IttdWII 


.23726*03 

9.91478*04 


kim'EXM 


8.00006  01  3.82432-04 


mWlWJI 

niir.T.ui.ii 


SUVJ-a1!! 

rvn.iJ3iU 


1 . 23288-0 
1.23828-01 


.24661-0 

1.23008-01 


uriiaoi 

mnnu 


9,396 
9.63868-01 


irnw-MH 
IrM 


.70401-01 

9.730ie*01 


inn  D'ii 

uniflni 


9.30006  01  4.40198-04  1.23816-01  9.7933G-01 


IWl'W'n 

MM.1-K71 


LTT+M.l-Q'iJ 

Monan 


1.10008  02  4,09268-04 


HIU'L'll'll 


ltwimu'X.3 

miiicB 


Hvn-u'i 

ini;  m 


oi 

1.26488-01 


imiwii 
ihubii 


wnno'ii 

■TTIWM 


nri3.ii 


2.18176  00  1.23006  02  4.21678-04  1.27068-01 


9.84396-01 


imrWC'M 
■mini 


16 

9.89236-01 


PAUSE  READY  PLOTTER 


Figure  D-9 1 , Volume  scattering  function  (sheet  3 of  3). 
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Figure  093.  Volume  wittering  function (ilieet  3 of  3). 
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Figure  D-96,  Volume  Ktttering  function  (iheet  3 of  3). 
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4>7,  Volume  Kettering  function  (iheet  1 of  3). 

D-226 


riijnn.i.  . "yir  r,».  a ,V  r."  ■ mn~  iir,  r 


I 

I 

I 

I 

I 

I 

I 


7T- — lib  rWOCT?  SWI  JTa.CATWT  - UA  Hid  a oh 

*3'  * • "P  i ■ . • • 

’ ' PATA  READ  IN  ' ''  .■Vl''.. 

• • v ' “ 

• ' • v>.  *•  '■  '• 1 ,3  • ; ■ 

-*  ITERATED  OAT A 

ANGLE  (DEG) 

SIGMa  1NSTA-0 

ANGLE  (OEG I 

SIGMA 

1 

0.17S0 

2.25006  02 

0 

0.  1730 

2. 12706  02 

2 

0.3500 

7.0000E  01 

0 

0.3500 

7.8313E  01 

1 

S'SiBW  /»T*T*  WHP&y v" 

«-A«OOB  01 

n 

n.7noo 

J.MMi  ni 

4 

10.00 

1.5349E-01 

0 

10.00 

1.53496-01 

5 

15.00 

6.6496E-02 

0 

15.00 

6.6496E-02 

4 

»o.on 

3.fl24AB-05  . 

n 

5n.no 

3.024AF-n5 

7 

40.00 

5.03296-03 

0 

40.00 

5.03296-03 

A 

50.00 

2 .60266-03 

0 

50.00 

2. 60266-03 

.9 

.60.00 

1.4922E-03 

. 0 

60.00 

1.4922E-03 

10 

70.00 

9.2919E-04 

0 

70.00 

9.29196-04 

11 

no. oo 

6.61446-04 

u 

60.00 

6.6144E-04 

1?. 

on. on 

3 .335/4B-04 

rt 

4fi  -fin 

4.3356B-fl4 

13 

100.0 

4.64952-04 

0 

100,0 

4.66956-04 

14 

110.0 

4. 54646-04 

0 

110,0 

4.3464P-04 

1 3 

150.0 

4.7  AA 1 A .04 

0 ... 

150.0 

4.7AA1B-04 

16 

130.0 

5 . 1 139E-04 

0 

130.0 

5.H39E-04 

17 

140.0 

5 ,39646-04 

0 

140.0 

5.39646-04 

1 A 

130.  Cl 

A.5315B-04 

...Q  ...  . 

1311.0 

6.5315P-04 

19 

160.0 

7.6257E-04 

0 

160.0 

7.62576-04 

20 

170.0 

1,10426-03 

0 

170.0 

1.10426-03 

21 

1 

1B0.0. 

1.2323P-Q3 

ALPHA-  0.2792 

S/ALPHA- 

0.628 

* 

S«  0.1733. 

A/ALPHA* 

0.372 

A-  0.1039 

e/s- 

0.019 

CORRECTED  ALPHA 

CQRR6CT10N-0.003 

?P  JUN  1976  1141.51 


L. 


TTSDCmwr  5U.dAT*il  " “sTT ,«2o  40*  , " RH  integral 

Hv.irii — PS8M=& — ilffTrir-w — r:mtrw — feiSSi=S — rr 

ISS*  - ga  a SI  SH  i\ 

V.iKilffi  i.Iu»i=M — i'.Miinr 


I.A4UE-S7" 

i.*i93e-oi 


nTTETO^r 

1,03786-01 


TT 

Si 


, Jfmj)  3.16236-01  9.06476  01  l. .1936-0*  rr 

81  — H 

r.ifnaB-o?  !*2e5?B  «}  lugoill-oi?  2 1 2^36-01 n 


1.35646-02 

T7WTW 


?.  ,19726-02 
2. 76626-02 
‘■j:4nnT-7tt 


7.94336-01 

T'.flWOOE  00 


1,2589?  00 
1.58496  JO 
' rTTOTT 


2.34736  01, 

1.4^176  (51 


1.1739E  01 
8.00546  00 
;.34g5'?"~^r 


3. 9064E-02  2.2295E-01  41 
T.4K7R-02  i.i524E-fii  t»T 
5. 39796-02  3.0792E-01  111 


6.25176-02 

7.14fclb-W 


3,56636-01. 


121 

TTT 


2.51196  00  3.52206  00  8.12716-02  4.6361E-01  141 

j^soljll  I-Soooe  fll  6,64946-02.  i i 4^SI=2L 

z.oouob  oi 


y.a?iE'6-gy  

4^36336-01  2.50006  01  t. 70786-02  3*446-01 

5 . 2360E-01  3.00006  01  1.0793^-02 XaAES|£tSJ 

-r1*  V fifths -ol ftsfiooFTt tTT^TIT-o'  iTn?.  eet-ol 

6*94136-01  4,00006  01  5.03296-03  1.646B6-01 


6.9413E-01  4,00006  01  S.W.'-t-uj  , o Ib»6B.01 

mm-  mm-mi  isfe  ft&hr 


9.59936-01  6.50006  01  1.746*6-03 


r 

0. 60976-01 

»iWTtiQ  1. 


246 

JLzJL 


2.18176  00 
2.26B96  0U 


1 ,25006  02 

1 ,30006  02 


2.26R9E _00 J . vc, 

- ihi’hra  1.35661  5T 

2.44356  00  1.40006  02 


4.95126-04 
5.1 1396-04 . 


1.73646-01 
1 .73865-01 

i.nofii-br 

1.74266-01 


1hKie=gr 

2.44336  00  1.40006  02  5.39646-04 
~ no  1.55006  02  6*fl2206—0^  1«7479b  Dl 

S - Jam  MttKtt  l:1«W=« 
*•)««; 

1 ,S‘i/!EI:-b3  " l.V536fc-ffF 


9.9175F-01 

*?f#53nF=crr 


it 1 


T9irw 

2.96716  00 
3.0943E  00 


1.75006  02 

I'.urlrw  i ,60006  fji' 


T5u 


PAUSE  READY  PLOTTER 


*;  r 


_ l t-i/TJ  fiA,.- 

1.13436 

00 

6*50006  01 
f . n n n n P n 1 

1.16666-03 

9.89196-04 

1. 69896-01 

1.7 0916.10  \ 

9,69136-01 

.,  .WglffiAL- 

256 

m 

- 

!"  ' ' 1.30906 

\ .397,36 

JjJ2 — 
00 
00 

IlsSobe  oi 

p.uoooe  01 

H. 50006  01 

7.76346-04 
6.81446-04 
6.  10666-07.  _ 

^.TOB66-0“ 
1,7 1256-0 1 

.444^41 

9,74656-01 

9.76P76-01 

266 

271 

lii  - • 

- 

TTTTMF 
l .65816 

r-a_ 

00 

9, 0000 B Ol 
9. 50006  01 

4 « 5ii6fc-ct4 
5.05016-04 
<,.64956-04_. 

1,71926-01 
1, 72216-01 

9. 80o9E-0 l 
9.R233E-01 
9. 8 3646—0 1 .... 

266 
, Sil_ 

1 . 7 ^ 5.16 

l . 83266 
1 .9199E 

■ l .QQ ilit 

J2I1  — 
oo 
00 
M 

1.05006  02 

1,10006  02 

4.50506-64  1,727 16-01 

4.547,46-09  1.7  2956-01 

4.65  L4E.-.Q4 LJaiMzflf- 

9.85236-01 

9.R658E-01 

StimtxSd- 

9.89226-01 

2 46 
301 

_2i2£_ 

311 

• 

Figure  D-97.  Volume  Mattering  function  (sheet  3 of  3), 


D-228 


?5  jon  X97<‘i  urging 


lKg  3fri 


520  1 ° JUL  7 5 *055  $1A,CAT#17  - SEA  MZO  20M 


CTTT 

AN0L6  1 u e n ) 

nfitmi 

SIGMA  INSTR-Q 

HAN0*1 

rrmrar 

ANGLE  inECO 

oaT  A 

SIGMA 

1 

o .1"7To 

4.9S00E  Hi 

C 

0 , l? jo 

4.64646  6i 

2 

0.3500 

1 ,75008  02 

0 

0.3500 

1.61596  02 

3 

. O.TOUO.  ... 

6c?oone  oi 

0 

0,7000 

A 

10.00 

2 ♦ B 77 5t-0 L 

0 

10.00 

2.87756-01 

5 

15.00 

1.21506-01 

0 

15.00 

1.21506-01 

A 

20.00 

9.93176-02 

0 

20.00 

5,93 1TE-02 

1 7 

25,00 

3.25526-02 

0 

23.00 

3.26526-02 

L 

..  30.00 

1.90886-02 

0 

..  30.no ... 

1,90686-02  ... 

f 

40,00 

8.5219E-03 

0 

50.00 

B.5219B-0? 

to 

50.00 

5.20966-03 

0 

50.00 

5.20986-03 

u 

7.5 1286—03 

0.  ... 

. 60. no 

2. 61 2 SP -03 . „ 

12 

70.00 

1 ,59786-03 

0 

70.00 

I.S978B-03 

13 

60.00 

1.06666-03 

0 

so. 00 

1,06606-03 

1 4 

9 0.00  ... 

.7.90136-05 

0 ... 

90.00 

7 .4A 1 

15 

100.0 

6 » 8R70E-04 

0 

too.o 

a. 116706-04 

15 

1 10.0 

6.5173K-05 

0 

1 to  .0 

6.41736-04 

, „L7 

_ . 120.0 

4.239*6-0.6 

n 

170.0 

6.2396E-05 

lfl 

130.0 

6.51756-06 

0 

130.0 

6.41756-04 

19 

150.0 

6,93536-06 

0 

150.0 

6.93536-04 

2n 

140.0 

7 * 

...  0 

. 150.0  

3. *55599— 05 

21 

lAO.t) 

9. 08906-05 

0 

160.0 

9,08906-04 

2? 

170.0 

1.31966-03 

0 

170.0 

1 .3 19BS-03 

2 3 

...  1 

...  160.0 

1.47126-03 

ALPHA*  0.5067 

5/ ALPHA* 

0.722 

S*  .0,3.479. 

A /.ALPHA* 

0.276 

A.  0.14U 

B/S* 

0.012 

cpwR'gcTBD  acppi eemugCTiOM.o.enT 


ALPHA*  0.5133 iZik£M!L_JUIU 


5*  0.3675 
A*  0,1481 

A/ALPHA* 

S/S* 

U.2K7 

0.012 

Slf.MAI  0.0  DEGREES)* 
str.MAi  0.1  oeodiPsr* 
SLOPE ( 3 M 1 l L] RAD  1 • 

1069. 

815.  1 

“1.6.21 

S UP  TO  0.1  UBGREES* 

H. 95406 

-03 

NORMAL  1260*  2.436676-02 

. ML  . . 

RADIANS  

DEGREES 

•iijfl  IAN 

0.9991 

0,42246-01 

2.420 

(■CAN  1 

0.9541 

0.3040 

17.42 

VARIANCE 

0.1921  

MEAN  2 

0.147* 

8 .451 

RMS 

0.3641 

19.72 

RMS  2 

0.3109 

. .17.81 

t 

KAPPA* 

0.2017 

kappa  1 ■ 2. 

5774E-03 

r 

THITHta  BAP  3.92O5E-02  RADIANS**2 


,1 


Figure  D-98,  Volume  icattering  function  (iheet  2 of  3). 
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Figure  D-9B.  Volume  Kettering  function  (iheet  3 of  3). 
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Figure  D-99.  Volume  Kettering  function  (iheet  .1  of  3), 
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figure  D-102,  Volume  wittering  function  (sheet  3 of  3). 
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Figure  0108.  Volume  wittering  function  (theet  2 of  3). 
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Figure  D-108.  Volume  wittering  function  (iheet  3 of  3). 
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Figure  0)09.  Volume  Kettering  function  (iheet  2 of  3). 


0-263 


I 


1 


t 

t 

i 


t 


! 


i 

i 

!i 


I 


J 

fi 


j 


I 


73  JOM  1-7-  US'*.  10 


•TTfTTTjTI  IT  TfTtTT f STc  IfTTi  ft  -i  Fa"  m JiT 

ANGLE  ( i<  Mi  | Ai-.GLEIUEGI  SIGMA 

1.7*336-03  1.00006-01  b.TBOBE  02 


7**28 -oi  P,£  3*96-01 

3, *32*6-03  1.99538  -01 


TTTgrmnr 

*.70006  02 
5. <*69  16  02 

i;»nnr 

1.95*06  02 

7. 2*3*8  01 
9.271*26  01 

•1  .'Him  ft 

2.71931  01 
1.90036  01 

Tnnwrpr 

|^mn  jjtj 


10M 

INTEGRAL 

7.19196-03 
•"  1 

1.35*56-02 

2.12196-02 


V.Iha,  INTEGRAL 

1.  R '•<!  16-02 

~rrr  fin-j: — 

3.69026-02 
9.301 OE-02 

"HiTWBB'-M — 


1 

t r 
21 
31 


T’sloJI-Ol  1.16236-01 
*.|*jRi|6-Ojj  3.961^fi-0| 

ill  111 20 -02  A, 30966-01 
1.19**6-02  T.9-33H-01 

Tmrnrro — rvrragTTr 

2,19726-02  1.29696  00 

2. 76626-02  l ,96696  00 

■■■Ti'TTiyvE-oz 


-T.TnyT'^r* 

3.93936-02 

mm 

b, *2126-02 
7,91  qqu-ti; 

1. 139*6-01 
1.39966-01 

1.79216-01 
2.0*136-0 
Tv 
2.M**6«ri 


*r 


6.93*26-02  91 

1.79276-01  61 

1 • "a*"3-01 -') 

T7WTBl=ot  rrrr 

2.93*02,-01  111 

3.33336-01  121 

•nr-TTWffi — rrr 

<•.*7316-01  1*1 

sw=a-wMai m- 

6.32'<46“0I  171 

•mm  i m- 

7,6  2*76*  ,1 1 20  1 

waasft  j-ft 

9. ’7*16-01  21* 

JiHHEfr  Hi 

9,912*1-01  23). 

* ,*•» « » E-  )< l \n 

y.7ibfr^ui  Kt 


Trwnrw 

2.91196  00 


*.83*16  *R)2 


3.7-736-02  9«i'l  1*6  GO  2.030*6  DO 

mm 4 must 

1.7*936-01  1.000N6  01  *.l*3i!6«!)l 

mm  ■ mm mm 


L 44*4=4. 

3.13*26-01 

•i.**«*6-ni 

3,6  19S6-fll 


*.3*336-01  2.90006  01 

MDBHft WH- 


3.2*626-02 

-hW*t- 

S. 03296-03 

HI5*7S-0? 


3.719,96-01 

3.3  9106-01 
3.37**3-01 

xirraf^r 


*,93136-01 
7. 33* 0^-0 

* . 1 " 


*,00006  01 

*. 


si  -mm 


9.99936-01 

1.0*726  no 

3, 90006  01 

a. 00006  01 

3.06608-03 

2,2*7*5-03 

3.«12*H-I)1 

1.92***-.C1. 

9,765*6-0 1 

2** 

— iiii. . 

1.13*9ii  no 
1.22176 
1.30906  IV. 

*,90009  01 
7,00008  01 
.-7.40001  01. 

l.Hiei-oi 

1 .33131-03 
1.10JO6-C3 

3, <!'!*:*  i?-oi 
3.9*226-01 
2,9*f,tt»-01.  ... 

Ati . : c -5  r 

n.ticoi-ol 

<;>*» 

* ... 

1,3*416  ft 6 

1. *11336  00 
1,9706“  00 

F,  Moo  « ot 
6. 90 JOB  01 

6,4i*<JI!-4i 
7, *3036-0* 
«S.*OO«E-0< 

3 , 939*6  *01 
3.9*2*6-01 

MiM-M 

9.ACJ96-01 
°»  *9036-01 

mi 

27b 

. m .. 

TTtrrtt  ^ 

1.7*976  I", 

1 . 1 32  * i!  0 V 

i'.OOO'M:  (12 
1. "91)0“  02. 

*, 1 79AE-U* 
9 .7*796-0* 
9.  *02*8-0* 

S.MiiMP 

.1 .9*926-0  l 
. J.°72  lij-01 

-.  MS-11 

A.n07;:-Jl 

*. *i 4 . 

i,  r 

2vt 

?'l  •, 

TTSTS^ST  *n 
? .0071  £ 0" 

2. 0***6  -X 

1 . 10006  02 
1,19006  02 
1.20006  02 

TvSaiiS-t* 

D.19*T£-0* 
J , 1 ft  28  E-Q* 

3 , 97  30e*0 l 
J.977GE-01 
3 .“P026-01 

9.9,?1  6-0t 
J,92'  -6»')l 

301 

V 

3M 

rTTW^T 

2.2*396  On 
2.39*26  00 

t.Hwf  it 

1 .30003  0? 
1.35006  02 

9,2*996-0* 
9 .3*906-0* 

3. 4 R2SB-L '1 
3. ‘'0*96-01 
..  3.9*708*01 

• **9?8-'.'l 

lyg.-'U  . 

3 in 

3t  1 

32* 

tfuisrw" 

2.93076  00 
2.91606  00 

1, *00(11  of 
1.4900B  02 

9 .4t33E-0* 
3.97566-0* 
*.51696-0*  . 

T.««vr-(^r 

?, *9106-01 

Kwnt-TT~ 

9 , 96 1 98-01 
,:l. 

V3T 

3 3o 
3-1 

HTTOTTr 
2,79236  00 
. . 2,37931  no  . 

1**0006  02 

1.*  3.006  ..0-2 

ft.iS*ir-w 

1. 19276-03 
1.91*26-03 

?,*4*7?-ftl 
3 ,996(16-0  l 

^5.  Ml 'E-o  T 
9. 97* IE-01 

a*  b 
331 

,_JiO 

zTsSTTa  *i0  1.70606  0* 

3.09*36  00  1.73001  02 

i.iipoi-oi 
*.3*398-03 
. <r.eflTgftB.T.qJL— 

*l¥JT3FEr 

* , 0052E-01 

9. 99716-01 
1.OQ0OB  00 

•)bb 

_iii 

Figure  D-109.  Volume  Kettering  function  (sheet  3 of  3), 
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Figure  D-l  10,  Volume  wittering  function  (iheet  2 of  3). 
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Figure  D-l  10.  Volume  Mattering  function  (iheet  3 of  3). 
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Figure  0-111.  Volume  Kettering  function  (elieet  2 of  3). 
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Figure  D-  111.  Volume  Kettering  function  (sheet  3 of  3). 
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Figure  D-l  12.  Volume  Kettering  function  (iheet  3 of  3). 
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Figure  D-l  1 5 . Volume  scattering  function  (sheet  2 of  3), 


D7281 


25  .I'JN  1976  1200.19 


5io"2  3jUL75 
ANCLE ( UAH ) 
1,7*536-03 

2 i z!)  STA  .CAT>21  - 56  A HO'J 
ANGLEIOEGI  SIGMA 

1,00006-01  5.5*716  02 

2M 

INTEGRAL 

5.97316-03 

NORM,  INTEGRAL 
1.85126-02 

1 

? . 197z6"03 
2.7A62S-03 
3.482*6-03 

1725996-01 

1.58*96-01 

1,99338-01 

3.94896  02 
2.99996  02 

8,  i4Huttt“U3 
1. 27905-02 
1.76106-02 

2.?Sz2B-0t 

3.9A07E-02 

5.4576E-02 

1 1 - 
21 
31 

i 5.51926-03 

A. 9**36-03 

2.5114  f-7fl 
3,16236 -01 
3,9811  ML 

2.S4S9E  02" 
1.69196  02 
,1,2595?  02 

2.33*96-02 

3.01616-02 

...3,82*66-02 

T«l -0 2 
9.3473&-02 
1 • 18}  J(f“0  ^ 

41 

51 

61 

1,10126 -02 
t .3*548-02 

s.  a i Ilf -or 

6.30966-01 

7.9*336-01 

T;i3Sff  Jr 

7.04388  01 
5.28596  01 

*.TfT7J?-5j 

5.92326-02 

UUH-Ot 

1.83576-01 

2.25*16-01 

H 

81 

91 

1 » 7 A 531 
2.19728-02 
f 2.766.26-02 

i.Ooooe  bo 
1.25996  00 
1.59*96  00 

i.763it  Ul 

2,59276  01 
1.T27S6  01 

fl.aseSe-oz 

. 1.05756-01 
1.24416-01 

2.?2e5l:-01 

3.27746-01 

3.8556B-01 

101  1 
111, 
121 

3 .*"2*c  —.>* 
<■,39*18-02 
3.5192--02 

u^bjd  6u 
2.51196  00 

" r.Tnrr  jit 

7.08206  00 
*.<.0236  00 

1 »*31‘ ZS-U  1 
1.63516-01 
1.83096-01 

* .‘'*5  fit! -0 1 
5.0673K-U1 

5,673*6-01 

~m 

1*1 

m 

. oiro^r 

9.7*738-02 
V.  1,10126-01 

J.96T1  BTT8 
5.01196  00 
..  A .20966  00 

2.94*86  d)U 
1 .*1686  00 

2.C193I-01 

2.19836-01 

5.33326-01 

6.Z5C16-U1 
6.78176-01 
7.23101-01  . 

iaI 

171 
181  . 

1 « v^nie  -;ii 
1 ,7*? 36  -01 
2 , At  a(’6-01 

7,?*v:lc  )3 
l.onooe  oi 

1.50006  01 

2.3*356-01 

1.1*958-01 

2,*?  U8-'Jl 
2, 5717c -Cl 
2. 77166-01 

7.'lUlE-Ml 
7 , 97006-0 1 
8,59958-01 

Ivl 
2o  1 

Jlfb  ... 

S.*oaT6  -01 

*.36338-01 

5.23606-01 

i.WOhfl  01 

2.5000c:  01 
. 3,00006  01 

Tl’i'&'AST-OT 
2.7 1236-U2 
lv7S97?-02 

2.6937E-01 
2.97676-01 
- 3.03146-01 

W • ^ 834E-0 1 
9.2253E-01 
. 9.39406-01 

zii 

216 

221 

i.iOfUc-oI 

.S.9913e-01 

..7..,.  6.5  ftOft  -0L 

3.6000e  01 
*,00006  01 
vww,  ? 

i.o1o*b-5? 

*.28906-03 

5;5Tr5?-5t 

' ,09736-01 
}>  1.1.706.-0.1.. 

9.4l736-Ul 

9.59406-01 

9.6*018-01 

226 

231 

7..i.0 

*!  . 72  iit  -01 
9.59936-01 
1 .0*726  00 

s.oo66b  ul 

5.50006  01 
_ 6.00008.01... 

i ii-ol 

2. *9556-03 
1.9. 93?  e- 03 

2.1 aioS-o i 
3.14576-01 

...  3,15576-01. 

4.7i59b?-OL 
9,74886-01 
. 9.779*6-01 

Til — 

246 

25) 

1.13*3:  00 
1.22176  00 
1.30°(H  no 

S. 5 0008  01 
7.00006  01 
7.50006  01 

1 , i 1336-03 
t. 23 196-03 
1.00016-03 

3.  r«Vt-4i 

3.17096-01 

i.ilosil-oi 

4,82696-01 
.ti8*4  B2.ru?. 

Ut 

261 

1.39631  oo 

l.*830e  00 
1,57006  00 

(1.6M61  61 

9.50006  01 
. o . 00008  01. 

97T!Tl6-a* 

7. 0*308-06 
9,30776-0*. 

3.10  Isc-ol 
3. 1P36E-01 

^•§Ti9PE*“6 1 
9.87266-01 
9,88391-01 

in 

276 

281 

1.65  A1!  C-0 
1,7*336  00 
1 . <•'6 

<?!hzOr  «n 

l.  non  ns  02 
1.05006  .02 

5 .8531)6-0* 
5.8A**C-0* 
n.<.17  ?£-0-’ 

3 • lu  ZAc'-O l 
3.11376-01 

....  3.  ie? 66-.". 

9.89*28-01 

9.90386-01 

4.91293-01 

7 i!  A 
?*1 
790 

1.91996  00 
2.0071c  00 
2. 09**6  00 

1.10006  02 
1.19006  0? 

9.29306-0* 
5.1906  6-0** 
5,1(13  76-0*. 

3.201*8-01 

3.20*06-01 

3.20666-01 

9.92166-01 
9,9294  6—01 
9..4375E-01 

?01 

306 

311 

2 . 19 i 76  00 
2.26H9E  00 
2 • 35.62.8  n0 

tildOOS  02 
1 .3000 6 02 
! .3 3006... 02 

5.09028-0* 
3. 168*6-0* 
5. 36*56-0* 

3.20896-11 
3 . 2 1 ’ 1 E -0 1 
3.ai»2E-0> 

9.9**68-01 

9.95178-01 

9.95836-01 

3 16 
821 

S.AhSBI  ort 
2.53076  00 
2,61*08  00 

l ,*006662 
’.*3006  02 
1 .5000?  3?, 

4:H758-84 

6.03926-04 

6,5022.6-0* 

j.llsSB-ffr 

3,21726-01 
..  3 ,..2.19 1 E *0  L 

9,96461-01 

9.9VO66-01 

.9,9764.8-01. 

ill 

336 

341 

2 . 7 o J 3 6 00 
2.79258  00 
2.979P6  00 

1.55006  02 
1.60006  02 
'..630.16  .02 

7,1*636-0* 
8.13796-P* 
1. 0157«-O3 

3.22086-01 

3.2224E-01 

3.22396-C’, 

4,98176-01 
9. 98666-01 

4,o'H2E-01 

346 

351 

356 

2 .9*7  it  00 
3.05*36  00 
3.1*166  00 

1.70006  02 
1.75006  02 
1 .80006  02 

1.39908-03 

1.578*6-03 

1,62916-03 

3,27336-01 

3.22A3E-01 

3,22676-01 

9,49556-01 
9. 99886-01 
1.00006  00 

36’. 

366 

371 

PAUSE  9 6 A 0 V 

PLOTTER 

Figure  D>  1 15.  Volume  scattering  function  (sheet  3 of  3). 
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APPENDIX  E 


GENERAL  DESCRIPTION  OF  UPLINK  PROGRAM 


Program  UPLINK  is  used  to  reduce  the  experiment  uplink  data.  This  data  was  recorded  aboard  an 
aircraft  using  a receiver  looking  at  an  underwater  laser  pointed  upward.  The  receiver  data  for  each  received 
pulse  is  digitized  and  recorded  on  magnetic  tape.  A data  sample  from  one  of  these  tapes  is  shown  in 
table  E-l.  Program  UPLINK  reads  this  data  from  tape  along  with  control  information  from  cards.  The  tape 
is  searched  for  a specified  rnn  number.  The  data  for  this  run  number  is  read  in,  filtered,  calibrated,  and 
plots  of  radiance  loss  versus  zenith  angle  are  generated. 

The  data  is  calibrated  using  the  values  for  receiver  sensitivity  and  laser  power  given  in  Volume  II, 
Section  S.  The  program  initially  assumes  that  the  receiver  was  operated  at  its  most  sensitive  range.  This 
value  is  then  multiplied  by  100  if  switch  5 is  zero,  indicating  that  the  1 percent  transmission  Alter  was  in 
place,  and  by  10  if  switch  4 indicates  that  the  10  percent  transmission  filter  was  in  place. 

The  data  for  one  run  are  read  and  the  receiver  output,  roll  angle,  pitch  angle,  and  value  of  switch  1 2 
are  saved.  Switch  12  is  an  indication  of  whether  or  not  the  receiver  was  getting  sufficient  signal  for  auto- 
matic tracking.  When  all  data  for  one  run  have  been  read  from  tape;  the  signal,  roll  angle,  and  pitch  angle 
are  Altered.  The  signal  is  multiplied  by  the  calibration  factor,  and  the  zenith  angle  is  calculated  from  the 
roll  and  pitch  angle,  Two  plots  of  the  signal  are  generated  (see  figure  E-l  A and  figure  E-IB);  one  is  the 
uncorrected  signal  and  the  second  is  multiplied  by  the  secant  squared  of  the  zenith  angle  to  correct  for  vari- 
able range  to  the  aircraft.  The  plots  are  drawn  aide  by  side  on  an  11  inch  by  17  inch  page. 

The  digital  Alter  applied  to  the  data  is  deflned  as  follows.  Let  the  digital  sequence  to  be  Altered  be 
denoted  by 

{sj},  1*  1,  n.  (E-l) 

then  deAne 

uj  - s, 

U1  " ul-l  "N^i-T'i5  2*“  (E*2) 

where  A'  is  a constant  wilh  N > 1.  This  is  essentially  a low  pass  filter  with  a time  constant  of  (N-l)  At 
where  At  is  the  sample  interval  for  the  original  sequence  (Sj),  The  sequence  (uj)  is  a smoothed  and  time 
delayed  Image  of  {sj}. 

Now  define 
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TABLE  E-l.  UPLINK  PROGRAM  DATA  SAMPLE. 
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22  JUL  1975 
RUN  NO.  12 

DEPTH  - 9.1  METERS  ( 30  FEET  1 
A/C  ALTITUDE  - 610  METERS  ( 2000  FEET  ) 
LRSER  ANGLE  « 42.5  DEGREES 
PERCENT  TRACKING  - 74. 6 _ 

g PEAK  « 4.2  X 10"8  METER-8  AT  54.3  DEGREES 
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Figure  B-l  A,  Uncor routed  uplink  plot. 


( ,V w* tv« ! I Wf  1 Vrt MlWi  i tt-y-u 


> ij  I’Afl  i4  W WiJrtif  /i 1 


LOSS,  OB 

-55.00  -50,00 


* • 
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Tills  is  really  the  same  filter  applied  in  the  opposite  direction,  so  that  { v } has  a time  advance  rather 
than  a time  delay.  Finally,  define 

, *i  ■ j(Uj  + Vj)  1-1,  n.  (E-4) 

Then  {zj  > has  neither  an  advance  or  delay.  Any  desired  degree  of  smoothing  can  be  accomplished 
with  this  filter  by  choosing  different  values  of  the  filter  constant,  N,  Figures  E-2A  through  E»2C  show  the 
results  of  filtering  the  same  raw  data  with  various  values  of  N, 

INPUT  DATA.  The  raw  data  to  be  processed  is  read  from  magnetic  tape.  Control  information  is 
read  from  cards  using  NAMELIST  with  name  INPUT.  Variables  to  be  input  from  cards  are  defined 
as  follows. 

SIZE  - A real  variable  allowing  plots  of  different  size  to  be  generated.  See  DOWNLINK  for 
details.  Default  value  is  1,0. 

NAVG  - An  Integer  variable  specifying  the  filter  time  constant.  Same  as  N in  the  discussion  of 
the  filter.  Default  value  is  NAVG  ■ 6 which  gives  a filter  time  constant  of  0.25  seconds. 

RUN  - An  Integer  variable  specifying  the  run  number  of  the  data  to  be  processed.  The  tape  is 
searched  until  the  proper  run  number  is  found.  Data  is  then  read  in  and  processed  until 
the  run  number  changes. 

DATE  - An  Integer  array  of  dimension  3.  Used  for  the  date  as  a character  string;  i.e.,  DATE  - 
‘24  JUL  75’.  Default  value  is  all  blanks. 

LASANG  - A real  variable  specifying  the  angle  of  the  laser  in  degrees  from  the  zenith. 

0 < LASANG  < 48.75  (the  critical  angle), 

DEPTH  - A real  variable  specifying  the  dopth  of  the  laser  in  meters. 

ALT  - A real  variable  specifying  the  altitude  of  the  receiver  (aircraft)  in  feet. 

OUTPUT  DATA.  All  output  from  the  program  is  either  to  the  plotter  or  to  the  printer.  The  track* 
ing  signal  is  plotted  on  the  graphs  just  above  the  zenith  angle  axis. 

EXTERNAL  SUBROUTINES  REQUIRED; 

AXISM,  COPY,  HEDING,  PHIL1NE  - See  discussion  of  these  routines  in  program  DOWNLINK. 

OPSAT  - A PL/I  subroutine  to  read  and  decode  the  raw  data  tapes.  For  each  coll  to  OPSAT,  one 
logical  record  is  returned,  A logical  record  consists  of  the  time,  the  settings  of  the  12 
switches,  and  the  values  of  the  6 data  channels.  OPSAT  also  returns  the  length  of  the 
block  read  from  tape  and  flags  for  end  of  record,  end  of  data  set,  and  error. 

Figure  E*3  represents  a program  listing  on  the  following  14  pages  of  the  SATCOM  UPLINK  data 
reduction. 
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Figure  E-2A.  Uplink  data  filtered  with  time  comtant  » 0 woondi. 
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Flyura  B-2B.  Uplink  data  fllttnd  with  thna  Qoutul  ■ AS  raoondi. 
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Figure  E-2C,  Uplink  data  filtered  with  time  conitint  ■ ,95  Mcomli, 
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Ftpite  E-3.  SATCOM  data  rcdnctMO  program  listing  (sheet  5 of  I4>. 


<# 

If 

$ 


I 

i 


it 

« 

* 

t 

* 


0 

1 


J 


i 


n/ii  * ♦ . 

**V*  •♦ 

**  |h  • 

xi.  ^-*2 

-*  II  Xj  00  u 
X ti'+W'*  - 

^S^-oSn 

li  «B  -u.*-* 


m 


?.i 


*a  m\ 

"*l>  *U-JJ  > 


rg-j£^*t«a 


ouuo 


U..iQtLLK^U.*qi 

jMtk.O *0iJt'  “-■“** Co 

: 3i3B 


**te  2"teJ3  - 

s'*it«*i  rfesB-is 

HdSi  s 


<9  o 

$ ! 


§§ 


SSWUR  *8 


du  Hh««or 
UMwaiCixir  c 


I 


ns 

P*5of"S“ 


B'S  B 


552*5 


aisaa^l 


»3<2; 


E-14 


uilMU  ^l/^.olt.d.llUJliJ'ij'ilUttrtS^x  t&SkXUlha!’ W*UMliiliLl&U\A  i/CidLuAAib,]',.,,!.,  ,n„w/, ■.<,..  -.-. 


a.iiui.:ai  mamma  tw  u I'liiifilfM^ili^liatflilfl 


Figure  E-3.  SATCOM  upfat  data  redaction  pmpm  Kithig  (sheet  6 of  14). 


RgmeE-3.  SATCOMopfakdai*  redaction  program  listing  (sheet  7 of  14). 
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Figure  E-3.  SATCOM  upSak  data  redaction  program  fifing  (sheet  10  of  14). 
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E-3.  SATCCM  upink  data  ictiuclun  progm  fatiag  (slut  13  of  14). 


APPENDIX  F 


GENERAL  DESCRIPTION  OF  DOWNLINK  PROGRAM 


Progrum  DOWNLINK  is  used  to  reduce  the  SATCOM  downlink  data  which  was  recorded  In  the 
quick  scan  mode.  In  this  mode  of  operution,  a 60  by  60  point  scan  was  made  or  the  photo  cathode  surface 
and  the  resulting  3,600  points,  plus  seven  words  of  additional  information  were  recorded  on  magnetic  tape 
in  one  record.  At  each  depth,  a specified  number  of  scans  (records)  were  written  onto  the  tape  with  the 
last  record  followed  by  a tape  mark,  The  first  record  of  each  file  contains  the  seven  words  of  environmental 
dutu.  The  original  raw  dutu  tupes  were  copied  to  prevent  loss  of  datu  in  the  event  of  dumnge  to  u tape.  All 
tiles  from  a given  original  tupc  were  written  Into  a single  file  of  the  copy.  In  addition  to  the  original  datu, 
each  record  of  the  copy  contains  four  words  of  information  giving  the  tape  number,  file  number,  record 
number,  ond  number  of  data  words  from  the  original  tape, 

Program  DOWNLINK  reads  certain  control  Information  from  cards;  included  ir  this  are  the  tape  and 
tile  number  to  be  processed.  The  input  tape  is  then  searched  for  tills  data.  All  datu  records  from  the 
Indicated  tape  and  file  are  read  in  and  uverngod,  This  averaged  data  is  then  calibrated  und  contour  plots  of 
the  Image  are  generated. 

Calibration  of  the  dutu  follows  the  procedure  outlined  in  Volume  II,  Section  C.  Let  Sx  y be  the 
ruw  dutu  value  recorded  on  tape  for  cumeru  coordinates  (x,  y).  Then  the  calibrated  value  ZX(y  Is  calculated 
us 


Z . (S>.ytl4)(60°)  (93)  iN 
'■x.y  r:  u d a 


zx,y  - 0 


a CJD  F px,y 
for  Sx>y  ■ 0, 


for  S*  w 0 


Where, 


(F-l) 


G is  the  camera  gain  factor, 

Cijj  is  the  dwell  gain  fuctor, 

Gd-  1 + 11.23D1-036, 

D Is  the  number  of  turns  on  the  dwell  potentiometer, 

F is  the  surface  flux  In  watt-em”2 

is  a unit  spectral  radiance  factor  as  defined  in  Volume  II,  Section  3.  It  has  the  value 
AN* « 2.68  X 10“8  r2  < 898704  (F-2) 


P-1 


AN*  • 2.68  X l<T8  + 2X  10‘l0^ty~  - 45^  r2  > 89B7Q4, 
where  r“  * x“  + y*-  and  x,  y ur»;  the  camera  coordinates  of  each  point. 


The  factor  6Q0/PXty  corrects  for  the  variable  response  over  the  photo  cathode  surface,  PX)y  is  a 
least  squares  fit  to  the  response  of  the  photo  cathode  surface  with  uniform  illumination. 

The  factor  93  accounts  for  the  93  A bandwidth  of  the  filter. 

The  term  14  added  to  the  raw  data  corrects  for  the  dc  bias  in  the  raw  data. 

The  calibrated  data  is  searched  to  locate  the  peak,  then  Integrated  over  circles  centered  on  the  peak. 
These  circles  ere  opened  up  from  5s  to  90°  In  5s  steps.  The  resulting  values  are  plotted  to  give  integrated 
loss  vs.  field  of  view.  The  values  are  also  output  in  tabular  form  on  the  printer. 

A sample  plot  generated  by  this  program  is  shown  in  figures  P-1  A and  F-1B  A sample  of  the 
printed  output  appears  at  the  end  of  the  program  listing. 

INPUT  DATA.  The  data  to  be  processed  is  read  from  magnetic  tape  as  described  above.  Control 
information  is  read  from  cards  using  NAMELIST  with  name  INPUT.  The  variables  to  be  input  from 
cards  are  as  follows: 

TAPE  - An  Integer  variable  specifying  the  number  of  the  original  raw  data  tape. 

PILE  - An  integer  variable  specifying  the  file  number  on  the  original  raw  data  tape. 

GAIN  - An  integer  variable  specifying  the  camera  gain  for  the  data  to  be  processed.  This  is  a 

coded  value  and  must  be  either  0.  1*2,  or  3. 

DWELL  - A real  variable  specifying  the  dwell  setting  for  the  data  to  be  processed.  This  is  given  as 
the  number  of  turns  of  the  dwell  potentiometer.  0 < DWELL  < 10. 

FLUX  - A real  variable  specifying  the  surface  flux  at  the  timo  the  data  was  recorded.  Units  are 
watt*cm“2. 

SIZE  - A real  variable  which  allows  plots  of  variable  size  to  be  generated.  The  default  value  is 
SIZE  ■ 1.0.  Tills  produces  plots  on  an  1 1 Inch  by  17  Inch  pBge.  SIZE  ■ 0.5  would 
reduce  the  plots  to  a 5 1/2  Inch  by  8 1/2  inch  page.  If  SIZE  Is  given  a value  greater  than 
1.0,  then  a request  must  be  made  for  30  inch  plotter  paper.  SIZE  should  not  be  given  a 
value  grater  than  2.72  since  this  would  attempt  to  generate  a plot  page  greater  than  30 
inches  high. 

C TIME  - An  integer  array  of  dimension  5 which  may  be  used  to  correct  the  tape  time  for  cases  in 
which  the  clock  was  in  error  when  the  data  were  recorded.  Default  values  for  this  array 
are  all  zeros  which  assumes  the  clock  was  correct. 

C TIME  ( 1 ) ■ day  correction 
C TIME  (2)  ■ hour  correction 
C TIME  (3)  ■ minute  correction 
C TIME  (4)  ■ seconds  correction 
C TIME  (5)  ■ milliseconds  correction 

The  correction  is  added  to  the  tape  time  so  that  the  correction  for  a clock  1 hour  slow 
would  be  C TIME  - 0,  1, 0,0,0. 

DEPTH  - A real  variable  specifying  the  depth  of  the  camera  ip  meters. 


F-2 


When  processing  is  completed  for  one  set  of  input  data,  another  set  is  read,  and  processing  con- 
tinues. Any  number  of  sets  of  input  data  may  be  stucked  for  a given  job.  No  special  code  is  required  for 
the  last  set;  the  progrum  terminates  normally  when  no  more  curds  are  found.  Since  the  input  tape  is 
rewound  ufier  each  set  of  data  is  processed,  it  is  not  necessary  to  process  the  datu  in  the  order  in  which  it 
appears  on  the  tape;  and  if  desired,  the  same  data  may  be  processed  more  than  once. 

OUTPUT  DATA.  All  output  from  the  program  is  either  to  the  plotter  or  the  printer  and  should  be 
self  explanatory  Note  that  the  printer  output  and  plots  may  be  matched  by  the  ID  number  which  appears 
along  the  left  hand  edge  of  each  plot  and  at  the  top  of  each  printer  page.  This  ID  number  is  merely  the 
date  and  time  at  which  the  job  begs.i. 

EXTERNAL  SUBROUTINES  REQUIRED. 

ALMNAC  - This  subroutine  calculates  the  solar  declination  and  the  equation  of  time  for  the  date 
and  time  given.  This  information  is  used  by  the  main  program  to  calculate  the  zenith 
angle  and  azimuth  of  the  sun  at  the  time  the  data  were  recorded. 

AXISM  - This  routine  generates  the  coordinate  axes  for  the  plots.  It  is  a slightly  modified 
version  of  the  CALCOMP  routine  AXIS. 

BCDBIN  - This  routine  converts  the  tape  time  from  three  words  of  binary  coded  decimal  infor- 
mation to  five  Integer  words  giving  time  in  days,  hours,  minutes,  seconds,  and  milli- 
seconds. It  also  produces  a character  string  giving  date  and  time  of  day  for  easy  output 
to  the  plotter. 

CONTUR  - This  routine  generates  the  contour  plots  from  the  60  X 60  array  of  tape  data. 

COPY  - This  routine  copies  the  input  cards  to  the  printer  and  to  a scratch  data  set  for  subse- 
quent reading  by  the  main  program. 

HEDING  - This  routine  keeps  track  of  lines  output  to  the  printer  and  prints  a new  heading  at  the 
top  of  the  next  page  whenever  the  current  page  is  filled. 

PHL1NE  - Tills  routine  generates  the  line  plot  of  Integrated  loss  versus  field  of  view,  It  is  a 
modified  version  of  the  CALCOMP  routine, 

Figure  F-2  represents  a program  listing  on  the  following  21  pages  of  the  DOWNLINK  dato 
reduction. 
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Figure  F-2.  Downlink  data  reduction  program  Biting  (ieet  6 of  21). 


Figure  F-2.  Downlink  data  reduction  program  listing  (sheet  7 of  21). 
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Figure  F-2.  Downlink  data  reduction  program  listing  (sheet  IS  of  21) 
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Figure  F-2.  Downlink  data  redaction  pragom  fisting  (Aeet  16  of  21). 
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figure  F-2.  DombSi*  date  redacttotp^of^  fating  (iwet  18  of  21). 
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Enpuc  F-2.  DovnfinkdaU  reduction  program  Estiog  (sheet  19  of  21). 
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RgpreF-2.  Doniiak  data  icdoctioa  ptognm  fisting  (Aret  20  of  21). 
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Rgaic  F-2.  Dowrfak  data  iedncoaa  program  fisting  (iheet  21  of  21). 


APPENDIX  0 


GENERAL  DESCRIPTION  OF  SLICE  PROGRAM 


Program  SLICE  Is  similar  to  DOWNLINK  and  processes  the  same  data.  The  difference  is  that 
instead  of  a contour  plot  of  the  data  array,  SLICE  generates  a plot  of  a cross  section  of  the  contours  along  a 
specified  azimuth  angle.  The  data  are  read  from  tape  and  averaged  exactly  as  in  DOWNLINK.  The  azimuth 
of  each  point  in  the  array  is  calculated  and  compared  to  specified  limits.  If  it  falls  within  these  limits  the 
value  is  calibrated,  its  zenith  angle  calculated,  and  the  calibrated  value  and  zenith  angle  are  saved.  If  the 
azimuth  of  the  point  falls  outside  the  specified  limits  it  is  ignored.  When  all  points  have  been  examined, 
those  saved  are  sorted  into  increasing  zenith  angle  order  and  a plot  of  radiance  loss  vs.  zenith  angle  is 
generated. 

It  is  necessary  to  specify  limits  around  the  desired  azimuth  angle,  since  it  is  probable  that  no  point 
will  have  exactly  the  desired  azimuth.  A more  elegant  method  would  have  been  to  utilize  an  interpolation 
scheme  to  obtain  values  at  exactly  the  required  azimuth,  but  this  seemed  to  be  an  unnecessary  complica- 
tion. Azimuth  limits  of  any  width  may  be  specified,  but  experience  suggests  that  the  desired  value  ±2  1/2 
degrees  produces  good  results. 

A sample  plot  generated  by  the  program  is  shown  in  figure  G-l  and  a sample  of  the  printed  output 
appears  at  the  end  of  the  program  listing. 

INPUT  DATA.  The  camera  data  to  be  processed  is  read  from  magnetic  tape  as  described  for  pro- 
gram DOWNLINK.  Control  information  is  read  from  cards.  All  variables  defined  for  DOWNLINK  also 
apply  to  program  SLICE. 

In  addition,  the  following  variables  are  required: 

AZL1 , AZL2,  AZL3,  AZL4  real  variables  specifying  the  azimuth  limits  in  degrees.  The  restrictions 
are  AZL1  < AZL2,  AZL3  < AZL4,  and  0 < AZL1,  AZL2,  AZL3,  AZL4  < 360.  Example:  to  obtain  a 
cross  section  along  the  azimuth  2 1 S degrees,  specify 

AZL1  ■ 212.5,  AZL2-  217.5,  AZL3-  32.5,  AZL4  - 37.5, 

The  values  AZL1  - 32.5,  AZL2  - 37.5,  AZL3  - 212.5,  AZL4  ■ 217.5,  would  produce  the  same 
results  except  that  the  resulting  cross  section  plot  would  be  reversed  left-to-right. 

As  a general  rule  use, 

AZL1  - A-2.S,  AZL2  - A +2.5,  AZL3  - AZL1  ± 180,  AZL4  ■ AZL2  ±180.. 

All  cards  are  read  from  NAMELIST  name  INPUT.  Any  number  of  sets  of  data  may  be  stacked 
for  a given  job.  Processing  continues  until  no  more  cards  are  found. 

OUTPUT  DATA.  All  output  is  either  to  the  printer  or  the  plotter.  Examples  have  been  cited  above 
and  should  be  self-explanatory, 
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EXTERNAL  SUBROUTINES  REQUIRED.  Same  us  lor  program  DOWNLINK  except  for  CONTUR 
which  is  not  required. 

Figure  G-2  represents  a program  listing  on  the  following  14  pages  of  the  SLICE  datu  reduction. 


. SLICE  data  reduction  prtpan  luting  (sheet  2 of  14). 


Figure  G-2.  SLICE  data  redaction  program  Kiting  (sheet  3 of  14). 


Figure  G-2.  SLICE  dxU  redDctHxr  program  fisting  (street  4 of  14). 


figure  G-2.  SLJCE  dat*  reduction  program  fisting  (sheet  5 of  14). 
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FjpneG-2.  SLICE  data  redaction  program  fating  (sheet  10  of  14). 


figure  G-2.  SLICE  data  redaction  pnyau  Itxling  (sLrct  12  of  14). 
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APPENDIX  H 

GENERAL  DESCRIPTION  OF 
AUTOMATIC  HEMISPHERICAL  SCAN  PROGRAM 


This  program  is  similar  to  DOWNLINK,  but  is  used  to  reduce  SATCOM  data  recorded  in  the  full  scan 
mode.  In  this  mode  of  operation,  the  photo  cathode  was  divided  into  three  different  regions  which  were 
scanned  with  a different  point  density  in  each  region.  Figure  H-l  illustrates  the  scan  areas.  The  center 
squuro  was  sampled  with  AX,  AY  ■ 20,  the  outer  square  with  AX,  AY  • 80,  Since  the  contour  plotting 
routine  requires  a rectangular  array  of  duta  with  constant  AX  and  constant  AY,  it  Is  necessary  to  interpolate 
values  to  produce  u constant  AX,  AY  * 20.  The  result  is  a 185  by  185  point  array.  Points  outside  the 
circular  urea  actually  scanned  are  set  to  zero. 

Except  for  the  differences  in  bookkeeping  to  read  the  data  and  the  interpolation  scheme  required, 
the  program  is  the  same  as  DOWNLINK.  The  same  input  is  required  and  the  same  output  is  produced.  A 
sample  plot  is  shown  in  figures  H-2A  and  H-2B,  and  the  printer  output  appears  at  the  end  of  the  program 
listing. 

INPUT  DATA.  Same  as  program  DOWNLINK. 

OUTPUT  DATA,  Same  as  program  DOWNLINK. 

EXTERNAL  SUBROUTINES  REQUIRED.  Automatic  Hemispherical  Scan  uses  all  the  subroutines 
used  by  DOWNLINK  plus  the  following; 

DENNIS  - This  subroutine  handles  the  bookkeeping  required  in  reading  the  camera  data  from 
tape  unci  loading  it  into  the  data  matrix  in  the  proper  location.  There  arc  actually 
three  versions  of  this  routine  to  handle  three  different  modes  of  recording  the  camera 
data,  In  some  cases  of  data  recording,  each  of  the  three  regions  of  the  image  were 
scanned  four  times;  once  at  each  of  the  four  possible  camera  gain  settings.  Thus,  the 
data  are  recorded  as  four  scans  of  region  1 followed  by  four  scans  of  region  2,  followed 
by  four  scans  of  region  3.  The  version  of  DENNIS  for  this  case  will  retain  for 
processing  only  one  scan  of  each  region.  The  scan  processed  is  the  highest  gain  which 
did  not  cause  saturation  of  the  A-D  converter,  If  the  lowest  gain  saturated  the  con- 
verter, that  data  is  kept  and  a message  is  written  to  the  printer  to  warn  of  the  satura- 
tion condition,  At  other  times,  only  two  scans  were  made  of  each  region,  and  at  still 
other  timeR,  only  one  scan  of  each  region.  Separate  versions  of  DENNIS  handle  each 
of  these  cases  and  some  care  is  required  to  Insure  that  the  correct  version  i«  used  for 
the  data  to  be  processed. 


H-l 
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1NTERP  - This  routine  calibrates  the  data  read  in  by  DENNIS,  and  interpolates  values  in  the  two 
outer  regions  so  that  the  resulting  data  array  has  the  same  density  of  points  in  these 
regions  as  in  region  1.  Calibration  follows  the  procedure  described  in  DOWNLINK. 

UWCAM  - This  is  a PL/I  subroutine  which  reads  the  raw  data  tape,  checks  for  errors  and  deter- 
mines the  length  of  each  record  read.  The  data,  the  length  of  the  record,  an  error  flag, 
and  an  end  of  data  set  flag  are  returned  to  the  calling  program. 

Figure  H-3  represents  a program  listing  on  8 pages  of  the  Automatic  Hemispherical  Scan  data 
reduction. 
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Figure  H-3 , Automatic  hcmiyhericai  scan  data  reduction  prognau  listing  (ahcct  4 of  8). 
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FjpueH-3.  Antnrotir  hemupheiicai  acm  data  redaction  program  listen  (Aeet  5 of  8). 
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APPENDIX  I 


GENERAL  DESCRIPTION  OF  DOWNMODEL 


This  program  is  used  to  evaluate  two  of  the  Integrals  in  Appendix  A.  The  first  section  of  the  pro- 
-am calculates  and  stores  functions  which  are  common  to  both  integrals.  The  different  integrals  are  evalu- 
ated in  subsequent  sections  of  the  program.  Evaluation  of  either  or  both  integrals  is  selected  by  program 
input.  The  first  of  the  integrals  which  may  be  evaluated  is 

I(x,y,z,7x.7y)-fl  +f[  ((l-e“«-A)AI  + AAI|p^Jdx0.dy0  W 

where  AI  is  defined  in  equations  (43)  through  (46)  of  Appendix  A. 

The  term  f|  is  given  by 

1 T fr’x-?*)2  + (7y-^>21 

2ir(fl}+^)  [ v x y 2(oy+*0)  J 


where, 


ai"[(i)  +(t)  *(' 

ej  is  the  angular  diameter  of  the  source, 

«2  i»  the  angular  diameter  of  the  receiver  field  of  vlow, 


■ variance  of  source  distribution  at  the  surface. 


The  term  A is  defined  as 


2.3 1 


A ■ ,255(sz  )exp 


itiij , ki.VitmA  I .;£uj^.u  SMUlUjj 


The  notation  means  that  this  integration  is  performed  with  A and  AI  evaluated  at  an  effective 

value  of  0 2.  This  is  accomplished  by  multiplying  the  original  62  by 


expfM'/y/  1 + 9/sz’)-  1 
exp(sz  '/V  1 + fii/sz')  - 1 


(1*5) 


The  calculations  are  further  generalized  to  allow  the  water  properties  a,  s,  and  0^  to  be  functions  of  the 
depth.  This  is  done  by  dividing  the  water  into  n layers  of  arbitrary  thickness.  The  top  of  the  i™  layer  is  at 
depth  Z|,.|  and  the  bottom  of  this  layer  is  at  depth  zp  The  waterin  this  layer  is  assigned  a scattering  coef- 
ficient sp  an  attenuation  coefficient  Bp  and  a variance  in  f(0)  of  Of,  z'  is  defined  by 


z'  - [z l + (x  - x0)2  + (y  - y0)2] 1 12  . 

Then  the  terms  uz',  sz\  (a  + s)z\  sPz',  and  sPz'^  in  the  original  equations  are  replaced  by  the  fol- 
lowing terms  respectively 


n 


where 

dZj » zj  - zpj  . (1-10) 


■v  *• 
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* t 


Also  evaluated  by  this  program  is  the  integral 


ll\.y,z,fl)“  //>-“«  A)A1*  + AAI*L  dx0dy 


(Ml) 


where  Al*  is  the  integral  In  equation  (52)  of  Appendix  A.  In  both  integrals,  the  region  of  Integration  R is 
defined  by 


- 5 ii  in  < x0  < S/3  zn 

- 4/3  < y0  < 4/3  zn  . (M2) 

Tlte  program  has  the  option  of  blanking  out  a region  in  the  x0  • y0  plane  corresponding  to  the  location  of 
the  barge  used  in  the  experiment.  The  shadow  of  the  pipe  holding  the  camera  platform  will  ulso  be  con* 
sidered  if  this  option  is  selected. 

As  part  ot  the  input  data,  the  program  requires  the  date  and  time  (PDT).  This  information  Is  used 
to  compute  the  azimuth  and  zonith  angle  of  the  sun,  and  the  x0,  y0  coordinate  system  is  oriented  so  that 
the  +xq  axis  is  directed  toward  the  sun.  The  known  orientation  of  the  barge,  camera,  and  sun  are  used  to 
define  a coordinate  transformation  from  x0  • y0  to  x”  - y"  with  the  x"  - y"  coordinates  having  the  origin 
at  the  comer  of  the  barge,  and  -x"  axis  along  one  end  with  -y"  axis  along  one  side  of  the  barge,  This  trans- 
formation is  used  to  simplify  the  determination  of  whether  or  not  points  in  R fall  on  the  barge  'or  not.  The 
relationship  of  the  barge,  camera,  and  various  coordinate  systems  are  Illustrated  in  figure  M . 


INPUT  DATA.  All  data  for  the  program  are  read  from  cards  using  NAMELIST  name  INPUT. 

The  following  variables  are  defined: 

ZL  -a  real  urray  of  dimension  50  giving  the  depths  in  meters  of  the  water 

layers,  ZL(I)  is  depth  of  top  of  layer  I.  ZL(I  + 1)  is  the  depth  of  the 
bottom  of  the  layer. 

AL  -u  real  urray  of  dimension  50  giving  the  scattering  coefficient  in  each 

layer  of  water 

S - a real  array  of  dimension  50  giving  the  scattering  coefficient  in  each 

layer  of  water  __ 

THBRSQ  - a real  urray  giving  the  vulue  of  0^  in  each  layer  of  water. 

LAYERS  - the  number  of  values  of  ZL,  i.e,,  1 more  than  the  number  of  water 

layers. 

NXGD  - un  Integer  vuriable  giving  the  number  of  grid  points  in  R parallel  to  the 

x0  uxls.  Default  value  Is  NXGD  ■ 90.  Values  less  than  90  may  be  input, 
90  is  the  maximum  value  unless  the  program  is  modified  to  increase  the 
dimension  of  ail  arrays  that  depend  on  tills  number. 

NYGD  number  of  grid  points  in  R parallel  to  the  y0  axis,  Default  (und  maxi- 

mum) vulue  is  NYGD  ■ 60. 

VARR  - u real  variable  specifying  the  value  of  varfR), 
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AZIMUTH  OF  SUN  - <* 


Figure  M . Coordinate  system  for  program  DNMODEL. 
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F0V1EW 

GAMP1.GAMP2,  STEF1  - 

GAMP3,  GAMP4,  STEP2  - 

THZERO 

GPY 

BARGE 

TIME 


DEL1.DEL2,  STEPD  - 


BETAX,  BETAY 


XPEAK,  YPEAK 


IRPROF 


INTPWWR 


a rea’  variable  specifying  the  angular  diameter  in  degrees  of  the  detector 
field  of  view.  Default  value  is  1°. 

real  variables  specifying  the  range  and  increment  for  7'  in  the  evaluation 
of  equation  (1-1).  7'x  takes  values  from  GAMP  l to  GAmP2  in  steps  of 
STEP1.  All  units  are  degrees.  Default  values  are  -90,90,  and  10. 

real  variables  which  allow  a different  step  size  in  7'  over  a subinterval  of 
the  range  GAMP1.GAMP2. 

a real  variable  specifying  a value  in  degrees  for  0O,  the  variance  of  the 
source  distribution  at  the  surface. 

a real  variable  specifying  a value  for  7y  in  degrees. 

a logical  variable  specifying  whether  or  not  to  consider  the  barge  in  the 
calculations.  BARGE  ■ .FALSE,  will  cause  the  barge  to  be  ignored. 
Default  value  is  .TRUE. 

An  integer  array  of  dimension  four  for  specifying  the  date  and  time. 
TIME(l)-  month, 

TIME(2)  * day  of  month, 

TIME(3)  ■ hour  of  day  (PDT), 

TIME(4)  - minute. 

real  variables  specifying  the  range  and  increment  for  values  of  A in  the 
evaluation  of  equation  (M 1).  Units  are  degrees  and  default  values  are 
5, 90,  and  5. 

real  variables  specifying  values  for  the  terms  (1  - n/n')Rx  and 
(1  - n/n')Hy  which  appear  in  the  evaluation  of  equation  (1-1 1).  See 
equation  (f>3)  of  Appendix  A. 

real  variables  specifying  the  values  of  7X  and  7I,  at  which  the  irradiance 
distribution  has  its  peak.  This  point  is  the  origin  for  A in  equation 
(1-1 1).  Default  values  are  0, 0 if  equation  (1-1)  is  not  evaluated.  If 
equation  (1-1)  is  evaluated,  tire  actual  location  of  the  peak  is  used  and 
any  values  read  in  are  ignored. 

a logical  variable  specifying  whether  or  not  to  evaluate  equation  (1-1). 
Default  value  is  .TRUE. 

a logical  variable  specifying  whether  or  not  to  evaluate  equction  (1-1 1). 
Default  value  is  .TRUE, 


OUTPUT  DATA.  A sample  of  the  printed  output  from  this  program  follows  the  program  listing. 
The  first  page  is  a copy  of  the  input  catjh.  Page  2 is  a list  of  parameters  for  the  first  set  of  input  data,  and 
includes  the  depth  profiles  of  a,s,  and  (fZ.  Page  3 contains  the  results  of  evaluating  equation  (1-1).  The 
columns  labeled  GPX  and  GPY  are  the  values  of  7I  and  7 y.  The  column  labeled  SUM  is  the  value  of  equa- 
tion (1-1)  at  this  point.  The  column  labeled  FI  is  the  value  of  f j equation  (1-2).  The  column  labeled  SUM1 
is  the  term 


if  (l-e“8Z'-A)Aldx0dy0, 


’.v ..  ;v., 
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and  the  column  labeled  SUM 2 is  the  term 


J/AAIi;3.5y,od!'°'  <M4> 

The  columns  labeled  zenith  and  azimuth  are  just  y'xl  7!,  converted  to  polar  coordinates  with  the  000° 
azimuth  taken  in  the  Xq  - Z plane.  Page  4 of  the  printed  output  is  the  result  of  evaluating  equation  (1-11). 
The  column  labeled  DELTA  is  the  radius  of  the  field  of  view  in  degrees.  The  column  labeled  UPPER 
BOUND  is  equation  (1-1 1)  integrated  over  the  outer  square  of  figure  12  in  Appendix  A.  The  column  labeled 
LOWER  BOUND  is  the  value  of  the  integral  over  the  inner  square  of  the  same  figure.  The  column  labeled 
ASMPTOTE  is  the  theoretical  value  approached  as  A goes  to  090°.  This  value  is  obtained  by  setting  G ■ 1 
in  equation  (S3)  of  Appendix  A.  A plot  is  also  generated  for  equation  (M);  a sample  is  shown  in 
figure  1-2. 

Figure  1-3  represents  a program  listing  on  19  pages  of  the  downlink  model. 
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Figure  1-3.  ProgrMB  fating  for  dowafink  modd  (sheet  1 of  19). 
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Figure  13.  Program  biting  for  downlink  modd  (sheet  3 of  19). 
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Rpw-U  Program  fating  f«  downfek  model  (sfacct  6 of  19). 
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figure  1-3.  Progm&stiaefwdomnaknodei(dxct7of  19). 


Fipirc  !-3.  Piogi jui  Kstjog  for  dowafak  mode!  (&eet  8 of  19). 
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Fieri  re  |-3.  Program  listing  fot  downlmk  model  (sheet  10  of  19 J. 


Fjjbjc  13.  f^c*r«Ssli^fafdo*«6okKSdd  (sheet  11  of  19), 
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1-3.  Program  feting  fw  dowafcnk  Model  (sheet  15  of  19). 
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Figure  1-3.  Program  tisimg  for  downbak  model  (sheet  16  of  19). 
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fipnI-3.  Program  fistiag  for  doamGnk  model  (sheet  17  of  19). 
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APPENDIX  J 


GENERAL  DESCRIPTION  OF  UPMODEL 


1 


Is 


Program  UPMODEL  is  designed  to  Integrate  equation  (63)  of  Appendix  A,  The  function  evaluated 


1(7')-  //{<!  -e"«z'-A)f(7'(7)  + Af(7'.7)|[  dx  dy  . 


(J-1) 


A modified  form  of  (M)  is  also  computed.  The  modification  is  the  substitution  of  f*(7\y)  for  fty',  7), 
where  f*  is  defined  by 


f*(7\  7)  “ f(7*.  7)exp(-az'sec^) , 


where 


' n_.  n n_  » ,^m,V 


1/2 


(J-2) 


(J-3) 


INPUT  DATA.  All  input  to  program  UPMODEL  is  from  cards  using  NAMELIST  name  INPUr. 

The  variables  to  be  Input  are  as  follows; 

ZL,  A.  S.  THBRSQ,  ^AYERS,  NXGD.  NYDG,  FOV1EW,  GAMP1 , GAMP2,  STEP1 . THZERO, 

GPY  - these  are  identical  to  the  input  for  DNMODEL;  see  the  discussion  there  for  definitions. 

ADATE  - an  Integer  urray  of  dimension  3 for  the  date  in  character  string  form.  Example, 

ADATE  - .’.2  JUN  1975.  Default  value  is  ail  blanks. 

AT1ME  - an  Integer  array  of  dimension  2 for  the  time  of  day  in  character  string  form.  The 
default  value  Is  "bbbbbPDT'. 

OUTPUT  DATA.  A sample  of  the  printer  output  is  shown  following  the  program  listing.  Page  l is 
a list  of  all  input  data.  Page  2 is  a list  of  parameters  for  the  first  set  of  input  data.  Page  3 is  the  tubular  list- 
ing of  the  results  of  evaluating  the  integrals.  The  columns  labeled  GPX,  GPY,  ZENITH,  and  AZIMUTH  are 
the  same  as  in  the  output  for  program  DNMODEL,  The  column  labeled  SUM1  is  the  value  of  equation  (J-1) 


I 


J-1 


A 


directly.  The  column  labeled  SUM2  is  the  value  of  equation  (J-l)  with  f*  substituted  fnr  f.  The  column 
labeled  SUM  is  the  term 


ff  1 -e*,z'-A)f(y,7)dxdy  . (J-4) 

R 

The  column  labeled  SUME  is  the  term 


A KV,  7) 


dx  dy . 


eff 


(J-5) 


The  columns  labeled  SUMS  and  SUMSE  are  the  terms  (J-4)  and  (J-5)  with  f*  substituted  for  f. 

A sample  plot  generated  by  UPMODEL  is  shown  in  figures  J-l  A and  IB.  The  first  curve  is  a plot  of 
SUM1  vs.  zenith  angle  and  the  second  curve  ia  the  plot  of  SUM2  vs.  zenith  angle. 

EXTERNAL  SUBROUTINES  REQUIRED.  COPY,  HEDING.  AXISM,  PHLINE  - See  discussion  of 
these  in  DNMODEL. 

FIT  - This  is  a subroutine  to  locate  the  peak  of  the  data  and  the  points  at  which  the  data  falls  to  a 
value  of  e“*'2  of  the  peak. 

Figure  J-2  represents  a program  listing  on  12  pages  of  the  uplink  model. 
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RpueJ-2.  Program  iang  for  apfak  nodd  (Aeet  2 of  12). 


Figure  J-2.  Program  fating  ter  uptink  mode}  (Aeet  4 of  12). 


Figure  J-2-  Program  fating  for  upfink  modd  (sheet  6 of  12). 


Program  listing  for  uptink  model  (dtect 


FjgnieJ-2.  Program  fating  for  uplink  model  (sheet  8 of  12). 
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Pjpm.g-1-  Prop*  for  ntgaace  sewn  caBtatiaa  (Sheet  3 of  39). 
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Figure  X-l.  Program  for  radiance  scanner  caabratjon  (Sheet  6 of  39). 


Program  for  radiance  scanner  calibration  (Sheet  7 of  39). 
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Figure  K-I  Program  For  radiance  scanner  c •J'-ra.  m (Ebeet  8 of  39). 


Figure  K-l.  Program  for  radiance  scanner  calibration  (Sheet  10  of  39). 


Figure  K l.  Propam  for  aSmce mudh aUbratian  (Sbeet  1 1 of  39). 


F^stK-I.  f’ropaci  for  r?  dance  snancr  dirfxatton  (ShceJ  13  of  39). 
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RpueK-I.  Program  for  iadmicc  vaaner  caBbntioo  (Sheet  14  of  39). 
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APPENDIX  L 

PLOTS  OF  1RRAD1ANCE  LOSS  AT  GIVEN  DEPTH 
AS  FUNCTION  OF  THE  SUN'S  ZENITH  ANGLE 


'.■a 

! 

.1 

.1 


10  FEET 


100 


—k.ifW 


tf"  3 


50  FEET 


1 


QP'sscnaaivyDaiNi 


figure  L-5.  Iiwfant  losw-acft 


Seel 


a § s I 

h Is 


« « a a 

HM 


HU 


it  ii  ii  ii 


© 


oioio: 


© 


® 


© 


ZENITH  ANGLE 


130  FEET 


140  FEET 


APPENDIX  M 


F(ff)  RECEIVER  DATA 


The  F(0)  receiver  was  deployed  to  make  an  independent  meaiurement  of  the  radiant  impulse 
response  as  described  in  Appendix  A,  equation  (7).  Two  sets  of  these  measurements  were  taken  and  are 
plotted  in  figures  M*i  and  M*2,  The  conditions  for  these  measurements  were  as  follows:  the  User  was  in* 
dined  at  an  angle  12.5°  off  the  zenith.  The  receiver,  mounted  alongside  the  barge,  was  moved  roughly  per* 
pendicular  to  the  direction  of  propagation.  At  each  point  (measured  in  feet),  the  receiver  was  rotated  and  a 
radiant  pattern  recorded.  In  figure  M*1 , the  position  of  the  receiver  was  measured  from  the  point  on  the 
surface  where  the  undistorted  beam  would  normally  intercept,  thus  recording  plus/minus  distances.  In 
figure  M»2,  the  distances  were  measured  from  the  location  of  the  laser  and  therefore  only  positive  distances 
occur. 
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Figure  M-2B.  Radiance  profile  through  F(8)  receiver  (July  1975). 
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